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Abstract We review the possible mechanisms for the generation of cosmo- 
logical magnetic fields, discuss their evolution in an expanding Universe filled 
with the cosmic plasma and provide a critical review of the literature on the 
subject. We put special emphasis on the prospects for observational tests of 
the proposed cosmological magnetogenesis scenarios using radio and gamma- 
ray astronomy and ultra high energy cosmic rays. We argue that primordial 
magnetic fields are observationally testable. They lead to magnetic fields in 
the intergalactic medium with magnetic field strength and correlation length 
in a well defined range. 

We also state the unsolved questions in this fascinating open problem of 
cosmology and propose future observations to address them. 

Keywords Cosmology • Magnetic fields • Early Universe • Cosmic Microwave 
Background ■ Gamma Rays 



1 Introduction 



Magnetic fields are ubiquitous in the Universe. Wherever we have the means of 

obser ving them, they are pr esent: in our solar s ystem, in stars (IDonati and Landstreet 
2009), in the Milky Way (jWielebinskil 120051) in other low (jKronber d Il99 i . 
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FletcheiJlioilt lBeckll2012h and hi gh redshift dKronberg et allll992l:lBernet et all 
2008) galaxies, in galaxy clusters (IClarke et a ll200lt Bonafede et a ll2010HFerettr et all 
20121 ). in superclu sters (IXu et a ill2006h and even in voids of the Large Scale 
Structure (LSS) (Neronov and Vovk 2010 : Dolag et all 2011 ; iTavecchio et a" ' 



Tavecchio et alll2010bl : IVovk et aj|2012b iTavlor et alll201ll: iDermer et all 
20111 ). Interestingly, the magnetic field strength in galaxies is typ ically of the 
order of a few-to-tens of y xGauss independent of the galaxy redshift ( Kronberg et a, 
1992t iBernet et all 2008 ) . Also the magnetic fields in clusters are of t he order 
of ^Gauss (jClarke et aj|200lt iBonafede et~all2010l: iFeretti et alll2012h . 

According to a well accepted paradigm, magnetic fields in astronomical 
structures of different sizes, from stars (sizes R ~ 10 11 cm) up to galaxy clus- 
ters (R ~ 10 24 cm) are produced by ampli fication of pre-existing weaker mag 



netic fields via different types of dynamo ( Parker 1955 : Ruzmaikin et all 1988t 



Kulsrudll 19991 : Brandenburg and Subramanianl2005l ; Kulsrud and Zweibell 2008 ) 
and via flux-conserving compression during gravitational collapse accompany- 
ing structure formation. On short distance scales, magnetic fields dissipate 
their energy into turbulent and thermal motions of astrophysical plasmas, so 
that a continuous re-generation of the field is needed on the time scales shorter 
than the life time of the astronomical object carrying the field. This is the case 
for e.g. the magnetic fields of the Earth and Sun and other stars and planets. 
This is also partially true for the galactic magnetic fields, including the field of 
our own Milky Way galaxy. Weak magnetic fields on the largest distance scales, 
from 10~ 2 to 1 Mpc, from the large scale fields in the galaxies to those in galaxy 
clusters, might not have enough time to dissipate their energy into plasma mo- 
tions. Once amplified by dynamo and compression mechanisms, they conserve 
their strength on time scales comparable to the age of the Universe. 

The dynamo and compression amplification mechanisms can act only if 
a non-zero magnetic field is present. This "seed" field for the amplification 
might be tiny, but it has to be generated by a different mechanism, which 
pre-dates the structure formation epoch or operates at the onset of struc- 
ture formation. The uncertainty of the strength and of the origin of this ini- 
tial seed field c onstitutes the long-standing problem of the origin of cosmic 
magnetic fields (Kronbergl|l994 : Grasso and Rubinstein! 2001 : Widrow 2002 : 
Kulsrud and ZweibefcOOa iKandus et alll201lH Widrow et alll2012l ). Two broad 
classes of models for the origin of the seed fields are discussed. One possibility 
is that the weak seed fields are produced in the early universe, during epochs 
preceding the structure formation. Another possibility is that the process of 
generation of the seed fields accompanies the gravitational collapse leading to 
structure formation. 

The existing data on magnetic fields in galaxies and galaxy clusters cannot 
provide direct constraints on the properties and origin of the seed fields. This is 
related to uncertainties of the details of the dynamo mechanisms operating in 
galaxies and clusters on the one hand and, on the other hand, to the numerous 
saturation effects which drive the galactic and cluster field strengths to fixed 
values largely independent of the properties of the initial seed fields. 



Cosmological Magnetic Fields: Their Generation, Evolution and Observation 



3 



The only potential opportunity for understanding the nature of the initial 
seed fields is to search for places in the Universe where these fields might exist 
in their original form, not distorted by the complicated plasma and magneto- 
hydrodynamics (MHD) processes. 

The only places where such "primordial" magnetic fields might reside is 
the intcrgalactic medium (IGM) , more precisely, the voids of large scale struc- 
ture (LSS). If weak magnetic fields were indeed present in the Universe before 
the onset of structure formation, they did not suffer much amplification be- 
cause of the absence of the dynamo and compression of the IGM in the voids. 
Cosmologically produced magnetic fields might passively evolve (be diluted by 
the expansion of the Universe) still today. Potential measurements of Inter- 
galactic Magnetic Fields (IGMF) using available observational techniques of 
radio, microwave and 7-ray astronomy might, therefore, provide an important 
clue on the origin of the seed fields. This idea is the prime motivation for the 
numerous efforts to detect the IGMF. 

If successful, detection and measurement of the properties of primordial 
magnetic fields in the voids of LSS will provide an extremely important source 
of cosmological data. Typical scenarios for generation of magnetic fields in the 
early universe concentrate on possibilities of field production via charge sepa- 
ration and/or generation of vortical currents at the moments of cosmological 
phase transitions: the electroweak and the QCD phase transitions, and the 
moments of photon decoupling and recombination. Another possibility is the 
quantum generation of very long wavelength photons during inflation which 
then are converted into magnetic fields at reheating. In most of the models the 
moment of cosmological " magnctogenesis" pre-dates the epochs of formation 
of Cosmic Microwave Background (CMB) signal and the Big Bang Nucleosyn- 
thesis (BBN). This means that the detection of relic magnetic fields might 
provide observational data on very early physical processes in the hot Uni- 
verse with temperature in the range above 100 MeV. If they stem from the 
electroweak phase transition or from inflation, they may even probe physics 
beyond the standard model. 

It is not obvious a-priori that weak magnetic fields which reside in the 
voids of LSS are primordial. Alternatively, they could be produced at the late 
stages of evolution of the Universe (at redshifts z < 10) by outflows from al- 
ready formed galaxies. These outflows can b e galactic winds generated by the 
star formation activity ( Bertone et alll2006h and/or relativistic outflows gen- 
erated by the Acti ve Galactic Nuclei (AGN) (|Reeslll987l:lDalv and Loeblll990l 
Ensslin et all 1997 ) . Both types of outflows are essential elements of the struc- 



ture formation process. They are responsible for washing o ut the baryon con- 
tent of galaxies leading to the "missing bary ons" problem ( Cen and Ostriker 
1999) and for metal enrichments of the IGM (jAguirre et afcOOll) . If these out- 



flows are (a) strongly magnetized and (b) able to spread into the voids of the 
LSS, they can result in "pollution" of the voids with magnetic fields which are 
much stronger than the relic magnetic fields of primordial origin. The pres- 
ence of magnetic fields spread by galactic wind into the IGM may prevent 
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measurement of the relic cosmological fields. At the same time, detection and 
measurement of IGMF spread by winds can constrain the properties of galactic 
winds and in this way shed light into the physics of the "feedback" provided 
by the winds on the formation and evolution of galaxies. Winds are thought 
to be responsible for the regulati on of star formation activit y of galaxies at 
different stages of their evolution ( Kennicutt and Evans! 2012 ). 



In this review we summarize and critically asses the current knowledge 
of the weakest magnetic fields in the Universe. Besides describing the present 
status of observations of magnetic fields in the IGM, there are two fundamental 
questions which arise and which we shall address: 

1 . Are IGMF primordial, in the sense that they have been present before the 
galaxy formation process took place, with all its complicated non-linear 
and non-gravitational physics, or have they been formed during galaxy 
formation e.g. in star formation and AGN activity and then spilled out 
into the galaxy and into intergalactic space? 

2. If they are generated in the early Universe how do they evolve? Are they 
just decaying with redshift z like l/(z + l) 2 as flux conservation would 
demand and maybe damped on small scales by diffusion (what are these 
'small scales') or can MHD processes move power from small to larger 
scales? 

To address the first question, we review the mechanisms to generate primor- 
dial magnetic fields. They fall into three broad classes: inflationary magnetic 
field generation, the generation of magnetic fields during phase transitions and 
magnetic fields from second order cosmological perturbation theory. 

Once the fields are produced at a certain " magnetogencsis" epoch in the 
early universe, they evolve interacting with different types of plasma of charged 
particles. Such plasma is present in the Universe both before the moment of 
recombination and after the epoch of reionization. A judgement of the primor- 
dial nature of IGMF is not possible without an understanding of the evolution 
of the field from the event of production until the present. 

To address the second question, we review current understanding of the 
evolution processes and describe evolutionary tracks of magnetic fields in the 
(B,Xb) parameter plane, where B and Xb are the most important integral 
(i.e. distance scale-averaged) characteristics of magnetic field: its strength and 
correlation length. 

Finally, we review present observational constraints on the IGMF and dis- 
cuss possible ways how to distinguish whether the IGMF comes from primor- 
dial fields or from fields produced by galactic outflows. 

Our review is organized in the following way. In Section [2] we summarize 
the basic equations governing the evolution of magnetic fields and charged 
particle plasma in the Universe. Next, in Section [3] we review previously pro- 
posed mechanisms of generation of magnetic fields in the early universe. In 
Section 2] we discuss the evolution of primordial magnetic fields. We discuss 
different damping and amplification mechanisms, determine their character- 
istic scales as function of cosmic time and we study their effects on different 
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types of primordial magnetic fields. We also investigate the possibility of an in- 
verse cascade. In Section [5] we review the observational situation and compare 
different constraints on IGMF. We also discuss future capabilities of several 
different observational strategies. In Section [5] we re-discuss different cosmo- 
logical magnetogenesis models to single out models which can be tested by 
observations. We show that relic cosmological fields occupy a distinct region 
in the (B, Xb) parameter plane. This region is different from that expected for 
magnetic fields from galaxy outflows. We argue that this opens a possibility 
to distinguish between the primordial and galactic outflow produced fields ob- 
servationally. In each section we also point out the main open problems which 
still have to be addressed in future research. In Section [7] we conclude. 



2 Basics of magnetic fields in an expanding Universe 

2.1 Notations and definitions 

In this review we consider magnetic fields evolving in a flat expanding universe 
described by the metric 

ds 2 = a 2 {t) [-dt 2 + 6 tJ dx i dx 3 ] , (1) 

where a is the scale factor which we normalize to unity today such that it is 
related to the cosmological redshift by 1 + z = 1/a. Here, t denotes conformal 
time which is related to physical time r by adt = dr. The Hubble parameter 
is 

da/dr a „ , . 
H=—t — = — =n/a. (2) 
a a z 

An overdot denotes derivative wrt. conformal time t, and T-L is the conformal 
Hubble parameter. This focus on conformal variables is useful since electro- 
magnetism is conformally invariant so that, as we shall see, with suitable 
re-scaling, the equations for a magnetic field in expanding space are identical 
to those in Minkowski space. 

We denote spacetime indices by greek letters and 3d spatial indices by latin 
letters, 3d vectors are denoted in boldface. We use natural units in which the 
Planck constant and the speed of light are unity: K = c = 1. The reduced 
Planck mass is denoted by Mp, such that 8ttG = 1/Mp. We also set k\, = 1 so 
that temperature is a measure of energy. More precisely, IK = 0.86 x 10 _4 eV. 

The electric and magnetic field arc defined in the reference frame co- 
moving with the coordinate system in which the space-time metric has the 
form (TT]), with the time a xis directed a l ong th e vectors of the Hubble flow 
(u^) = o _1 (l, 0,0,0). As in lBarrow et~a] (|2007t ). we define 



= F^u v , Bp = , (3) 

such that 

F^ v = u^Ey - u v E^ L + e^B 1 . (4) 
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Here e M „ 7 = u a e a ^ v ^ is the totally antisymmetric tensor on the 3-space normal 
to , while e a ^ UJ is the totally antisymmetric tensor in 4 dimensions with 
e oi23 = y/— 9- The indices of t^ vl are are raised and lowered with h a p = 
gap + u a up. Note that both B a and E a are normal to the four-velocity . 

These definitions are useful also for a generic 4-velocity which need 
not even by hypcrsurface-orthogonal. However, since in a perturbed Fried- 
mann Universe, both, the deviation of matter / plasma velocities m m from the 
Hubble flow and the electromagnetic field are small, it suffices to consider 
the background velocity in the definition of and B^, so that the elec- 
tric and magnetic fields are parameterized by the 3d vectors {E^) = a(0,E), 
(B/i) = a(0,B). Maxwell's equations in terms of E and B are given in Ap- 
pendix |Bl eqs. (|184l) to (|187[) . They are simply re-scaled versions of the equa- 
tions in Minkowski space. Here and in the following we denote 3d vectors in 
boldface. 

Homogeneity and isotropy of the Universe implies that statistically the 
spatial structure of magnetic fields is the same at any location in the Universe. 
In the view of this, it is often convenient to study the properties of magnetic 
field in terms of its Fourier components, 

B(k,t) = J AB(x,t)e !k - x . (5) 

Statistical homogeneity and isotropy then imply that expectation values can- 
not depend on any vector except k and on any tensor except Sij and Cy- TO as 
well as combinations of these. The spectrum of the magnetic field therefore is 
of the form 

a 4 (A(M)5*(k',i)) = (27r) 3 <5(k - k') [(<% - fe^uC*) - ie ijm k m P aB {k)\ , 

( 6 ) 

where k = k/fc and k = |k|. The bracket ( ) denotes an ensemble average, i.e. 
an average over many realizations of the stochastic magnetic field. In observa- 
tions we of course always only measure one realization, but it is usually justficd 
to assume an 'crgodic hypothesis', namely that spatial average over many in- 
dependent patches of size L 3> 2n/k is a good approximation to the ensemble 
average, especially, if the size L is larger than the cosmological horizon at the 
time when the magnetic field was generated. 

Pb and P a B are the symmetric and antisymmetric part of the magnetic 
field spectrum. The (5-function is a consequence of spatial homogeneity and the 
tensor structure comes from the homogeneous Maxwell equation V • B = 0. 
In terms of a right handed orthonormal system (e*- 1 ', e^ 2 \ k) with = 
(e« +ie^)/2 and e<~> = (e« - ie^)/2, we have 

B(k,f) = B {+) (k,t)e (+) + .B(_)(jfe,i)e (-) and (7) 
(B (+) (k)S ( * +) (k')) + (B ( _) (k)S ( *_ 5 (k')) = (2tt) 3 * (k - k')P B (fc)/a 4 (8) 
<2? (+) (k)£ ( * +) (k')> - (B ( _ ) (k,f)B ( *_ ) (k')) = (2^(k-k')P QB (fc)/a 4 .(9) 
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Again, the angular brackets signify "ansemble averaging". Eqs. © and ^ 
are obtained easily with the help of the identities (e( + ') 2 = (e^^) 2 = and 
e (+) . e (-) = 1 together with k A e (+ ' = —ie^ and k A e ( ~) = ie^K Parity 
transforms into B^~> and vice versa, hence Pb is even and P a B is odd 

under parity. Eqs. ^ and © also imply that Pb > |P b|- Equality is reached 
if one of the hclicity modes vanishes completely; such a field is called totally 
helical. 

Th e energy density of the magnetic field (in Hcavyside-Lorentz units see lJackson 
is given by 



a 4 PB =a 2 J dkp B (k) = J ^k 3 P B (k) . (10) 

Here ps is the ensemble average of the magnetic field energy density, which is 
independent of position. Eq. (|10[) is obtained by noting that 



p B = ^(B(x)B(x)) = -i-^ J d 3 fcd 3 fc'(B(M)B*(k',t))cxp( 4 x. (k-k')) 
a~ 4 f ,^ , „ , , > a A f dk , q „ . 

= (^J dkPB{k) = ^jT kP ^- 

Note that B(x) 2 has the dimension of energy density, hence B(k)B(k') has the 
dimension of energy density x length 6 and P(k) has the dimension of energy 
density x length 3 which is the dimension of the Fourier transform of the energy 
density as required. Hence PB{k) = dps/dk = k 2 Ps(k) / '(2n 2 ) is the energy 
density per unit k interval. We shall sometimes also employ the energy density 
"per log interval": dp B /d\og(k) = k 3 P B (k) / '(2tt 2 ). 

It is also convenient to introduce the " characteristic" magnetic field strength 
at scale A = 2w/k, 



dlog(fc) 



and B = ^2p B , (11) 

fc=27r/A 

where log denotes the natural logarithm. We shall systematically use the field 
strength on scale A and the scale-averaged field strength B in the following 
sections. 

Below we shall also need the power spectrum of the fluid velocity field, we 
therefore also introduce it here. 

v(k,f) =y e ?;kx v(x,t)d 3 x (12) 

and 

(vi(\s.,t)vj{k',t)) = (27r) 3 <5(k - k') 5ijP sK - hk P vK - ie ijm k m P aK (k) 

(13) 

Like the index B for "magnetic", the index K stands for "kinetic". Note that 
for P s k = P V K the velocity field is divergence free. The fact that (vjV*) > 
implies P s k > PvK- 
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Similarly to the energy of magnetic field, given by Eq. f| 10[) . we also intro- 
duce the spectral kinetic energy 

PK = i / ^ ( 3P ^( fc) _ P -K(k)) dk = J p K (k)dk , (14) 

where px (fc) is kinetic energy per unit wave number interval. Note that, con- 
trary to pb, Pk is dimensionless and has to be multiplied by the plasma energy 
density p to yield the true kinetic energy. We shall however use the customary 
language and refer to pk as kinetic energy density. 

We assign a characteristic correlation length to the stochastic magnetic 
and velocity fields, defined by 

A B = 2-Kpg 1 J p B {k)k~ 1 dk , 

\ K = 2np K 1 J pui^k^dk. (15) 

This is the characteristic scale, also called 'integral scale' of the magnetic field 
and the velocity field. We shall sometimes also call it the 'correlation scale' 
even though this is not strictly correct in a statistical sens^E 

The magnetic field and the velocity field are said to be in equipartition 
if pb/p — ( v a)/2 = Pk — ( w2 )/2- Here va is the Alfvcn speed defined by 
v\ = B 2 /(2p). If Ps(k) ~ a 4 pPftr/2 we speak of detailed equipartition or 
equipartition on all scales. 



2.2 Helicity 

The magnetic helicity is the volume integral 

H(V) = I A ■ Bdv , (16) 
Jv 

over a volume through the boundary of which no magnetic field lines cross. The 
3d vector A is the magnetic vector potential. The above volume can also be 
infinite if the magnetic field decays sufficiently rapidly at infinity. The helicity 
is gauge independent, since under a gauge transformation, A — > A + Va 

H(V) -> H(V) + I Va ■ Bdv =) -> H(V) + <f aB ■ nds = H(V) . (17) 

JV JdV 

Here n is the normal to the boundary dV and we have assumed that B • n = 0. 
Magnetic helicity has a simple topological interpretation in terms of linking 

1 In statistical mechanics correlations decay exponentially on scales larger than the cor- 
relation scale while our correlations usually decay like a power law. Therefore, even though, 
most of the magnetic/kinetic field energy is concentrated on scales close to As respectively 
\k, this is not true for all its cumulants. 
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and twist of isolated flux tubes, see Brandenburg and Subramanianl ( 2005 ) 



The hclicity density of the magnetic field is then given by h = A ■ B. 

Using a gauge in which A is transverse, k • A = 0, we have in Fourier space 
kA = ik A B which yields 

/dk dh 1 f dk , o „ , . 

Td^) = ^jT kP ^). (is) 

Hence fc 3 P a s(fc)/(27r 2 ) is the helicity density per log-k interval. 

In a Friedmann-Lemaitre universe, the electromagnetic Lagrangian, 

L = Iv=S V" = \if a V uP F^F aP 

is independent of the scale factor a(t). Hence a freely propagating electromag- 
netic field F^ a is independent of a. This is simply a manifestation of conformal 
invariance of electromagnetism in 4-dimensions. This implies that Bi oc 1/a 
and B l oc 1/a 3 such that B 2 oc a~ 4 . We have taken out this trivial conformal 
scaling in the power spectra in Eqs. (|6l8p and ©. 

Below we see that this scaling remains true when interactions with the 
cosmic plasma are relevant in the special but cosmologically most relevant 
interacting case of the magnetic hydrodynamic (MHD) limit due to flux con- 
servation. 

If we write B = .B l e, for the orthonormal basis = a -1 ^, the scaling of 
the components is B l oc a~ 2 . 



2.3 Co-evolution of the magnetic field and the cosmic plasma 

Dynamical equations for the evolution of the interacting matter and electro- 
magnetic fields in the expanding Universe are derived starting from the law of 
conservation of stress-energy tensor, = 0, for the stress-energy tensor con- 
sisting of the electromagnetic and plasma (fluid) contributions. In the simplest 
case of an ideal fluid, its stress-energy tensor is 

Tjf = (p + V )u»u v - pg» v (19) 

where p and p are the energy density and the pressure of the fluid and w M 
is its four-velocity. We assume that the different components (photons, elec- 
tron/positrons gas etc.) of the dominant rclativistic particles are sufficiently 
strongly coupled so that we can consider them as one fluid. (For T ~ lMcV 
this means that we neglect the neutrinos in our qualitative considerations). 
The stress-energy tensor of electromagnetic field is, see Appendix IB1 

T^ m) = F^F\ - X - gilv F x °F Xa (20) 
= ^{E 2 +B 2 )u^u v + ^(E 2 +B 2 )h^-E^E v -B^B v + P^u v +u^P v . (21) 
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where is the Poynting vector. Following Brandenburg et al (1996), we in- 
troduce the following conveniently rescalcd quantities, 



~ 4 ~ 4 7-1? 

p = a p, p = a p, B 



a 2 B\ E i = a 2 E\ J 1 = a 3 J\ 



(22) 



As the present value of the scale factor is unity, this implies that the tilde- 
quantities correspond to their values scaled to today. In the rest of the review, 
when ever we indicate B or p we mean the value of the magnetic field or of 
the energy density scaled to today. 

Using Maxwell's equation, see Appendix [Bj we can write the conservation 
equations, (Tp U + T^ M ). V = in the form 



W/3)] + |v(jiv) = -J-E 



4 / dv dp _. _ ,_ , 

3 ^+v ¥ +p(vV)v + ,v ( V.v, 



-Vp + J A B . 



(23) 
(24) 



The derivatives are wrt. conformal time t and comoving coordinates x. Here we 
have neglected terms which are of third order in the perturbed quantities like v, 
B, dtp etc. Even though we consider a relativistic fluid with p = p/3, peculiar 
(bulk) velocities are small. Nevertheless, we want to keep quadratic terms in 
order to be able to describe non-linearities which can provoke modifications in 
the spectrum like an inverse cascade. 

In the early Unive rse, conductiv i ty is v ery high, for relativ i stic electrons 
we typ ically have, see lEnqvist et a] (jl995t) ; lArnold et all (|2000h : lArnold et all 
(|2003l) and Appendix |A] 

° - ] 7 _ n • (25) 
a log(a 1 ) 

where a is the fine structure constant. It therefore makes sense to work in the 
ideal MHD limit where 



E = -v A B and J = V A B 



(26) 



to lowest order. The first equation is simply the condition that the Lorentz 
force on charged particles vanish. The second equation follows from Ampere's 
law using E <C B, see Appendix [Bj In this limit E is already of quadratic 
order and we can consistently neglect the 3rd order term J • E in Eq. (|23)l . In 
this approximation, up to first order dtp = — (4/3)pV • v. 

However, viscosity can become significant and we want to take it into 
account. We also take into account the damping of the magnetic field due 
to Ohmic losses. Rcscaling also shear viscosity, v = u/a and the conductivity 
a = aa, including dissipation in Eqs. (|23I24[) they become, see lBaneried (|2002l) 



dp 
dt 



+ V ((p + p)v) = , 



(27) 



^ + (vV)v 



dp 



Vp B A (V A B) 



(p + p)dt (p + p) 



(P + p) 
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i>(v 2 v+iv(V-v)j , (28) 

^--VA(vAB) = -V 2 B. (29) 
ot a 

In addition , we consider a r adiation dominated equation of state p = p/3. With 



respect to iBaneried (|2002l ) we have neglected 'heat losses' as in the radiation 
dominated era the photons are part of the plasma and their energy density is 
included in p. 

The first equation is just the continuity equation for the cosmic plasma. 
The second equation is the Euler equation with the dissipation term on the 
right hand side. The rescaled shear viscosity is of the order of the comoving 
mean free path of the plasma, v ~ A m f p /5, see Appendix El 

On small scales A <C A m f p /5, we have to replace diffusion damping by 
damping due to free streaming. This can be done by replacing the dissipation 
term by — av, where a oc A~L (see Baneriee and Jedamzikl ( 20041 )). A more 



rigorous treatment would require to solve the Boltzmann equation. Neverthe- 
less, qualitatively we expect fluctuations to be damped by diffusion and by 
free streaming on scales A ~ A m f p . 

The third equation is the magnetic induction equation, with the Ohmic 
dissipation term on the right hand side. 

If the dissipation terms are subdominant to the non-linear terms, MHD 
turbulence develops. This is controlled by the Reynolds numbers, 

R-fc(fc) = 7-7 (kinetic Reynolds number), (30) 



kv 

Rm(fc) = — — (magnetic Reynolds number), (31) 
k 

R 

P,„ = — — = av (Prandl number). (32) 

Rfc 



Here k is some comoving wave number and Vk = \J (|v| 2 ) = \j Pfcifyk 3 /(2ir 2 ). 
In Appendix lAl we compute these numbers for k = ks = 2tt/Xb rsp. k = kx 
as functions of the temperature and show that for T < 100 GeV the Prandl 
number is much larger than one, so that we may neglect magnetic diffusion 
with respect to the kinetic one which is much faster. 

On scales where the Reynolds numbers are large, the quadratic terms (v • 
V)v ~ kv 2 and V A (v A B) ~ kvB dominate over the damping terms ~ vk 2 v 
(or k 2 B/a for small Prandl number) and turbulence develops. 

On small scales, k > kd — v^/v, the damping term dominates and the ve- 
locity field is damped. Due to the coupling to the magnetic field, this damping- 
is only a power law, but nevertheless very rapid (on the timescale td ~ k^ 1 . 
Once the velocity field is essentially damped away, the quadratic term in the 
induction equation drops and the magnetic field remains frozen. Later on, 
when the viscosity scale becomes smaller, the magnetic field re-generates a 
velocity field. 
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If (which is not the the case in the situations we are interested in) k > kd — 
va > v/v, the magnetic field is damped and non-magnetized fluid turbulence 
remains. 

If the fluid is incompressible, dtp = 0, Eq. ([27)) implies V ■ v = and 
the flu i d motion is purely vortical. It h as been argued, see e.g. iJedamzik et al 
(|l998l ): iBaneriee and Jedamzikl (|2004[ ). that this is the case for cosmological 
magnetic fields. However, even though we know that cosmological density fluc- 
tuations are small on large scales, this need not be the case on small scales. 
Especially, since the energy density p also contains kinetic energy, we expect 
its fluctuations to be at least of order v 2 . This means that compressible terms 
in the MHD equations are, in general, as important as the incompressible ones. 
We shall see in Section 01 that this is relevant for the evolution of the magnetic 
fiel d spectrum. 

( Bovarskv et al 2012af ) have shown that an additional effective degree of 
freedom (chiral asymmetry) should be added to these equations. Its origin is 
a subtle quantum effect - the chiral anomaly - that couples the change in the 
number of left and right-chiral particles with the change of the helicity of the 
magnetic field. Taking into account this degree of freedom and its interaction 
with electromagnetic fields significantly changes the evolution in the case of 
the strong helical magnetic fields at temperatures above a few McV. 

Let us summarize the situation as follows: From equation (|28[) we see that 
magnetic field sources the velocity field. Thus, the process which leads to the 
production of magnetic fields simultaneously sets the plasma in motion. The 
nonlinear form of the Eulcr equation assures that the plasma motions are 
turbulent, so that the process of generation of magnetic fields is inevitably ac- 
companied by the excitation of plasma turbulence. Furthermore, the turbulent 
velocity field couples back to the magnetic field via the term V A (vAB) which 
leads to turbulence also in the magnetic field. Thus, a consistent description 
of the co-evolution of magnetic field and plasma in the radiation dominated 
Universe has to be described in the "language" of MHD turbulence (see e.g. 
the books bv lBiskamd (|2003l) ; iTsvtovichl (|l977h V 



3 Generation of primordial magnetic fields 

3.1 Inflationary magnetic field production 

The electromagnetic field is conformally coupled and does not 'feel' the ex- 
pansion of the Universe. Therefore, in order to generate magnetic fields during 
inflation, one has either to couple the electromagnetic field to the inflaton or 
to introduce another coupling which breaks conformal invariance, e.g. a term 
Lint oc R^va/jF^ F a ^ or even break gauge invariance, li ke R l ,, v A fl A L ' for ex- 



ample . These possib i lities have first been investigated bv iTurner and Widrowl 

( 19881) and bvlRatral d 19921) and later been revisited by many authors (|Martin and Yokovamal 



2008 



2012 



Subramanianll2010l; lKunzell2010l ; iKandus et"all201ll; iMotta and Caldwelll 



Jain and Slothll2012l) to cite a few recent accounts. Another possibility 
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is that during inflation gauge symmetry is broken a nd the gauge fields be- 
come massive, which also breaks conformal symmetry (jEnqvist et alll2004f ) . In 
Section [6] we discuss possible observational signatures of inflationary magnetic 
fields. Below we indicate some of the constraints which are summarized in 

Fig.mi 



3.1.1 Standard inflaton coupling 



We consider the Langangian 

L = 



(33) 



Adopting Coulomb gauge Aq( k, t) = 0, djA^{x.,t) = and following the no- 
tation of Subramanian ( 201df ). Maxwell's equations, [pF^]^ = 0, lead to 
an evolution equatio n for the space com ponents Aj(x, t). In a cosmological 
background it reads ( Subramanian! 201dt ) 



At 



f ■ 

l Ai 



AAj = 



(34) 



where A is the comoving spatial Laplacian. For a Fourier mode k, we simply 
have A = —k 2 . The time evolution of the vector potential depends on the 
coupling function f(ip). One may adopt, at least for a sh ort time, a simple 
power law in conformal time (jMartin and Yokovamall2008l ): 



/(*) = h 



(35) 



For example for power law inflation with an exponential potential this cor- 
responds to a coupling of the form / cx exp(—a<j)/M). For this coupling the 
damping term is simply 2j/t 2 and in Fourier space eq. (|34p can be solved in 
terms of Bessel functions. Setting for the electromagnetic potential in Fourier 
space AW(k,t) = A(k,t)e ( - ± ' ) /a we obtain 



A(k,t) 



Ci(7)J 7 _ 1/2 (a;) +C 2 {j)J_ y+1/2 (x) 



(36) 



where x = \kt\ = —kt, J v denotes the Bessel function of order v, and C±, C% 
are 7 dependent coefficients which are fixed as usual by imposing v acuum 
initial condition on sub-horizon scales, — kt — > 00 1 Subramanian 201fJ ). Note 
that during inflation conformal time t is negative. For t he symmetric magnetic 
and electric field spectra a nd their correlator we obtain ( Martin and Yokovamal 
20081 : ISubramanianllioioh 



Pp. 



Peb = 4tt 




Pf 



A*{k,t) 



4tt 



A(k,t) 



and 



(37) 



(38) 



14 



Ruth Durrer, Andrii Neronov 



In this case no anti-symmetric part is generated, P a B = 0. 

On super-horizon scales, x <C 1, we can approximate the Besscl functions 
by power laws so that these spectra become 

M ' ] P \(l-2 7 ) 2 |c 2 | 2 x^if 7 >-l/2 1 ' 

„ ^k ["M^ +1 if 7 < -1/2 

PEB(k,t) = —l (2 7 -l)cJc 2 if - l/2< 7< 1/2 

J l(27-l)|c 2 | 2 .T 1 - 2 Tif7>l/2. 

(41) 

The coefficients Ci are 7-dependent but of order unity. We want to discuss the 
dependence of these spectra on 7. First of all, in order to avoid an infrared 
singularity in this simple model we must require — 2 < 7 < 2. At the boundary 
the divergence is logarithmic and can be removed in a way which depends 
only very weekly on the cutoff. Also, when 7 < the magnetic energy density 
dominates while for 7 > the electric energy density dominates. For 7 = — 2 
the magnetic power spectrum is scale invariant and we obtain 

^-- 7 M ! I<P^^/(« 2 t 2 ). (42) 
dlogfc 7T J{ if a 4 

The condition < P4, is required such that we can neglect the effects 

of the magnetic energy density on inflationary expansion (backreaction). The 
same condition has to be satisfied independently by pe- Normalizing e 2 such 
that / = /1 = 1 after inflation, we obtain in all cases the ratio 



Pb 1 / Hh 



p rad {t iai M P ) 2 V m p 



(43) 



With a suitable choice of / 2 1\ it is then easy to obtain magnetic fields of 
the order of e.g. 10 _9 Gauss/a 2 on all cosmologically relevant scales in the 
scale invariant case, 7 ~ —2, while maintaining the condition f\t\a < Mp 1 
during all of inflation, in order to prevent back reaction. After inflation, the 
conductivity of the cosmic plasma is very high and the electric field is rapidly 
damped. The inflaton is frozen and the function f(4>) — > 1. This scenario 
has one serious problem: For 7 ~ — 2 / is a rapidly growing function during 
inflation. On the other hand, the electron field does not couple to A but to 
the canonically normalized electromagnetic potential, \fjA-. The charge of 
the electron is therefore ej\f] rapidly decreasing. To arrive at e 2 = 1/137 at 
the end of inflation, e 2 / / must have been much larger than 1 during most of 
inflation. The electron field becomes strongly coupled and we cannot trust our 
perturbative quantum field theory calculation anymore. This problem has been 
noted first bv lDemozzi et a 1 (2009). We cannot solve it by simply changing 
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to y/jAp since such a coupling violates gauge invariance. This is actually the 
only way to save this model, to violate gauge invariance. The consequences of 
this, e.g. the generation of electron-positron pairs due to this coupling to the 
inflaton have not yet been studied. 

This problem is avoided if 7 > 0, hence / is decreasing. But then, since the 
magnetic field power spectrum is a power law with spectral index n s = 1 + 27 
(if 7 < 1/2) and n s = 3 — 27 (if 7 > 1/2) the magnetic field spectrum is very 
blue. Let us denote by fc ma x the smallest scale on which B is still generated 
during inflation, i.e. the scale that exits the horizon briefly before the end of 
inflation. If on this scale, the magnetic field energy density is a fraction e of 
the radiation density after inflation, on some other scale fci we then have 



l B 2 {k) d PB (k) 



2 v ' dlogfc 



fn 1 (fci/fc max ) 4+27 ,7<l/2 UA , 

k=k 1 l(fcl/femax) ,7>l/2- 



This magnetic field spectrum for the case 7 = is shown in figure [TH] for two 
different values of the inflation scale. 

Knowledge of the initial spectrum of the magnetic field allows to make 
predictions for the expected "relic" magnetic field which might survive until 
the present epoch. We shall in the following term 'naive evolution of the power 
spectrum' and evolution where Ps(k) does not change on large scales and is 
simply damped away beyond a certain small damping scale fcdamp(t) which 
may depend on time. This 'naive evolution' assumes that beyond the damping 
scale, the magnetic field just scales like B ex 1/a 2 which is required by flux 
conservation. We shall show later, in Section |4j that evolution is usually more 
complicated. Assuming naive evolution, Eq. (|44j) is nearly time independent 
since both, p rad and ps scale as a -4 (apart from the changes in the number 
of relativistic degrees of freedom which we do not expect to account for more 
that one or two orders of magnitude). The present radiation density is given 

by 



(3 x 10 _6 G) 

p rad (t ) ~ 2 x l(T 15 (eV) 4 ~ 4.66 x l(r 34 g/cm 3 ~ ^ '- . (45) 

07T 



2 



Inserting this in eq. (j44j) we obtain 
dp B {k) 



d log k 



k—k\ 



= e| *j_\'"(8»10-O) \ m 

tax / 87T 



If inflation happens at high energy with £/i n f ~ E? ni /Mp with E ln { ~ 
10 15 GcV, this yields 

fcmax ^ 1/icnd ~ Y+l~ d Mp~ { lO^GeV j ' 

where t cn d is the comoving itme at the end of Inflation and To is the present 
temperature of CMB, T Q ~ 2.3 x 10" 4 cV. On a scale of say k^ 1 - IMpc 
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3 x 10 24 cm, cq. (|4^)l then yields a tiny left over field of 
dp B {k) 



dlogfc 



fc 1= l MpC 87F V V ^inf 



Here we have set 7 ~ to obtain the most optimistic value, which is still 
devastatingly small. 

Lowering the inflation scale helps somewhat but even when setting it to the 
electroweak scale, E- m { ^200 GeV, we obtain only fields of 10 _23 Gauss on Mpc 
scales, assuming e ~ 1. However, the natural normalization of the magnetic 
field energy spectrum is ps ~ fc* ax ~ H* nl ~, while the energy density of 
the Universe scales as p ~ MpH 2 . This means that typical model calculation 
results in e ~ H? nf /Mp <C 1 for H ln f <C Mp. Thus, in fact, lowering the energy 
scale of inflation gencrically results in weakening of magnetic fields. 

Very generically we shall see that if we want to generate magnetic fields 
early and if we want to have reasonably large fields also on large scales, in 
order for the small scale fields not to over-close the Universe their spectrum 
should not be very blue. We must have either n s ~ —3, or an evolution which 
raises the magnetic field power on large scales by some plasma processes. The 
latter is called an 'inverse cascade'. 

Furthermore, from inflation we expect e = ps / p ra d — Hf nf /Mp. This value 
is indicated by the thick solid line in Fig. If 6I T0 obtain much larger amplitudes, 
like e.g. equipartition, e ~ f as assumed in the dashed line annotated by Ei n { = 
200GeV, we need in addition a dynamo mechanism e.g. during reheating which 
rapidly amplifies the magnetic field to equipartition. This possibility is not 
excluded but also not confirmed by any detailed study. 

3.1.2 Coupling to curvature 

We want to discuss biefly also another possibility, namely that the electromag- 
netic field is coupled to curvature. We consider the Langrangian 



—— + - ( F^ F 4- — B^ va ^ F F a 
2k z 4 V m 



Varying the action with respect to we find 



(48) 



a 
m 

In a Friedmann Universe we have 

FL Fi? 

—8) and R lJ tm = — 
a z J a 



d»[F flv + ^R fl „ af3 F ap >)=0. (49) 



B°\ 3 = ^8} and R ij lm = ^ - ^) with (50) 



FL 2 /a 2 = p/(3M P ) and Ufa 2 = {p+ 3P)/(6M'p) . (51) 

During perfect de Sitter expansion, -R AI „ Q ' 3 is constant and the curvature term 
does not affect the equations of motion. This is also true for the other possible 
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curvature terms, RF^F^ and R jlv F^ la F v a . Therefore, in this case magnetic 
field production is suppressed by the slow roll parameters. 

We now assume p = wp with —1 ~ w. This is the case of power law 
inflation, where the scale factor and the energy density behave like 

2 -6(l + m) 

acxt 1 + 3 "' and poet i+a™ , so that (52) 
(H/a 2 ) = (l + '3w)H 2 /a 2 , (H/a 2 )- = -6(1 + w)H 2 /(a 2 t) . (53) 

Inserting this in the equation of motion for A in Coulomb gauge we find 



2a ,„ „ V-A 'j 12a s „ , 1 V , 9 / , 2a _ 

t m z 



(1 - — (1 + 3w)H 2 )A + —(1 + w)H 2 -A + k 2 (l - -^H 2 )A = . (54) 



Typically, the relevant mass scale is the electron mass, m ~ m e and the 
curvature terms dominate in the early universe when p > Mpm 2 . Terms 
of this form do actually occur in 1 loop vacuum polarization calculations, 



Drummond and Hathrell ( 1980l ). This is the situation we want to consider. 



see 

We therefore neglect the standard term and obtain 

\ _ ^i±4li, + e^—L = . (55) 
(1 + 3w) t (1 + 3w) y ' 

On large scales, \kt\ <^ 1 there is an uninteresting constant mode and a mode 
behaving like 

~ Al oc + ™ OC a (7+9-)/2 _ 

For w > —7/9 this is a growing mode. The general solution is again given 
in terms of Bessel functions with coefficients which are determined by the 
initial conditions. At early times, \kt\ 1, we may neglect the non-standard 
first derivative term and start from the Minkowski vacuum. Using gain the 
variable x = —kt we find in te rms of Hankel functions of the second kind (see 
Abramowitz and Stegunl 1972 ) 



On large scales, x <^ 1, the magnetic field spectrum is given by 

i Z>\ ;:;:i ,+3 " ^ 

If v < 0, i.e. —7/9 < w < — 1/3 we can obtain a red spectrum. Actually, 
we must ask that w < —4/7 in order to avoid an infrared divergence, i.e., to 
obtain n s > —3. 

On the other hand, the scalar spectral index n of C MB fluctuations in 
power law inflation is then given by (see e.g. |Purrerll2008l) 

n _ i= 60+^) i_( n _i)/ 6 ^_ 

l + 3u; l-(n-l)/2 v " v 1 
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With present data iKomatsu et all ( 2011 ). which requires 1 — n < 0.05 we 
cannot reach w > —7/9. However, this argument is not entirely solid since 
w varies slowly during inflation. We know that it has been very close to -1 
when the CMB scales of order several lOOMpc exited the Hubble scale, but it 
may have been larger later, when e.g. the scale of IMpc, relevant for primordial 
magnetic fields exits the horizon. Hence a running spectral index with n ~ 0.96 
at lOOMpc and n ~ 2 at IMpc such that if ~ —0.6 at the time when IMpc 
exits the horizon might be marginally possible. Even though simple running 
is also strongly c onstrained by rin/rflo gfc = —0.022 ± 0.02 at the pivot scale 
of about lOOMpc lKomatsu et all (|201lh . 

However, the CMB results are certainly not compatible with power law 
inflation at constant w and w > —7/9. 

The maximal amplitude is again such that 



PB 
Prad 



H 



inf 



Mf 



(59) 



and m n s ^ —3, a dynamo mechanism after inflation is needed to obtain fields 
with observable amplitude, see figure [T^l 



3.1.3 Helical inflaton coupling 

We can also add a term AL = y r ^gf((f))FF to the Lagrangian, where 

is the Hodge dual of the 2- form F = \F a pdx a A dx^ . In terms of electric and 
magnetic components we have FF = — 4B ■ E. Along the same lines as above 
one can now derive the equation of motion for the gauge potential. In Coulomb 
gauge, writing the Fourier component of the vector potential in the helicity 
basis, A(k) = i (+) (k)e (+) (k) + i ( _)(k)e( _ )(k) we obtain 

2 (±) + [k 2 ± kf]A {±) = . (60) 

It is interesting to compare this equation with (|34[) . The main difference is 
the factor k which replaces here one time derivative. This comes from the fact 
that the mixed term B • E has one time and one spatial derivative. Another 
very important difference is of course the different sign for the two helicities. 
There is always one helicity which will be enhanced and the other which will 
be supressed. 

The new term kfA(±\ is much smaller than the k 2 A(±^ on sub-horizon 

scales and much smaller than the term Ai±\ on super horizon scales. Only at 
horizon crossing in can be relevant. 
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Choosing / such that / = fisr/t with a roughly constant pr e- factor /jy, 
Eq. (|60l) can be solved exactly in terms of Coulomb wave functions (jDurrer et al 



2011 



(x) 



[G (TfN/2,x)+iF (Tf N /2,x)} 



\kt\. 



(61) 



Here Go & n d Fq are the irregular and regular Coulomb wave functions (jAbramowitz and Stegun 
19721 ). and the pre- factors can be obtained by requiring vacuum initial condi- 
tions. From the asymptotics of these functions for small x one finds on super 
Hubble scales 



P B (k) 



sinh(7r/jv) 



N 



P aB {k) = k 



cosh(7r/Ar) - 1 



7T/ 



1 



N 



(62) 

The amplitude of this spectrum can become very large if Jn 3> 1, however, 
then the new interaction Lagrangian dominates over the standard term and 
it is not clear that the perturbative approach adopted here is still valid. The 
spectral energy density grows like k and is dominated by the upper cutoff, 



dp fc 4 sinh(7r/j\r) 



d log k 



4^/ 



N 



Ps(icnd) — H\ 



4 sinh(7r/]v) 



16tt3/, 



N 



Here we use that the Hubble parameter is approximately constant during 
inflation, i? cn d ~ Hini ■ 



PB 
Prad 



^s(^cnd) 



sinh(7r/Ar) / H inf 



58tt3/. 



N 



M P 



~Mp~ 



From this we first conclude that if /n is not too large, back reaction is unim- 
portant since i?i n f <C Mp. However, since this energy density is dominated by 
the contribution at the high-fc end, k ~ H m f , and since the spectrum is again 
blue, n s = 1, we have to draw the same conclusion as in Section 13.1.11 

In this case, however, helicity conservation requires an inverse cascade 
which alleviates the constraints somewhat. An analysis using the ev o lution 
of the spectru m during an i nverse cascade as proposed in ICampanellil (|2007l) 
is presented in lDurrer et all (|2011l) . see also Sections @] and [5] of this work. 

In conclusion we retain: inflation usually leads to a blue spectrum of mag- 
netic fields. It can generate a scale-invariant spectrum only if either the spec- 
trum of scalar inflaton fluctuations is very blue n ~ 1.8 below about IMpc 
(curvature coupling) or if the charge of the electron becomes very large dur- 
ing inflation (inflaton coupling) . This latter conclusion can be evaded if gauge 
invariance is broken during inflation. A helical coupling to the inflation gener- 
ically leads maximally helical fields with spectral index n s = n a = 1. If simply 
scaled to today, such a blue spectrum from the early universe has far too little 
power on Mpc scale to account for the magnetic fields in galaxies, clusters and 
voids. 
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The typical amplitude expected from the quantum generation of magnetic 
fields during inflation is like the one of gravitational waves given by 



P 



(Hi 



inf 



\M P 



VA 



inf 



M P 



(63) 



This fraction can in principle be amplified by dynamo action after inflation, 
e.g. duri ng reheating to n ear equipartition. A possibility which has been pro- 
posed by Sigl et a] (1997) for fields generated at first order phase transitions, 
but which may as well be realized at reheating. 



3.2 Magnetic fields from cosmological phase transitions 



Let us now investigate another possibility, namely that magnetic fields are 
generated during a phase transition. 

Even if the electroweak phase transition is very weak, of second order or 
only a cross-over, magnetic fields with correlation length at the phase transi- 
tion of the order of aA» ~ 1/T ~ 1/mw can form (As be fore, A deno t es co - 
mo ving scales hence the physic al correla tion length is aA*). Vachaspatil (|1991l ) 
and Enqvist and Olesen ( 19931 ) (see also Grasso and Riotto ( 1998f )) have esti- 
mated that these fields have an amplitude of the order of B ~ m^. However, 
the above correlation scale is smaller than the mean free path of particles in 
the plasma and is of the order of the inter-particle distance in the plasma, 
so that one can hardly speak about a persistent magnetic field on time scales 
larger than 1/T in this case. Furthermore, the Ohmic dissipation time on this 
distance scales, rohmic ~ (aA*) 2 er ~ T _1 is many orders of magnitude shorter 
than the Hubble time, H^ 1 ~ Mp/T 2 at the electroweak phase transition 
see Appendix [21 Therefore, even if such fields are generated, they are rapidly 
damped away. 

Magnetic fields with corr elation length significantl y large r than T _1 can be 
generated in the scenario of Joyce and Shaposhnikovl (1997). There the corre- 
lation scale is enhanced by the left-right asymmetry in the leptonic sector. The 
scale of magnetic fields then becomes of the order of aA* ~ ~ (T//i)T _1 , 
w here [i < T is the chemica l poten tial for right-handed leptons. The scenario 
of iJovce and Shaposhnikovl (|1997l ) operates in the temperature range much 
above the electroweak scale, T > 80 TeV. 

At the QCD phase transition, the situation is somewhat different. For 
100 GcV> T > 1 McV, damping by viscosity (see ICaprini et all l2009a and 
Appendix [A| is dominated by the neutrinos with the mean free path aA m f p ~ 
{ZGpT 5 )' 1 - [100 GcV]" 1 (100 GeV/T) 5 , while the magnetic diffusivity is of 
the order a(T)/T, which assures that Ohmic dissipation still damps away the 
fields on small scales aA* ~ 1/T <C aA m f p . Here Gf — l/(292GeV) 2 is the 
Fermi scale and a{T) is the fine structure constant at energy T. Magnetic 
fields on scales shorter than a\ m { p (T) but larger than the Ohmic dissipation 
scale arc frozen in. 
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As long as the electrons are relativistic, both, the correlation length and 
the magnetic diffusivity scale grow like the scale factor. Hence the relation 
fx <C aA„ is maintained. However, as we shall see in Sections [4] and [6l also 
these fields subsequently decay. 

In this section we discuss in some detail the magnetic field spectrum which 
may result from a first order phase transition. A realistic value for the correla- 
tion length, somewhere in-between the extreme values aA* ~ 1 /T and A* ~ £h, 
where £h = l/W denotes the comoving Hubble scale, is still a matter of some 
debate. For a second order phase transition we do expect it to be of the order 
of 1/T. However, if the transition is first order, we expect a correlation scale 
which is of the order of the size of the largest bubbles at coalescence, which 
are of the order of A* ~ O.Olljj. This result has been obtained with numeri cal 
simulations, see iKamionkowski et all (|1994l) ; iHuber and Konstandinl (|2008l ). 

A first order phase transition proceeds via bubble nucleation which is a 
very violent event likely to lead to turbulence in the cosmic plasma. In a 
highly conducting cosmic plasma, turbulence is usually MHD turbulence, and 
a turbulent flow generates both, eddies and magnetic fields in the plasma. A 
de tailed account o f the fascinating field of MHD turbulence can be found e.g. 
in Biskamij ( 2003 ). In this section we shall not enter into any details of MHD 



turbulence but just discuss some generic aspects which will already allow us 
to make very strong statements. 

First, we just note that the two known transitions of the standard model, 
the electroweak transition and the QCD transition are both not first or- 
der. In fact, they a r e not e ve n true phase t ransit ion but just crossovers (see 
e.g. Kaiantie et all (jl996bl laT): Csikor et all (1998) for the electroweak tran- 

sition, if the Higgs mass i s ttih ~ 80GeV and iRoberge and Weissl ( 1986 ): 
de Forcrand and Philipsen (2003). for the QCD transition at vanishing chem- 



ical potential). 



2005 


); 


Huber et all ( 


2007) 



It also has been suggested, that the QCD transition can be first order if the 
neut rinos have a sufficiently l arge, but cosmologically allowed chemical poten- 
tial ( Schwarz and Stukell2009l) . Such a potentia l is eve n required if dark matter 



is to be a sterile neutrino, see iBovarskv et a j (120091) . With this in mind, we 



summarize that taking into account present experimental constraints, it is still 
possible for both, the electroweak and the QCD phase transitions in cosmology 
to be of first order. In this case they lead to the generation of the magnetic 
fields which we now study. 

As has been discussed by IShaposhnikovl (Il987l ) and iTurok and Zadroznv 
( 1990l ). if the electroweak phase transition is first order, it can also explain the 
baryon asymmetry of the Universe. Interestingly, the electromagnetic part of 
the Chcrn-Simons number which determines the net baryon number generated 
at the transition is simply the helicity. This relates the helicity of the mag- 
netic field generat ed at the transition to the baryon number as worked out in 



Vachaspatil (l200ll ). 
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Let us now consider a first order phase transition where correlation lengths 
can diverge. This divergence is of course obtained in a static, thermodynamical 
context where all the modes are in thermal equilibrium. In cosmology the fact 
that the Universe is expanding leads to an effective (comoving) maximal length 
scale A max = t over which correlations can extend. In other words, arbitrary 
correlations £ which are generated in cosmology after inflation satisfy 



£(x, x', f) = if |x-x'|=r>* 



(64) 



Here x and x' arc comoving coordinates, and if the process that generates the 
correlations is statistically homogeneous and isotropic, the correlation function 
£ is a function of r and t only. 

The power spectrum of such causal correlations, which is the Fourier trans- 
form of the correlation function, is therefore analytic for kt ~ 1. For the case 
of magnetic fields this implies that both, (Sij — kikj)PB{k,t) and k m P a B(k,t) 
are analytic at small k. Hence, Pb ~ k ns , P a s ~ k na where n s > 2 is an 
even integer and n a > 1 is an odd integer. The fact that Pb > \P<ib\ even 
requires n a > 3. For more details about this conditions which are simply 
a c onsequence of causa l ity to gether with the fact that B is divergence free, 
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After the phase transition, the magnetic field spectrum has roughly the 
following form: 



Pr, 



for k < fc* 



P aB * f32ir 2 BX 



j^j for k < k„ < k < k d (t) 
for k d (t) < k 

for k < k„ 

for k < k* < k < kd{t) 
for k d (t) < k 



(65) 



(66) 







Here fc* is the correlation scale, k* — 27r/A», with A» < t*, and is the (confor- 
mal) time of the phase transition and ft denotes the helicity fraction. Typically, 
A* ~ t»/100 is of the size of the largest bubbles which form during the phase 
transition b efore coalescence. This is a typ ical number found in numerical sim- 
ulations bv iHuber and Konstand in (200 a), but it depends sensitively on the 
strength of the phase transition (|Caprini et al 20081 : Espinosa et al 2010() . 

In numerical simulations it has been found that h elical magnetic fields be- 
come totally helical soon after the phase transition (see Baneriee and Jedamzikl 
20041: ICampanellill2007l) . so that soon after the phase transition either fj = 1 
or helicity vanishes, /3 = 0. 

The spectral index a which is attained in the so called inertial range is not 
really certain. We shall assume, t hat the turbulence is fully d eveloped and we 
obtain a Kolmogorov spectrum (jLandau and Lifschitd[l990h with a = oi = 
11/3, see Section H] for more details. 
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The energy density in the magnetic field, as given in Eq. (JTUJ), is dominated 
by its value at the correlation scale Xb — 2ir/k*, 

For the density parameter we then obtain with (j45[) 

c PB „ ( 2 \ 1/3 / B, \ 2 

Prad \9es(t*)J U X 10- 6 GaUSsj ' 1 ' 

Here we have taken into account the change in the relativistic number of 
degrees of freedom. We also have assumed that today all neutrinos are massive, 
i.e., m„ > To ~ 2.3x 10 _4 eV for all types of neutrinos so that go = 2. Assuming 
adiabatic expansion one requires a constant entropy, g e gT 3 a 3 = goT§ and 
therefore a 2 p ra( j cx g c gT 4 a 2 behaves like g^ 3 '. 

In the radiation dominated era, the relation between conformal time and 
temperature is given by 



\/3Mp \ 3 x 10 5 sec /100GeV\ 



3 x lCT 3 pc /100GeV\ 



(69) 



Similarly to the case of the inflation-generated magnetic fields, one can 
estimate the strength of the relic fields surviving until the present on large 
scales. If there is no inverse cascade and the magnetic fields evolve passively, 
the field strength at k\ = IMpc" 1 ~ 10 _14 sec _1 is of the order of 



B 



= (fc5 i F s (fc 1 )/(4^ 2 )) 1/2 ~ 3 x 1(T 6 Gauss e^ikx/Kf' 2 

k 1= i Mpc 



Hence passively evolving magnetic fields from the electrowcak phase transition 
can at best amount to about 10 -29 Gauss on Mpc scales while those from the 
QCD phase transition at T* ~ 100 MeV can amount to about 10 -23 Gauss. 

As we shall discuss in Section [SJ magnetic fields generate a spectrum 
of anisotropic stresses which induces a cosmological background of gravi- 
tational waves. The spectrum of this background peaks at frequency = 
fc*/(27r) ~ 100/f*. Interestingly, for the electroweak symmetry breaking scale, 
T* ~ 100 GeV, this corresponds to mili-Hcrtz frequencies which are in the op- 
timal sensitivity range of the pl anned space antenna eLISA (european Laser 



Interferometric Space Antenna) (jBinetruv et alll2012h 



All the above mechanisms are related to some non-equilibrium processes 
(i.e. the phase transitio ns or to the relaxation of the initial conditions). It was 
recently demonstrated (jBovarskv et a 3 l2012bh that already in the Standard 
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Model, long-range magnetic fields can be spontaneously generated as part 
of the equilibrium state in the presence of matter-antimatter asymmetry. At 
finite baryon or lepton at finite density, quantum corrections due to parity- 
violating weak interactions induce a Chern-Simons term in the free energy of 
the electromagnetic field. This result is based on a subtle quantum effect that 
appears in the Standard Model at second order perturbation theory (at two 
loops). This effect can be relevant in both the symmetric phase and the Higgs 
phase. 



3.3 Results from second order perturbation theory 



Within first order cosmological perturbation theory and within the strong 
coupling limit of electrons and protons, no magnetic fields form due to the 
inhomogencities of the matter distribution of the Universe. For this to happen 
we need a current with non-vanishing vorticity. For such a current, J, Ampere's 
law gives 

AB = -— VAJ. (71) 

c 

Within linear cosmological perturbation theory J = e(n p — n e )v, so even if we 
go to 2nd order in the strong coupling limit so that electrons and protons are 
not perfectly coupled and n p ^ n e , since v is a scalar perturbation, hence a 
gradient, we have also to go to second order in the inhomogeneities to obtain 

V A J = eV(n p -n e )Av/0. 

Clearly, such second order perturbations are very small on cosmological 
scales. 

The full system of perturbation equations to second order taking into 
account the imperfect coupling of protons and electrons has been derived 
and st u died numer i cally in several pape r s by Ichiki et all ( 2007 ): Maeda et all 
(|2009h : iFenu et all (|201lh : iMaeda et al (|201lh . Even though the details of 
the results do not quite agree, they all obtain very small magnetic fields, 
B ~ 10 _24 Gauss, on the scales, k ~ 10/i/Mpc where they can calculate the 
field reliably. However, the spectrum is raising towards smaller scales and it is 
not clear whether higher resolution simulations which go up to say l O fe/kpc 
might not give more promising results. (Even though Rcf. Ilchiki et all (|2007j ) 
claim to have a result until k = 10 9 /Mpc, this is just an interpola tion of the 
result found at (1 - 10)/i/Mpc, which has been refuted later by Fenu et al 



(|201lh .1 

Even though these results are not fully under control yet, it seems therefore 
unlikely that straight forward second order perturbations without any initial 
seed fields can lead to sufficient magnetic fields on galactic and inter galactic 
scales. 

This situation can change if first order vector perturbations are present, 
then, the fact that electrons and protons are not perfectly coupled can lead 
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to a different vorticity in each of these fluid at first orde r. The generatio n of a 
magnetic field by this mechanism was first discussed by lHarrison (|l973l) . 

However, in standard inflation such perturbations are not generated, and 
eve n if they are g enerated they decay during the radiation dominated era (see, 
e.g. iDurrerl 120081 ). In order to have genuine, non-decaying vector perturba- 
tions, one either has to source them continually, e.g. with topological defects, 
or one has to modify gravity as, e.g. in the Aether theory of vector-tensor 
gravity. The pert urbative generation of m agnetic fields has been studied for 
both these cases. Hollenstein et all ( 2008 ) have shown that vorticity conser- 
vation prevents the transfer of vorticity by purely gravitational interactions 
which would be needed for the Harrison mechanism to work. Therefore, the 
vector perturbations of topological defects cannot help. Also within the aether 
theo ry, only very sm all magnetic fields of order B ~ 10 _22 G can be generated 



(sec 



Saga et alll2013D . 



4 Cosmological evolution of magnetic fields 



The generation of magnetic fields of strength by a process operating at 
the comoving time i* with comoving correlation scale A* , sets up the initial 
conditions for the subsequent evolution of the coupled magnetic field - pri- 
mordial plasma system from the moment of magnetogenesis up to the end 
of the radiation dominated era and the moment of decoupling / recombina- 
tion. The evolution continues also after recombination in a system where the 
charge density of the plasma is strongly reduced beca use most of the electrons 
and p rotons /nuclei have combined to neutral atoms ( Sethi and Subramanianl 
2005). Plasma effects on the evolution of magnetic fields increase again at the 
latest stages of evolution, when the density of primordial plasma grows again 
after the re-ionization at redshift z ~ 10. 

The qualitative picture of the evolution of magnetic fields and the pri- 
mordial plasma is governed by the MHD equations introduced in Section [2] 
In general, the nonlinear MHD equations are difficult (if not impossible) to 
solve, both analytically and numerically. Some general properties of the so- 
lutions can, nevertheless, be established based on relatively straightforward 
order-of- magnitude estimates of the imp o rtance of the different terms in th e 
MHD equations see Jedamzik et al (1998 ); Subramanian and Barrow dl998a ' 



Jl998at); 

Baneriee and Jedamzikl (l2004l);lBarrow et alM 2007t) : I Jedamzik and Sigll (|201lh ; 
Kahniashvili et aT l 2012 ) ; Saveliev et all ( 2012 ) . In the following subsections we 



summarize these general properties. 



4.1 Initial conditions for the evolution 

The discussion of mechanisms of generation of magnetic fields in Section[3]sug- 
gests that the magnetogenesis results in the production of a turbulent plasma 
and magnetic fields characterized by the power spectra Pg, P a B, PsK, PvK, PaK 
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(see Eqs. (H3J) for definitions), in different intervals of wave numbers. At 
the initial time, the large wavelength (small wavenumber) tail of the power 
spectra of the magnetic field and plasma on scales larger than the character- 
istic scale of magnetogenesis A* are given by 

P B (k,U) = P* B k n ° , P K (k,Q = P, K k nk . (72) 

where Pk scaling applies to P s k and/or P v k- The largest power on the longest 
scales is achieved for the minimal possible values of n s , n^. Formally, for t 
magnetic fields generated at the inflationary epoch, the requirement that the 
magnetic field energy is not infrared divergent restricts n s to n s > —3. How- 
ever, nearly all self-consistent mechanisms of field generation during inflation 
proposed so far satisfy stronger constraint, namely n s ~ 1. As we have argued 
in Section I3~2l see Eq. (|55|) . for causal field generation, e.g. a phase transition, 
we require n s = 2. Different possible spectra of the field at small k are shown 
in Fig. [TJ 

For the velocity spectrum, the situation is more complicated. First of all, 
there is no cosmic plasma present during inflation hence turbulent motions 
develop only after inflation and this in a causal way. This implies that the 
velocity power spectrum defined in Eq. (|13[) is always the Fourier transform of a 
function of compact support and therefore analytic on large scales. The leading 
term in the Taylor expansion of the kinetic power spectrum Pk at small k is 
determined by the nature of fluid motions. As mentioned above, considerations 
of hydrodynamic turbulence often adopt the assumption of incompressibility 
of fluid motions, which imposes a divergence-free velocity field, 

V • v = . (73) 

This condition is identical to the divergence-free condition satisfied by the 
magnetic field, and as mentioned above in this case P s k = P v k = Pk and 
causality requires the same asymptotics for the power spectrum at small k as 
for the magnetic field, 

Pic(k) ~ k 2 , k — > (incompressible fluid/plasma). (74) 

Assuming the validity of condition (fT3")l significantly simplifies the MHD equa- 
tions and facilitates numerical modeling of turbulence. This is why th i s cond i- 
tion is commonly adopted in turbulence modeling (see e.g. Biskamp ( 2003h ). 



In particular, it was adopted for the s t udy o f cos mological magnetic fi elds in 
the papers bv lBaneriee and Jedamzik (2004) and lCaprini et al ( 2009c ). 



However, there is no particular reason why the process of generation of 
magnetic fields in the early Universe would excite only incompressible fluid 
motions. Indeed, it is clear from the system of Eqs. ([27l - 129| that the term 
B A (V AB) in the Euler equation provides a source term for both, compressible 
and incompressible modes. Furthermore, the argument that the Mach number 
M = v/c s = \/3i> <g ; 1 does not suffice, since the additional necessary condition 
see iBiskamplliool . that time derivatives dt are much smaller than the term 
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Fig. 1 Possible spectral energy distributions of cosmological magnetic fields. At small k, the 
spectra of the fields generated at phase transitions in a causal way follow powerlaw with the 
slope n s = 2. Inflationary mechanisms typically result in the slope n s = 1. Inflation could 
in principle generate a scale-invariant spectrum with n s = —3. At large k all the spectra 
follow a universal slope formed by turbulence. 



v • V is not satisfied in our situation. For our relativistic plasma, we expect 
time derivatives which are of the same order as spatiai derivatives. 

Thus, in a generic situation, the power spectrum of the kinetic energy of 
the plasma motions is not restricted to have P s k = P v k- While analyticity 
requires P v k <x k 2 for small k, P s k has no non- analytic pre- factor and is in 
general white noise for small k. Then, the asymptotic of Pk — P s k is white 
noise, 



Px(k) ~ k , k — >• (compressible fluid/plasma). (75) 
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The evol ution of magneti c fields in compressible plasmas is called Burgers tur- 
bulence ( Tsvtovichl 19771) . In cosmological settings i t has been considered by 



Brandenburg et al t 1996 ) ; Jedamzik and Sigl (|201ll ) ; Kahniashvili et al ( 2012 ) . 

On scales k > k* the initial power spectrum of both the magnetic field and 
the kinetic energy of plasma is suppressed. The detailed shape of the initial 
spectra in this regime is usually irrelevant because turbulence establishes a 
"universal" slope of the power spectra on wave numbers k > ks with ks ~ fc* 
which is independent of the initial shape as we discuss in the next subsection. 



4.2 The regime of freely decaying turbulence 

Both the rescaled magnetic and kinetic energy densities, p~b,Pk and the cor- 
relation lengths A b , Aft- defined in Section [2] evolve with time. This evolution 
can in principle be obtained by solving the system of equations (|27H29p . for 
given initial conditions, B(fc,t = t*),v(k, t = t*) and p{k,t = <;*). However, 
the nonlincarity of the evolution equations renders the analytical or numerical 
solution very complicated, for any realistic set of initial conditions e.g. after 
the electroweak or QCD phase transitions. Because of this difficulty, a com- 
mon approach is to derive a qualitative picture of the evolution based on an 
order-of-magnitude analysis of the relative importance of the different terms 
in Eqs. 

The main process determining the evolution of the magnetic fields and the 
plasma at intermediate wavenumbcrs is the establishment of MHD turbulence. 
This process operates in the regime in which the dissipation terms on the 
right hand side of the Euler equation (|2"5|) and on the right hand side of the 
induction equation (j2"9"|) can be neglected. As discussed in Section [5] these are 
scales for which both the magnetic and kinetic Reynolds numbers are large. 
In this situation, the "mode coupling" terms containing spatial derivatives 
in the Euler and induction equations continuously generate larger k modes 
at the expense of the lower k modes and in this way transfer power from 
large to smaller scales. This process is extensively studied in various contexts 
of hydrodynamic and MHD turbulence, both analytically and num erically, see 
e.g. lBiskampI (|2003f ): lTsvtovichl (|l977f ): lLandau and Lifschit3 (|l990h and more. 



The result of these studies is that the power spectrum of the magnetic field 
and kinetic energy evolve to power law spectra given by 

P B oc k- aB , P K oc k~ aK , a, > (76) 

independently of the details of initial conditions. The most commonly known 
result is the Kolmogorov turbulence spectrum with 

a K = -11/3 (77) 

encountered in incompressible hydrodynamic turbulence. The same slope for 
the magnetic field power spectrum, typically as = Q-k is found in numerical 
models of MHD turbulence. The order of magnitude argument that leads to 
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this spec trum is originally due t o | Kolmogorovl (|1941r) . We follow the argumen- 
tation of Landau and Lifschit j ( 1990t) : 



We assume that turbulence is full developed on the scales under consid- 
eration. Energy is transferred from large to smaller scales and dissipated at 
some dissipation scale k^. On scales k < k^ on which turbulence is developed, 
the energy transferred to smaller scales per unit time, let us call it e, can only 
depend on the mean velocity on this scale, v\ = (Pk (fc)fc 3 ) 1 / 2 and on the scale 
itself, A = 2n/k. It must have the units of [e] = [v 2 /t] = [v 3 /X] = [v 3 k]. Setting 
e ~ v\k, this implies 



P„(fc)~£ 2 / 3 fc- n / 3 0cfc- 



with ax = —11/3. 



(78) 



Note that this scaling does not depend on the value of e. The only hypothesis 
used is that e is independent of k which is necessary for the situation to 
be stationary. Even though this is not evident from the above 'derivation', 
the Kolmogorov sp ectrum has also be en observed (numerically) to hold in 
relativistic plasmas ( Mueller et all 2007 ). The part of the magnetic field spectra 
formed by the free turbulence decay is shown in Fig. [1] It is important that 
the large k behaviour of the spectrum processed by the turbulence is largely 
independent of the initial spectrum of the field at the moment of generation. 

If the velocity field and the magnetic field are in equipartition on all scales, 
we expect also as = ax = 11/3 . 



Goldreich and Sridharl (119951 ) use a slightly different ar gument and obtain 
a s pectrum for the m agnetic field with as = 10/3 while llroshnikovl (|1964l ) 
and iKraichnanl (119651 ) propose as = 7/2. The latter values are obtained con- 



sidering the collision of Alfven-waves on a strong background field and are 
probably not relev ant here. Also, simulations are in good agreement with the 
Kolmogorov slope (jMuller and Grappinll2005l ). 

The values cited above are all quite close and the precise values of as , olk 
are, in fact, not important for the general understanding of the process and 
of the time evolution of the power spectra Pb,Pk- It is just important that 
a m < —3 such that the turbulent energy is concentrated on scales around kg 
and not at the damping scale. The evolution of the correlations scale and of 
the energy density can be understood, at least qualitatively, in the following 
way. 

Consider a moment of time t > t*. A generic property of the MHD tur- 
bulence is that it transfers energy from large to small scales. On the scale kd, 
the dissipative terms can no longer be neglected and the turbulent energy is 
lost into heating up the plasma (see below). On very large scales, turbulence 
did not have enough time to fully develop and the initial spectral slope is 
maintained. As time goes on the largest scale on which turbulence is devel- 
oped, the integral scale or correlation scale Xb, Xk grows and correspondingly 
ks = Zk/Xb and kx = 2tt/Xk decrease. Beyond ks the initial spectrum has 
been processes e.g. into a Kolmogorov slope. On k < ks } k < kx the spectrum 
still has its original slope. 
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We expect that the time dependence of these scales follows a power law, 

X B ~ t KB , Ajc ~ t KK , (79) 

with the indices k. > (where • is either B or K) which we now derive. 
The rescaled energy density in magnetic field and in turbulent motions of the 
plasma decreases with time due to dissipation 

pb ~ t~< B , ~ PK ~ r c * , (80) 

with C. > 0. 

Suppose that at a given moment of time the rescaled energy density of 
plasma motions is p~K (t) ■ This energy is associated to a characteristic velocity 
via the relation 

~PK= V \- (81) 

In the same way, a characteristic velocity can be associated to the energy 
density of magnetic field. This is the Alfven velocity given by 

^ = f- (82) 
P 2 

These velocities characterize the speed of the spread of changes in the configu- 
ration of the velocity field and the magnetic field. Typically the size of regions 
over which the fields can change in a coherent manner on a time scale t are 
A, ~ v m t. Unless the magnetic fields are generated in an a-causal way, e.g. 
during Inflation, the integral scale of the magnetic field and of the plasma mo- 
tions at the time t > cannot exceed v m t and, in general, this scale provides a 
reasonable estimate of the integral scales (also called " largest processed eddy" 
scales): 

X B ~ v A t , Xk ~ v K t . (83) 

Coupling between plasma motions and the magnetic field usually estab- 
lishes equipartition between plasma kinetic energy and magnetic field on a 
certain distance scale. In this case the conditions 

p B ~ Pk ~p, v K ~ va ~ v , X K ~ X B ~ A (84) 

are satisfied. We then expect already for reasons of dimensionality that 

A - vt. (85) 

where v is given by v 2 ~ 2p ~ P KlB /X - (A) { where n — n s or n = rik, 
depending on whether matter or magnetic field power spectrum initially dom- 
inates on the scale A. Here n characterizes the slope of the (dominant) unpro- 
cessed part of the power spectrum. Substituting this expression in Eq. (|85D . 
we arrive at the relation 

T - / * 2/5 > when n k = °> n s = 2 => n = 0, . . 

A ~ 1 \ ei\ when n k = 2, n s = 2 => n = 2. (8b) 
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The first case corresponds to compressible turbulence (Burgers turbulence), 
when the power spectrum of plasma motions is Pg ~ k° at small k and the 
power of plasma motions dominates over the magnetic field power Pb ~ k on 
large scales k — > 0. The second case corresponds to incompressible turbulence 
in which Pk ~ Pb ~ k 2 at small scales. 

The energy of plasma motions and of the energy of the magnetic field evolve 
in equipartition according to 



- 2 



Pb ~ Pk 



r 



t 6 / 5 , when n k = 0, n s = 2 . 
t- 10 / 7 , when n k = 2, n s = 2 . 



(87) 



T he characteristic magnetic field strength at t he scale A evolves as (jBaneriee and Jedamzik 



2004 ICampanellill2007t iJedamzik and Siglll201lh 
%b~ 



B 



(3 + ..) 




._. 3/2 

(Aj , when = 0, n s 



2, 



, when nfc = 2, n s = 2 . 



More generally, for an arbitrary n, the dependence of the magnetic field 
strength on the correlation length 



A" 



(89) 



holds with n = min(n s ,nfc). Here we assume that on the scales where tur- 
bulence is developed equipartition between the magnetic field energy and the 
kinetic energy is established. 

In a hypothetical case of the scale-invariant magnetic field generated during 
inflation with spectrum n s ~ —3, the above equation implies that magnetic 
field strength does not change in the course of cosmological evolution, while 
the correlation length grows as Xb oc t. In this case therefore, the correlation 
scale is always the same fracti on of the hor i zon s cale. A case quite similar to 
scaling topol ogical defects, see Durrer et all l)2002h . This scaling has also been 
argued for bv lChristensson and Hindmarshl (|1999l ). without however, realizing 
that it is only valid for scale invariant magnetic field spectra. 

We now consider the case of a maximally heli cal magnetic field, in a cos- 
mic plasma th at respects helicity conservation (jBiskamp and Mulled Il999t 
Biskamplliool . which we expect in the regime where T < m e (see below). 
The helicity density is of the order of H = (A • B) ~ \bPb and its conservation 
implies the relation 

B ~ - A s 1/2 . (90) 

This is equivalent to the substituti on n = —2 in the Eq. The time evolu- 
tion of B, Xb then follows the law ( Biskamp and Miiller 19991) 



A, 



.2/3 



,-1/3 



(91) 



It has been shown recently by iBovarskv et al ( 2012a ). that the chiral 
anomaly of the Standard Model plays an important role for T > m e even 
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though chirality flipping reactions are in thermal equilibrium. Even in the ho- 
mogeneous approximation (neglecting turbulent flows) it was shown that an 
inverse cascade develops solely due to this effect. In the process of this in- 
verse cascade the magnetic helicity is approximately conserved (changes very 
slowly). There f ore it i s an important future project to include the effect of 
Bovarskv et a 3 (l2012ah consistently in the standard MHD analysis. 
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Fig. 2 The evolution of the magnetic field spectrum. Top left: incompressible flow, top 
right: compressible flow, bottom: a fully helical field. The spectra evolve towards smaller k. 
Dashed lines show the kinetic energy spectrum. 



The evolution of causal magnetic field spectra for incompressible fluid mo- 
tions, compressible fluids and for the fully helical case are shown in Fig. [2] 
in the simple case where all components are in equipartition. In the incom- 
pressible case, even though the correlation scale is growing, there is no 'inverse 
cascade' in the sense that no power is transferred from small to larger scales. 
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The growth of As is simply due to the loss of small scale power by dissipation. 
On large scales, A > Xb the power spectrum is not affected. This is sometimes 
called 'passive growth'. This is different for the compressible and the helical 
cases. There, the power spectrum grows on scales A > Xb and decreases on 
scales X < Xb which is what we call an inverse cascade. 

Realistic situations might well be more complicated. For example if the field 
is not maximally helical, or if the scalar velocity mode is much smaller than the 
vector mode, |V • v| <C |V A v|, the field first decays like in the incompressible 
case until it becomes maximally helical, or until the vector amplitude is on the 
level of the scalar amplitude and only then it evolves according to the scheme 
shown in Fig. [5J 

The evolutionary tracks of the correlation length and average strength 
of magnetic field in the B,Xb diagram are shown in Fig. [3] The locus of the 
natural termination points of the tracks is the line corresponding to the largest 
processed eddy size at the end of radiation dominated epoch / recombination 
when the temperature drops to T rec ~ 0.3 eV 

B ~ 10 rrccT ° 3A ~ 1CT 8 ( —^—] G. (92) 
M P \ 1 Mpc J v ' 

The factor ~ 10 accounts for the numerical coefficients obtained when re- 
expressing t and B in Eq. (|83[) through T and Mp , using the Friedman equation 
and the identity 

1(GaUSs)2 = 1.9 x 10- 40 GcV 4 = 3.4 x 10 12 ^ 4 . 
8tt 

This linear relation between the magnetic field amplitude and the correlation 
sc ale, which is simply a consequ ence of Eq. ([85)1 . has first been pointed out 
bv lBaneriee and Jedamzik ( 2004 ). 



4.3 Evolution with strong viscous damping 

In principle, if free turbulent decay proceeds up to the end of recombination 
(when most of the charged plasma disappears and the evolution of the field 
changes), the integral characteristics of the magnetic field move with constant 
velocity along the lines shown in Fig. [3j However, the periods of free turbulent 
decay terminate at least during two regimes where the damping term in Eq. 
([25)1 becomes important at scales up to Xb- 

In the cosmological case, for T < 100 GcV the Prandl number is very large 
so that dissipation of magnetic fields proceeds by kinetic diffusion: the mag- 
netic fields generate velocity fields via the source term in the Eulcr equation 
and the velocity fields are then dissipated via kinetic viscosity. The kinetic 
viscosity is provided by the least coupled particle in the plasma and the vis- 
cosity coefficient in Eq. ([25)1 is of the order of v ~ A m f D /5, where A m f D i s the 



comoving mean free path of the least coupled particle ([Weinberg! 1 197ll ). see 
Appendix [2] for more details. 
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Fig. 3 The evolution of the magnetic field amplitude and integral scale for helical fields 
and for non helical compressible and incompressible flow. Both, the electroweak and QCD 
phase transitions are indicated. The line on which the tracks end is given by the relation 

v A = -^s/trcc, with tree ~ 200MpC. 



Depending on the temperature, the least coupled particles in the plasma 
are neutrinos (for the temperatures above T ~ 1 McV) or photons (for the 
temperatures below MeV down to decoupling at T ~ 0.3 cV). At the time 
of neutrino or photon decoupling, the mean free path of the particles grows 
beyond the horizon scale. The dissipation scale Xd = ^n/kd is the scale at 
which the Reynolds number, Eq. (|30[) . becomes of order unity 

27T u(T) X mip ^p 
k d v K bB 

In the last estimate we have substituted vk ~ va ~ B / -Jp. The growth of A m f p 
at the moments of decoupling of neutrinos and photons leads to the growth 
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of the viscous damping scale, up to values comparable to the integral scale 
Ad ~ Xb — Xk and the fluid enters the so called dissipativc regime. 

The magnetic field and the plasma no longer evolve according to free turbu- 
lent decay, the plasma motions are damped by viscosity already at the integra l 
scale Xk- This regime was studied in detail bv lBaneriee and Jedamzik" ( 2004 ). 



who found that damping at the integral scale suppresses plasma motions and 
removes the coupling of the plasma to the magnetic field (setting v — > in 
the induction equation (|29p ). Once the coupling between the plasma motions 
and magnetic field is removed, the rescaled magnetic field temporarily stops 
evolution, dB/dt ~ 0, so that B and As do not change. To the contrary, the 
kinetic integral scale Xk continues to grow and pk continues to decrease be- 
cause the turbulent kinetic energy of the plasma is efficiently dissipated into 
heating up the plasma. This happens, as long as A^ ~ X m i v /vA ~ Xk 3> A m f p . 
Therefore, this regime is relevant only if va <C 1. 

Once the mean free path of the least coupled particle becomes significantly 
larger than the typical scale of the system (integral scales Xk,Xb), the evo- 
lution changes once more. The least coupled particles are now too weakly 
coupled to the fluid to provide a true viscosity. This situation is called "free 
streaming" . In this regime the viscous damping term i n the Euler equation (1281) 



has t o be replaced by a friction term of the form av (IBaneriee and Jedamzik 



20041) : 

i>V 2 v -> <5v . (94) 
The coefficient a oc A~f p is an ordinary friction or 'drag force' term , anal ogous 



to Silk damping of baryon fluctuations during decoupling (Durrcr 2008). The 
proportionality factor is of order I as long as the electrons are relativistic and 
becomes p-y/ Pb after electron-positron annihilation. 

The growth of A m f p leads to the decrease of a and, as a consequence, 
to a suppression of this damping term in the Euler equation. This provokes 
the end of the dissipative regime and the restoration of turbulence soon after 
the onset of the free-streaming regime. Note, however, that for this to take 
place, it is important that there are stronger interactions (in our situation 
electromagnetic forces) with much shorter mean free path, A m f P 2 <C A m f p so 
that the fluid picture still applies, otherwise the Euler and continuity equations 
have to be replaced by a Boltzmann equation. When the next largest A m f P 2 
grows to the value A m f P 2/Ag ~ a, the original Navier-Stokes form of the 
damping term term is re-installed with v ~ A m f P 2 and the coupling to the 
weakest coupling particle species can be neglected. 

The kinetic energy of plasma motions is dissipated into heat in the dissi- 
pative regime. The energy contained in the magnetic field is constant, up to 
the trivial dilution due to the overall expansion of the Universe on scales much 
larger than the Ohmic dissipation scale given by 

M^ 1/2 



Aohmic ~ (oU)- 112 ~ V"l°g(« _1 ) ( if ) T o _1 ; (95) 
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the rescaled magnetic field is not evolving when the velocity field vanishes: 
dB/dt ~ 0. Once turbulence is restored after decoupling of the most weakly 
coupled particle , the Lorentz force term in the Euler equation ((28|) 



B A (V A B) (96) 

serves as a source term for plasma motions. The coupling between plasma and 
magnetic field restores also equipartition between the magnetic and the kinetic 
energy, so that a new cycle of free turbulent decay starts. 

The equations determining the correlation length A ~ Xb ~ A^ and the 
energy density / velocity scale v ~ vk ~ va in the free turbulence decay 
regime are largely insensitive to the details of the shapes of the kinetic and 
magnetic power spectra. The only property of the spectrum which matters is 
the total energy density, which determines the average velocity scale v and, 
as a consequence, the "eddy processing" time scale t ~ X/v. Since, the en- 
ergy density of the magnetic field right after the restoration of the turbulent 
regime is roughly the same as just before the end of the turbulent regime 
at transition to the viscous regime, further evolution of the system in the 
restored t urbulence regime proceeds " as if there was no episode of viscous 



damping (jBaneriee and Jedamzikll2004l) . In terms of the evolutionary diagram 
for B,Xb, shown in Fig. [3J the system always remains on the same track 
(shown by the lines with arrows). During free turbulent decay, the system 
moves along the track. During the viscous damping and free streaming (the 
dissipative regime) the system halts and remains at the same point of the track 
B = const, Xb = const until free turbulent decay is restored. 

This type of evolution is also observed in numerical simulation s both for in- 

comp ressible ( Baneriee and Jedamzikl2004f) and for compressible (iKahniashvili et all 



20121 ) 



2012) MHD and also in the helical case (jBaneriee and Jedamzikl 2004: Ka hniashvili et all 



4.4 Example: evolution of the field generated at electroweak phase transition 

As an example, let us consider magnetic field with initial comoving correlation 
length A* ~ 0.01^ (here £h = l/'H* = t* is the comoving Hubble scale at 
time i*) and magnetic energy density ps ~ O.lp produced at the electroweak 
phase transition at temperature T* ~ 100 GeV (see Section 1372]) . Wc assume 
that the evolution of the magnetic field is governed by the compressible MHD, 
so that the correlation length evolves with comoving time aa^ 

Xb=X.J -) ~ 0.01^- ( -) ~ 2/5 ~ 10 14 ( -^—) ~ 2/5 cm . (97) 

\tj t*t \tJ Vioocevy v ' 



2 Even if we have initially As = Xi((t — ti)/r) 2 / 5 after a few Hubble times this very 
turns into As = A* (t/t» ) — 2// ' 5 , where A, denotes the correlation scale at t». 
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In the temperature range 100GeV> T > lMcV the dominant contribution to 
the viscosity is provided by neutrinos with the comoving mean free path 



1 / T 



W = Gp5jr*HiooGSvJ (98) 



The velocity v ~ y pb/p evolves as 



, v -3/5 / j, \ 3/5 



^^0.3^-j -0.3^^) 

so that the Reynolds number at the scale Xb evolves as 



v aX b . 1n i2 ( T 



This means that the free turbulence decay terminates (i.e. Rk ~ 1) when the 
temperature reaches T ~ 0.1 GcV. Starting from this moment the rescaled 
magnetic field and Xb stop evolving. 

Once the mean free path of neutrinos becomes significantly larger than 
Xb, the damping term in the Euler equation changes to i>V 2 v — > av with a ~ 
A~f pi/ . This means that the Reynolds number (which expresses the relative 
importance of the damping term compared to the (v • V)v and/or B A (V A B 
terms) becomes 

Rf ., = ^ = 10 -.y_)-\ (101) 

Thus, at T ~ 8McV, the system returns into to free turbulent decay, corre- 
sponding this time to large R^ 00 -*. The fact that in this case the dissipative 
regime is relatively short comes from the fact that va is rather large and 
j^(free) ^p^-i becomes larger than unity after a damping regime which 
lasts only about a decade in temperature (and conformal time). 

At T ~ lMeV, the neutrinos finally decouple from the cosmic plasma After 
that time, the main contribution to the viscosity is provided by the photons. 
Assuming that the electrons are non relativ istic, their ab undance is suppressed 
by a factor rj b = n b /n 7 ~ 2.7 x W~ 8 Q b h 2 (jPurrerlbOoj) . 

The mean free path of photons then is 

Amf P , 7 = -— ^ ^ - 10 11 (t^tt;) cm , (102) 

a(T)aTn e rjbcrTT^To \lMeVy 

and the Reynolds number corresponding to the photon viscosity becomes 

v aX b ,„/ T 



Rk = A^^»(TMevJ • < 103 > 

so that the free turbulent decay terminates again at temperatures around 
T ~ 0.3 MeV because the mean free path of photons becomes comparable to 
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vaXb- Then the magnetic field correlation length stops evolving for a while, 
until the photon mean free path does become significantly larger than A^. 
Once this is the case, the damping term in the Euler equation takes its free- 
streaming form i>V 2 v — > av, but now with with 



{Pb+e± + Pb+e±)X m fp 3pb 

where pt>+ e ± is the density of the baryon + electron/positron fluid and p~, is the 



photon density. The last equality is valid when T < 0.5 MeV (jSubramanian and Barrow 



1998a) 



Similarly to the neutrino decoupling regime, the importance of the heat 
dissipation via photons now decreases and free turbulent decay starts again. 
However, this restart of turbulence is somewhat delayed by the large factor 
Pi I Pb ~ 10 6 (T/MeV). This time the viscosity is provided by the electrons, for 
which the growing Coulomb collisions cross-section 



(^)V (105) 



2 

/ ••"e 



assures that their comoving mean free path is shorter than that of photons. 
We obtain 

Amfp.c = {cr c n e a)~ l = — — ~10 n cm, (106) 

rj b a c T z T 

which is independent of the temperature. 

The natural termination point of this new free turbulence decay period is 
the moment of recombination, t rcc when most of the charged plasma disap- 
pears. This is the termination time adopted in Fig. [3J Xb/va = tree — 200Mpc. 

A similar discussion of the ma gnetic field evolution has first been presented 
bv lBaneriee and Jedamzik (|2004h . 



4.5 Late evolution in the matter dominated Universe 

Evolution of the magnetic field modes via MHD interactions with the charged 
plasma is, in principle, possible also at later stages of evolution of the Universe, 
i.e. in the matter dominated era. Indeed, some residual charged plasma is still 
present after recombination. Furthermore, the intergalactic medium becomes 
completely re-ionized again after reionization of the Universe at the redshift 
z r i ~ 10. This means that MHD processes, including turbulent and damped 
decay can operate also at the late stages of the evolution of the Universe. 

Processing of modes via MHD turbulence during the late time can be 
qualitatively understood in the same terms as during the radiation dominated 
era. The equation 

X B ~ v A t (107) 

for the size of the largest processed eddies does not depend on details of the 
evolution of the Universe and is applicable also during the matter and cos- 
mological constant dominated eras. The main qualitative difference with the 
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radiation dominated era comes from the time evolution of the velocity scale 
v A = v / 2p B /p- 

To understand this difference, it is convenient to consider first the most 
optimistic case when magnetic field energy density decreases only due to the 
expansion of the Universe (i.e. no processing via MHD turbulence occurs). 
During the radiation dominated era, ps and p = /9 ra( j evolve in the same 
way so that va — const. However, during the matter and /or cosmological 
constant dominated era, ps evolves like the radiation energy density, while 
p = pb oc a -3 . Here p^ is the baryon density. Dark matter is decoupled from 
the charged plasma and does not participate in MHD turbulence. This means 
that Pb/ P decreases with time and, as a consequence, 



Thus, no further growth of the comoving magnetic correlation scale occurs 
during the matter dominated era. If we still assume that the magnetic field 
strength decreases as a result of dissipative processes, va decreases even faster 
with comoving time and the size of the eddies which can be processed via 
MHD even decreases. 

It is interesting to note that Eq. (|107|) also applies to magnetic fields which 
might be generated in the late Universe (e.g. via Galactic winds, see Section[S] 
If relatively strong magnetic fields generated at rather short distance scales are 
ejected into the intcrgalactic medium, they drive turbulence in the intergalactic 
medium, via the Lorentz force term B A (V A B) in the Euler equation (|28[) . 
This leads to an processing of the magnetic eddies by MHD turbulence up to 
the scale given by Eq. (|107|) . The process of free turbulent decay then moves 
the magnetic field power to shorter scales where it is dissipated into heat. 
Thus, the strength of the magnetic field is reduced and its correlation length 
is increased until the condition 



is satisfied, where to is the present age of the Universe. Actually, the time in 
the above equation should be to — U where ti is the injection time, but we 
assume U -C to, hence injection at redshift z, > 1 so that we may neglect this 
correction in our order of magnitude estimate. This provides an upper bound 
on the fields injected at scale A into the intergalactic medium by any process 
(including relic fields produced in the early universe) 




Xb ~ constant. 



(108) 




(109) 



B < H^qI /2 M p X b ~ 10" 



8 



Xb 



G. 



(110) 



1 Mpc 



Occasionally, this limit is close to the limit (|9"!?|) on the comoving field strength 
at the moment of recombination. 
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The viscosity of non-rclativistic plasma (jLifschitz and Pitaievskilll983f ) is 
v ~ vt A m fp, where vt — (T/rrie) 1 / 2 is the thermal velocity of the plasma elec- 
trons. The Reynolds number is, therefore, Ra' = 5vA^B/{vT^mfp) ~ v\to/(VT^n 
In order for these magnetic fields to be in the turbulent regime we therefore 
have to require also 

V T Xmip < v\t . (Ill) 
The electron mean free path is dominated by the Coulomb collisions such that 

\ ~, (- „ \-l ~ l0T IGM m P „, 5 P C ( Tigm \ 2 
A mfp _ (acn.j _ aTm 2 nbH 2 M 2 - (1 + z)2 ^ 1Q 4 K; ) ■ 

where Tjgm ~ 10 4 K is the present-day temperature of the IGM, which got 
re-heated by the reionization at the redshift z r j ~ 10. Inserting this in the 
expression for the Reynolds number we find that the condition > 1 requires 

B>10- 12 G(^)~ 10/6 . (112) 

Thus, fields with correlation length larger than Xb ~ 0.1 kpc (see Eq. (JTTDJ) ) 
could have been processed by the turbulence in the IGM. 

Weaker fields with correlation length in the range 5 pc< Xb < 100 pc 
might avoid damping via turbulence at the late stages of evolution, because in 
this distance range IGM plasma velocities are damped by the viscosity due to 
Coulomb collisions. However, at the onset of reionization, when the tempera- 
ture of the IGM was much lower than ~ 10 4 K, the dampling distance scale 
might have been shorter, so that fields with shorter correlation length possibly 
excited turbulence. Detailed understanding of the IGM turbulence excited by 
the magnetic fields would require modeling of the reionization dynamics. 

At still shorter distance scales, Xb < A m f p , the IGM plasma is collisionless, 
so that its not appropriately described by the MHD equations. Instead , one has 



to res olve the Botzmann /Vlasov equations for particle distributions (jKulsrud 



19831 ). Turbulence could also develop in the collisionless plasma. The generic 
nature of the relation (|110[) suggests that it might be also applicable to the 
collisionless case, although a detailed investigation of the behaviour of the 
IGM in the presence of relatively strong short scale magnetic fields is needed 
to verify this. 

In general, magnetic fields which are initially stronger than the limit (|110[) 
could evolve along the evolutionary tracks outlined above toward this limit. 
Therefore, the line in the (B,Xb) plane, given by Eq. (|109|) is the locus of 
the "termination points" of the evolution tracks for all cosmological fields 
generated at short distance scales. 

Further pro cessing of the fie l ds by the MHD processes occurs duri ng struc- 



ture formation (|Rvu et alll200l ISchleicher et allBoiol ISur et alll2010h . In gen 



eral, gravitational collapse leading to formation of galaxies and galaxy clusters 
amplifies any pre-existing fields via straightforward magnetic flux conserving 
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compression and/or via action of various types of dynamos. In galaxies, flux 
conservation yields an amplification of 



where p mi pfi n are the initial and final average matter densities before and 
after the gravitational collapse. 

The amplification by dynamo action is much more uncertain, but it may 
be many orders of magnitude larger. In this case, the final characteristics of 
the fields in the gravitationally collapsed structures is largely independent of 
the initial conditions at the onset of the structure formation, but is given by 
some dynamo saturation amplitude. The fact that magnetic fields in galaxies 
are all roughly of the same amplitude hints that this may well be the case. 

The only place where the "relic" initial magnetic fields are preserved is then 
the intergalactic medium in the voids of the large scale structure. The pre- 
existing field in the voids is not processed by the MHD effects accompanying 
structure formation and, therefore, the field in the voids must still satisfy the 
relation (|110|) . Detecting magnetic fields in voids with correlations scale given 
by this relation would be a strong indication of their primordial nature. 

5 Observational constraints 

The range of the field strengths B\ and correlation lengths Xb, implied by 
the relation (|110[) is within the reach of available observational tools based on 
the methods of radio and 7-ray astronomy. Therefore, it appears reasonable 
to explore the possibility of observational detection of the relic fields which 
might be present in the voids of the LSS. In this Section we summarize the 
status of the searches of Intergalactic Magnetic Fields (IGMF), including the 
fields in the voids of the LSS. 

In the absence of positive detections, the discussion of this section is limited 
to the summary of observational constraints on the strength and correlation 
length B\,\b of IGMF. It is important to note that relic magnetic fields from 
the Early Universe are not the only magnetic fields which might populate 
the voids. Therefore, even a real measurement of of IGMF, which should be 
possible with future observational facilities, does not necessarily imply the 
measurement of the relic fields from the Early Universe. We discuss the possi- 
bility to distinguish the relic magnetic fields from the fields of different origin 
at the end of the Section. 

Similarly to the previous sections, we present the observational constraints 
on IGMF in the (B\,Xb) parameter space, see Fig. [3j This figure shows the 
evolution of the field strength and correlations scale throughout the history of 
the Universe. On the other hand, it can also be used to show constraints on 
present day fields, i.e., limits on the allows ranges for the parameters (B\, Xg) 
at z = 0. The white, unshaded area in Figs. [3] to [5] and similar figures shows 




(113) 
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Fig. 4 Theoretical constraints on the IGMF parameters in the present day Universe. 

the allowed range of parameters of IGMF in the present Universe. Each of the 
constraints (boundaries of the unshaded region) is explained in detail in this 
section. 



5.1 General theoretical bounds on B\, Xb 

The straightforward theoretical constraint on the present day strength and 
correlation length of cosmologically produced IGMF is given by Eq. (|110[) . 
Strong magnetic field injected at small distance scales would drive turbulence 
in the primordial plasma and later in the IGM. Eq. ()109[) provides an estimate 
of the size of the largest eddies which can be processed by the turbulence on a 
time scale comparable to the age of the Universe. Turbulence removes power 
from the short-scale modes of magnetic field. This leads to the increase of 
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Fig. 5 Constraints on the IGMF from Faraday rotation measurements. 



the field correlation length until the relation (jllOp is satisfied. The constraint 
(fTT0)l is shown with the label "MHD turbulent decay" in Fig.|4] 



There is no formal upper limit on the possible correlation length of IGMF. 
If generated during inflation, it might be even coherent on the scales larger 
than size of the visible part of the Universe. However we can never observe 
correlations on scales larger than the present Hubble scale, which is indicated 
by the line marked " Hubble radius" in Fig. [U A fields with correlations scale 
larger than the present Hubble scales would be perceived as a constant field 
throughout the Universe. Observable limits on such fields are not so strong, 
they are of the order of 10 _9 Gauss and mainly come from Faraday rotation 
in the CMB polarization, see Section 15751 A recent revised discussion on CMB 
limits for constant magnetic fields, taking into account free streaming neutri- 
nos, can be found in lAdamek et a 1 <l201lh . 
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5.2 Faraday rotation measurements 



Upper bounds on the strength of IGMF are imposed by the non-observation of 
Faraday rotation of the polarization plane of linearly polarized radio emission 
from distant quasars. 

Propagation of a linearly polarized electromagnetic wave of wavelength 
A through plasm a causes a rotation of the polarization vector by an angle 
(|Kronberdll994l ) 

& = RM A 2 (114) 

where the rotation measure (RM) is determined by the distance to the source 
d(z), the strength of the magnetic field component parallel to the line of sight, 
Bii and by the free electron density n e in the region through which the wave 



is propagating 



RM 



d ^ n e (z)B«(z 



-dx(z) . 



(115) 



In standard /1CDM cosmology, the distance element dx(z) is related to the 
redshift by 

dz 

dx{z) = = . (116) 

H (l + z)^/(2 m (l + z) 3 + n A 

Q m , f^A are the present matter and dark energy fractions. 

Given the distribution of free electrons in the IGM one can derive a mea- 
surement (or an upper limit) of £>| from the measurements of (or upper limits 
on) the Faraday rotation of polarized radio emission from distant cxtragalactic 
sources of (most of them are distant quasars). 

There are several challenges to the measurement of IGMF from Faraday 
rotation. First, the distribution of free electrons in the IGM along the lines 
of sight toward different quasars is uncertain. Depending on the assumptions 
about thi s distribution, different bounds on the IGMF have been rep orted in 
the past (|Kronberelll994llKronberg and Perrvlll982l: iBlasi et all["l999h . 

Furthermore, the effect of Faraday rotation due to the IGMF is small com- 
pared to that produced by the magnetic field of the Milky Way. Sensitive con- 
straints on the IGMF can be derived only after a proper characterization and 
subtraction of the effect of Galactic magnetic field. However, our knowledge of 
the Galactic magnetic field is rather limited (Han and Qia olll994t iBrown et all 



2007 : Pshirkov et all2011 : Jansson and Farrarl2012af b: Qppermann et a lll2012h 
and there still are large uncertainties in Galactic field models. This introduces 
uncertainties in the constraints on the IGMF derived from the Faraday rota- 
tion measurements. 

In the simplest approximation, an estimate of the free electron density in 
the IGM can be obtai ned from the known mean baryon density p b = ftbP, with 
n h ~ 0.02//i 2 ~ 0.04 dKomatsu et afeOllh . Here and in what follows we adopt 
the value Hq = /il00km/s/Mpc = 70km/s/Mpc for the Hubble constant. The 
IGM is almost completely ionized today, so that 

Q b p b 



(117) 
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is a good estimate of the average free electron density in the Universe (and in 
the voids of LSS) today, implied by the electric neutrality of the Universe. This 
estimate has been ad opted in the early studies of the Faraday rotation con- 
straints on the IGMF ( Rees and Reinhardtlll972 ; Kronberg and Simard-Normandinl 
19761: lKronberelll994h . 



The rcdshift dependence of n e depends on the ionization history of the 
interstellar medium. After z ~ 10 the IGM is mostly ionized so that n e is about 
the average baryon density n e (z) ~ n e o(l + z) 3 . Assuming that the magnetic 
field strength at the integral scale decreases only due to the expansion of the 
Universe, B\\(z) ~ B (l + z) 2 and substituting this redshift dependence of -By 
and n e into Eq. (|115p gives numerically 



-R^igmf — 10 



10" 8 G 



n e o 



10" 



y/f2 m (i + z)^ + n A -i 

firm, 



rad 
m 2 
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Hn 1 . For fields 



for magnetic fields coherent on the Hubble scale Xb ~ (-H — "o 
with smaller coherence scale, a suppression factor of {Xb / (-h) 1 ^ 2 is introduced 
to account for the randomness of the magnetic field direction. Constraints on 
the IGMF at the level of 



( -7— 

"h 



-1/2 



G 



(119) 



have been derived in this way from the Faraday rotation data by flRccs and Reinhardt 



1972UKronberg and Simard-Normandinlll97alKronberdl9"9llBlasi et aJ1 999'). 



(Note t hat stronger constraints are quoted in the o riginal literature. iRees and Reinhardt 
(<1972h : 1 ronberg and Simard-Normandinl (|1976l ). since there Qj, ~ 1 is as- 
sumed.) This constraint, properly rescaled to fib — 0.04, is shown by the dark 
blue shading in Fig. [3] 

A more elaborate analysis leading to a somewhat more accurate estimate 
of n e along the line of sight to distant quasars i s based on the account of 
information obtained from the Lyq f orest data ([Kronberg and Perrvl 11982 ; 
Blasi et a ll999h . lKronberg and Perrvl (|l982l) were first to notice that an excess 
of Faraday rotation in the signal of distant quasars is observed when the line 
of sight toward a quasar passes through absorption line systems. They have 
interpreted this as being due to the presence of magnetic fields in the " clouds" 
responsible for the absorption lines (these may be either intervening galaxies, 
or lower density Lya clouds which hav e not yet collaps e d to f orm galaxies). 

This idea was further developed bv lOren and Wolfe (1995) who also made 
an attempt to subtract the Galactic RM from the RM measurements of distant 
quasars, to search for the residual RM due to the IGMF and/or magnetic fields 
in the intervening clouds / galaxy systems. Their claim is that the contribution 
of the Galaxy to the RM of distant quasars can be determi ned with p r ecisio n 
better than ~ 15 — 20 rad/m 2 . This result has been used bv lBlasi et al (1999). 
who has introduced a detailed model of electron density distribution along the 
line of sight toward quasars, based on the statistics of the density distribution 
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in Lya clouds. IBlasi et a 1 (1999) also assumed that in higher density clouds 
the magnetic field is amplified by compression, to that the effect of the density 
increase in the clouds contributes twice to the Faraday rotation signal (|115[) : 
once directly through the increased n e and second time through the increase of 
-B|| . Stronger Faraday rotation signals from the clouds results in stronger con- 
str aints on th e unam plified IGMF outside the clouds. The constraint derived 
bv lBlasi et a 1 (1999) is shown as the light-blue shaded region in Fig. [S] 

Another application of the correlation of excess rotation measures with 
i nterv e ning structures alon g the line of sight is developed in iKronberg et al 
(|2008h : lBernet et all (|2008tl . They find strong correlation of the increased RM 
with Mgll absorbing systems with the equivalent width (EW) of the Mgll line 
of EW > 0.3 A. The excess RM introduced by the intervening systems is only 
RM ~ 140 rad/m 2 . Measurements of so small additional RM are possible be- 
cause of the redshift dependence of the additional RM due to the intervening 
Mgll absorption systems. The typical size of the Mgll absorption systems is 
~ 100 kpc, which indicates that these systems are, most pro bably, bubbles 
aroun d the star- forming galaxies produced by galactic winds ( Bordoloi et all 
201 lh . Based on the measurement of the RM and on t he estimates of the hy - 
drogen column den sities of the Mgll absorption halos, Kronberg et all ( 2008 ): 
Bernet et all ( 2008 ) derive an amplitude of the magnetic field in these 100 kpc 
scale galactic halos of B ~ 10 fiG, under the assumption that the coherence 
length of the magnetic field is comparable to the halo size. The field correlation 
length in the halos is most likely significantly shorter than the halo size. Pres- 
ence of the field reversals would boost the the estimate of the field strength 
by the square root of the rati o of the halo size to the correlation length (see 
Bhat and Subramanianl ( 2013 ) for further discussion). 

The main uncertainty in the measurements of magnetic fields in different 
components of the LSS, based on the Faraday Rotation technique comes from 
the uncertainty of the Galactic contribution to the RM signal. Substituting 
typical scale height of the Galactic disk -ffcai for d(z) and the free electron 
density in the interstellar medium n e ,Gai for n e (z) in Eq. (| 1 1 5[) one finds a 
Rotation Measure due to the Galactic magnetic field of the order of 



RM Ga 



10 z 



D. 



Gal 



10" 6 G 



0.1 cirr 



H, 



Gal 



1 kpc 



rad 
m 2 



(120) 



where £?Gai is the magnetic field strength in the interstellar medium. Although 
"typical" values of all the three quantities, -ffcai, £?Gai, "eGai are difficult to 
define (e.g. they arc widely different for the Galactic disk and halo), the order of 
magnitude estimate in Eq. (|120p shows that the Faraday rotation accumulated 
during the propagation of the radio beam through the Galaxy is much larger 
than the that accumulated during the propagation through the IGM. 



5.2.1 Prospects for IGMF measurement with next- generation radio telescopes. 



Recent accumula tion of a large da t a base of Rota tion Measures from extra- 
galactic sources (jTavlor et a 1 2009; IStil et all 120111) has enabled a significant 
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improveme nt of the knowledge of the global structure of the Galactic mag- 
netic field ( Jansson and Farrar 2012al lb; Qppermann et al 2012 ). This will, in 
principle, allow a better control of the Galactic RM and, as a consequence, 
lead to better constraints on the IGMF contribution to the RM. However, an 
order-of-magnitude improvement in the sensitivity of the Faraday Rotation 
measurements of the IGMF requires to shrink the error bars of RMq^i by 
a factor of 100, which in term requires the knowledge of the Galactic mag- 
netic field and free electron density with sub-percent precision. Taking into 
account the remaining large un c ertainties in modelin g the free electron distri- 
bution (jCordes and Lazia2002 ; Gaensler et all 20081) . as well as in the degen- 
eracy of model parameters o f the Galactic magnetic field (jJansson and Farrar 
2012aH bl; |Psnirkov et al|[201lh . it is not clear whether this precision can indeed 
be reached. 

The next qualitative improvement of our knowledge of the three-dimensional 
structure of the Galactic magnetic field and of the free electron distribution in 
the interstellar medium is expected from the next generation radio telescopes 
LOFAR (ht tp : //www. lof ar . org/) and SKA (http : //www . skatelescope . org/), 
(|Becldl201lh .' Qualitatively new survey capabilities of these facilities will fur- 
ther increase the sample of RM me asurements for extragalactic sources from 
~ 1 source/deg 2 ( Tavlor et all 20091 ) up to ~ 10 3 sources /de g 2 in the case o f 
an SKA all-sky survey with 1 hr exposure per field-of-view dGaensler 20061) . 
This will allow a much more detailed modeling of the Galactic magnetic field, 
thereby improving the sensitivity for the search of weak extragalactic contri- 
butions to the RM. 

An improved measurement of the three-dimensional structure of the Galac- 
tic magnetic field can be obtained by taking into account not only extragalactic 
sources, but also sources of linearly polarized radio emission inside the Milky 
Way. The most important class of polarized Galactic sources are pulsars. Mea- 
surements of dispersion and rotation measures of the pulsar emission provide 
constraints on both the magnetic field and free electron density in the inter- 
stellar medium. Pulsars can be found at different locations inside the Galaxy, 
so that they in principle allow for a three-dimensional "tomography" of the 
Galactic magnetic field (provided that sufficiently large num ber of pulsars ca n 
be found in thin distance slices and in different directions ( Han et al 20061) . 
Up to now some ~ 2 x 10 3 pulsars are known. 554 of them have been used by 
Han et al ( 2006 ) to study the structure of the Galactic magnetic field. SKA 
will provide a qualitative i mprovement due to a 10 times larger (~ 2 x 10 ) 
pulsar detection statistics (jSmits et a 3 l2009h . 

A qualitatively new possibility to distinguish the IGMF contribution to 
the RM from the Galactic contribution will also arise with some ~ 10 s extra- 
galactic sources in the SKA sky survey. This will be a possibility of a three- 
dimensional "RM tomography" of the Universe, i.e. study of the gradual ac- 
cumulation of the RM signal in thin redshift slices. The most straightforward 
effect expected from the IGMF contribution to the RM is the characteristic 
dependence of the signal on (1+z), sec Eq. (I118p . Detection of such dependence 
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in the RM(z) signal may provide a possibility to "bypass" the uncertainty of 
the Galactic contribution to the RM . 



5.3 Limits from CMB observations 



Magnetic fields interact with the primordial plasma in the early universe. The 
presence of sufficiently strong magnetic field therefore, affects the evolution of 
the plasma. The imprint of magnetic fields on the state of the plasma can po- 
tentially be revealed in the properties of cosmic microwave background (CMB), 
which encodes information on the state of primordial plasma at the epoch of 
photon decoupling and recombination, z rcc ~ 1100. The observed temperature 
fluctuations and polarization provide the most precise cosmological data set 
and it is therefore most interesting to study the effect of a primordial magnetic 
field on these data. 

A m agnetic field affects t he CMB anis o tropies and polarizatio n in many 



ways, see 



Barrow et al dl997h:lDurrerl (l2007MShaw and Lewis! (|2010h : |Paoletti and Finellil 



(|201l[ ); IShaw and Lewis! (j2012h : IPaoletti and Finellil (l2012t). The result s for 



constant magnetic field derived in the pioneering paper bv lBarrow et a 1 (1997) 
are actually invalid, since there the compensation by neutrino anisotropic 
stresses (see lAdamek et alll2011l) which isotropize the Universe are not taken 
into account. 

The energy momentum tensor of the magnetic field perturbs the geom- 
etry of the Universe which governs the geodesic motion of CMB pho- 
tons. This introdu ces scalar, vector and tensor perturbations in the CMB 
( Durrer et~all2000h . 

The evolution of the cosmic plasma is affected by the presen ce of a mag 



1. 



netic field which leads to fa st and slow magnetosonic waves ( Adams et al 



1996) and to Alfven waves (|Subramanian and Barrowlll998b 



Durrer et al 



1998) . The former introduce sl ight shifts in the acoustic peaks of the CMB 
(jKahniashvili and Ratra 2007 ). while the latter mainly lead to vector per- 



turbations ( LewisH2004ir 
Faraday r otation turns -E-polarization of the CMB partially into B-po- 
larization ( Seshadri and Subramanianl2001 ; Kahniashvili et a ll2009h . bmce 
Faraday rotation is frequency dependent, see Section I5T21 this can be sep- 
arated from the the effects under points [TJ and [2] which are 'achromatic'. 
Actually magnetic fields generate large vector modes which generate domi- 



nantly B-polariza tion, but with the usual thermal CMB spectrum fsee lSeshadri and Subramanian 



200 It [Lewis 2004) 



If the magnetic field is helical, its parity violation leads to correlations of 
the temperature anisotropy and of E-polarization with B-pol arization; cor- 
relati ons which are forbidden in a parity invariant universe ( Caprini et all 
20041) . 



5. The presence of a magnetic field affects recombination an d Silk damping. It 



there fore alters the damping tail of the CMB anisotropies (jJedamzik and Abel 
20111) . 
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Non-thermal dissipation of magnetic field energy into the energy of elec- 
trons/positrons during the recombination epoch can lead to distortion of 
blackbody CMB spe ctrum, mainly by introducing a chemical potential 
(jjedamzik et allll998h . 



Magnetic fields affect the formation of clusters. Their abundance is well 
deter mined by Sunyaey-Zerd ovich (SZ) decrement measurements in the 
CMB (|Shaw and Lewisll2012| ). 



8. Since the energy momentum tensor of the magnetic field and the Lorentz 
force in the MHD limit are quadratic in the field strength, magnetic fields, 
even if they are Gaussian, will introduce non-Gaussian CMB anisotropics 
and polarization. 

Interestingly, all these effects yield limits on magnetic fields of the order of 
nG. This is not so surprising as the fluctuations in the CMB arc of the order 
of 10~ 5 and the energy density in a cosmic magnetic field is 

^ 10 ~ 5 (iF^g)~^- (121) 

We therefore expect that magnetic fields of 10 -9 G leave an imprint of about 
1% on the CMB anisotropics and polarization, which is marginally detectable. 
The only effects on the CMB for which this argument is not valid are the 
points [5] and |B] of the above list. 

In the following we explain the above points in more detail. 



5.3.1 Limits from CMB angular power spectrum 

The energy momentum tensor of a stochastic magnetic field is fluctuating 
from point to point. Its power spectrum is given by the 4-point function of the 
magnetic field, e.g. 



« 4 



(T}f\k)T< B \V)) = (2^) 3 <5(k-k') J d 3 9 ( J B ! (k-q) J B,(q)i?,(q)i? m (k-q)) + . 

(122) 

In particular, even if the magnetic field distribution is Gaussian, its energy 
momentum tensor is not. But in this case the power spectrum of the energy 
momentum tensor can be expressed in terms of the m angetic field power spec- 
trum with the help of Wick's theorem as explained bv lDurrer et a f (|2000l ). For 
example, for the magnetic energy density spectrum we obtain 

PpB{k) ^J^Y J d3 1 P B(<l)PB(\k-q\). (123) 

Note the difference between this power spectrum which describes fluctuations 
in the magnetic energy density and the mean magnetic energy density which 
is its zero mode given by the dfc/fc-intcgral of Ps(k)k 3 . 
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The metric fluctuations from this source of energy and momentum are 
calculated via the first order Einstein equations and they enter the Boltz- 
mann equation of CMB fluctuations as a source term. For vector perturba- 
tions, which are absent in sta ndard cosmology, they hav e been imple mented 
in CAMBcode (developed bv lLewis and Bridle! l2002h by iLewisi (|2004 ). Even 
though the magnetic field is a vector field, its energy momentum tensor (11221) 
which is quadratic in the field, contains scalar, vector and tensor contributions 
of similar amplitude. 

For a magnetic field spectrum behaving like k" s on large scales, k < ks, the 
spectrum of the magnetic energy momentum tensor is dominated by the upper 
cutoff and behaves as white noise on large scales, k < ks if n s > —3/2. For 
n s < —3/2 the energy momentum tensor inherits the magnetic field spectral 
index n s . 

A magnetic field with super horizon scale correlations, A ~ @h(z), keeps 
the universe homo g eneous but renders i t anisotropic, a Bianchi I model, see 
Barrow et al (1997); Adamek et all ( 2011 ). At temperatures below lMeV, where 



neutrinos free stream, such a global anisotropy is however compensated by the 
induced neutrino anisotropic stress gene rated by the gravita tional effects of 
anisotropic relativistic free streaming, see Adamek et all ( 20 111) . For a stochas- 
tic magnet ic field this leads to a, supp ression by a factor (k/H) 2 on super hori- 
zon scales ( Bonvin and CaprimfeOlOl ) . This effect is very relevant especially for 
scale invariant magnetic field spectra. However, such spectra which can only 
come from inflation also generate, in addition to the compensated mode and 
to the passive mode present in causally gen erated magnetic fields, e.g, from a 
phase transition (see lShaw and Lewisll2010l) . a small constant mode which is 
not co mpensated and which can have a scale invariant spectrum (|Bonvin et al 

In addition to gravitational effects which dominate on large scales, inter- 
actions of electrons with the magnetic field and the CMB affects the CMB 
anisotropy spectrum also via a modification of the acoustic peaks. The sound 
speed in the presence of a magnetic field is enhanced, c 2 — >• c 2 + (k • B) 2 /p, 
which affects the position of the acoustic peaks in the CMB anisotropy and 
polarization. Furthermore, Alfven waves, i.e. vecto r perturbtio ns are gener- 
ated which lead to relatively strong B-polarisation (|Lewisl 120041 ) and the Silk 
damping tail is affected. 

For a given magnetic field spectrum, all these effects can b e incl uded in 
standard n Boltzmann codes like CAMB by iLewis and Bridle ( 2002 ). which 
treat the photon+baryon+dark matter system within linear perturbation the - 
ory to compute CMB anisotropies and polarization ( Shaw and Lewis! 2012 ). 
The magnetic field contributions just enter via a source term on the right hand 
side of the linear perturbation equations and by modifying the initial condi- 
tions. They lead to a so called 'passive mode' which has the same initial condi- 
tions as the inflationary mode. In this mode the magnetic field just intervenes 
by changing the evolution equation of the baryons due to the Lorentz force. In 
addition there is a 'compensated mode' where the gravitational effects of the 
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Fig. 6 Constraints on the cosmologically produced IGMF from CMB anisotropy measure- 
ments. 



magnetic field enter. On very large, super horizon scales, these are, however, 
compensated by initial fluid under densities and neutr ino magnetic stresses, a t 
least if the magnetic fields a re generated causally . See IShaw and Lewis! (|2010ft : 
Bonvin and Caprinil ( 2010h ; Adamek et all ( 2011 ) . If the magnetic field is gen- 
erated during inflation, the compensation mechanism after neutrino decoupling 
is still active, but an additional 'passi ve mode' due to the matching condition 
at the end of inflation is introduced ( Bonvin et all 2012 ). The relation of its 
amplitude to the late time magnetic field strength depends on the details of 
reheating. 

The non-observation of the large angular scale anisotropies of the CMB 
have led to an upper limit B < 4 x 10 -9 G for the fields with correlation 



scale of the order of the CMB scale, \ B ~ 10/i _1 Mpc (jBarrow et allll997l) 



Limits dependin g on the power law of the magnetic field spectrum have been 
derived in Refs. JPurrer et afeood: lGiovaiminil2009l IPaoletti and FinellillioTl ; 
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Yamazaki et all 2012 ). The envelope of the upper bounds on (B , Xb) for the 
range of the power law indices —3 < n s < 2 (jPaoletti and Finellill201 lh is the 
lower boundary of the light blue shaded region of (B\,Xb) parameter space 
in Fig. |H1 

Similar limits can also be derived from the absence of non-Gaussianitics 
in the observed CMB anisotropies. Non-Gaussianitics from magnetic fields are 
mainly of the local type. The limits on the CMB bi-spectrum, usually param- 
eterized in terms of f n i, therefore constrain a possible contribution to CMB 

anisotropics and polarization from magnetic fields. These co nstraints can be 

used to limit the magnetic field amplitude for a given spectrum ( Seshadri and Subramanianl 
20091: ICaprini et alll2009bllTrivedi et ajl201oT) . Recently also the tri-spectrum 
has been calculated (jTrivedi et alll2012l) . 



5.3.2 Spectral distortions 



Non-thermal dissipation of magnetic field energy into the energy distribution 
of electrons befo re recombination can lead to distortions of the blackbody 
CMB spectrum ( Jedamzik et all 19981) . This distort ion generates a non-z ero 
chemical potential /i, which has been calculated bv I Jedamzik et a j (l2000h in 
the form of a double integral. In order to obtain bounds on B, Xb, we consider 
two limiting cases in which the integral can be performed analytically, namely, 
the cases when the magnetic field correlation length is much smaller or much 
larger than the characteristic damping length scale: 



Xd = 



to 



2tt 



? 3 / 2 V 15n° e a T 



400 pc, 



(124) 



where time constant is in = 2.4 x 10 1 



is the characteristic redshift 



to = /.i x iu~~ s, z M 
of freeze-out from double- Compton scattering z M = 2.5 x 10 6 , n° e is electron 
density and <tt is Thomson cross section. 

Taking into account the co nstraint \n\ < 9 x 10 -5 at 95 % confidence 
level from COBE FIRAS data (jFixsen et allll996lh one can derive analytical 
limit s on magnetic fields i n the two limiting cases. In the case Ab < Ac one 
has ( Jedamzik et al 2000l) : 



B < 3.2 x 10~ 8 G 



1 / Ai 



K V 400 P c 



-(n a +3)/2 



(125) 



where K = 1.4T(n s /2 + 5/2)r(3n s /5 + 9/5)2-(" s+5 '/ 2 (6/5)(n s + 3) is a con- 
stant of order unity, K — 0.8 resp. 2.1 for n s — —2 rcsp. +1. In the opposite 
regime, Xb 3> Ad, the constraint becomes 



B < 3.2 x 10 _B G- 



1 



A; 



K 2 \400 pc 



(126) 



where Ki is another constant of order unity. The strongest limit on the field 
strength is obtained for the fields with correlation length Xb ~ Xd. For such 
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fields no convenient analytical approximation can b e found and instead a nu- 
merical integration of the expression given in Ref. (jJedamzik et al has 
to be performed. In Fig. [6] we show the bound on B, Xb implied by the anal- 
ysis of distortions of the CMB spectrum as extrapolations of the analytical 
approximations given by Eq. (|125p . (|126[) for the entire ranges Xb < Xd and 
Xb > Arj. Note that the dependence of the limits on B,Xb on the power-law 
index n s practically disappears in the case Xb ^S> Xd (only weak dependence 
remains in the constant K%). For A < A^>, we consider the softest power spec- 
trum slope n s = 2 corresponding to the fields causally produced at the phase 
transitions in the Early Universe (see Section [3]). 

Apart from producing a non-zero chemical potential, transfer of the mag- 
netic field energy to electrons/ positrons can result in non-zero Compton pa- 
rameter y. Taking into account restrictions y < 1.5 x 10~ 5 from COBE FIRAS , 
one finds a limit B < 3 x 10~ 8 G at A - 0.3 - 0.6 Mpc (jjedamzik et alll2000l) . 
Note, that limits Eq. (|125[) and Eq. (|126p constrain magnetic fields created at 
z > z M = 2.5 x 10 6 , while the limit following from restrictions on y applies 
for fields created before z > z e i ~ 2 x 10 4 where elastic Thomp son scattering 
drops out of thermal equilibrium, i.e. t e i > t, see iDurrer (<2008l) . 



5.3.3 Limits from CMB polarization 



Since Thompson scattering is anisotropic, temperature anisotropics in the 
CMB (more precisely the quadrupolc) induce a small net polarization of CMB 
photons. Depending on the polarization pattern this is called E-polarization 
(gradient field on the CM B sky) or B-po larization (rotational pattern on the 
CMB sky). For details see Durrer ( 2008t ). Scalar perturbations from inflation 
only generate E-polarization. Magnetic fields lea d to strong vector pertur- 
bations which gen erate significant B-polarization (jSeshadri and Subramanian 
20011: lLewisll2004 . 

Magnetic fields which the CMB photons encounter on their way from the 
last scattering surface into our antennas Faraday rotate this polarization pat- 
tern and thereby generate also B-polarization. Due to its dependence on the 
wavelength given in Eq. (|114[) . this effect can be easily distinguished from 



achromatic gravitational effects, e.g. primordial gravitational waves which also 
generate B-polarization. Constraints o n B, Xb stemming from non -observation 
of this effect were first discussed by iKosowskv and Lqeb and subse- 

quently updated using the 5-years data of WMAP by iHinshaw et all (|2009h . 
The limits coming from non-observation of Faraday rotation in CMB signal 
are, at present, weaker than the limits imposed by the rotatio n measures of dis- 
tant b lazars or limits from the CMB angular power spectrum (jKahniashvili et al 
2009h . 



5.3.4 Limits from magnetic field effects on LSS 



Cosmological magnetic fields present at the onset of LSS formation may af- 
fect the properties and the statistics of the matter perturbations. In this 
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way, strong enough field might affect the abundance galaxy clusters and 
there fore, modify the Sunyaev-Zeld ovich effect on CMB angular power spec- 
trum rtTashiro and Sugiyamal 12011 ). Using the data of the South Pole Tele 



scope, [Sh aw and Lewis! (|2012|) derived a limit B < 4 nG at Xb = 1 Mpc scale 
In a similar way, the presence of strong magnetic fields can distort the statis- 
tics of column densities of the Lya clouds. Non-observation of this effect also 
constrains the field strength. Combining the CMB angular power spectrum, 
Sunyaev-Zeldovich effect data, and L ya data, a joind constra int B < 1 nG at 
Mpc distance scale has be derived by IShaw and Lewis! (|2012j) . 

Strong magnetic fields provide an additional contribution to the pressure 
and thereby enhance the Jeans mass for gravitational collapse. The influ- 
ence of magnetic fields can be strong enough to partially suppress the col- 
lapse of small structures. This leads to the decrease of the amount of ioniz- 
ing radiation produced by the star formation in small galaxies. On the other 
hand, the Lorentz force induces additional perturbations in the baryons which 
are transferred to dark matter and enhance collapse of small scale structure 
above the magnetic Jeans mass dSethi and Subramanianll2005l ; iKim et alll996t 
Subramanian and Barrowlll998a). This can also i nfluen ce t he onset of reioniza- 
tion of the Universe, see lSethi and Subramanian and lTashiro and Sugivamal 

( 20061 ) . Based on this idea, Schleicher and Miniatil ( 2011 ) derive a scale- independent 
upper limit B < 3 nG on cosmological magnetic fields present at the redshift 
z ~ 7. This limit is comparable to the limits imposed by the CMB data. 

The effects of magnetic fields on structure formation and the thermal and 
ionization evolution of cosmic hydrogen also leave traces in t he 21cm signal 



2009h : ISchleicher et all 



which have been investigated in lSethi and Subramanian! ( 
(2009). It is found that future radio surveys like SKA could probe magnetic 
fields with correlation scales of O.lMpc to several Mpc down to amplitudes of 
10- 10 Gauss. 

The modification of LSS formation in the prese nce on magnetic fie l ds als o 
affects gravitational weal lensing, as discussed in IPandev and Sethil (|2012f ) . 
Again constraints of a few nGauss are obtained for a spectral index n s ~ —2.9 
of the magnetic field spectrum. 



5.3.5 Helical fields and the CMB 

If magnetic fields are helical, their induced gravitational fields arc so as well 
which leads to non-vanishing cross correlation spe ctra of temperatur e and 
B-polarization and wall as E- and B-polarization ( Caprini et all 2004 ). The 



non-detection of such cross-correlations can be used to place additional limits 
on helical magnetic fields. However, they are also of the order on nG as the 
other limits discussed above. 
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5.4 Limits from gravitational waves 

The energy momentum of magnetic fields always gives raise to anisotropic 
stresses. Their transverse traceless component are sources gravitational waves. 
More precisely 

hij + ZHfkj + k 2 h l3 = IGnGa^niP . (127) 

( B) 

Here h,j is the gravitational wave amplitude and 7I 8 - • is the transverse trace- 
less part of the anisotropic stress due to the magnetic field. On a scale A = 
2ir/k ~ "H" 1 within a Hubble tim e, gravitational waves with energy den- 



sity (|Caprini and Durrerll200lT |2006( ) 



*^ = «f(^) (128, 

^ ^rad \ J 'rad / 

are generated. The dominant contribution to the gravitational waves on scale 
A are generated when this scale crosses the horizon so that we obtain 



^rad \ ^rad 

Once generated, these gravitational waves propagate freely and do not inter- 
act anymore. Even when the magnetic field on these scales is damped away 
by fluid viscosity, the gravitational waves which it has produced remain. They 
contribute an additional relativistic component to the expansion of the Uni- 
verse. The strongest limits on t his come from combining nucleosyn thesis con- 
straints and CMB observations ( Steigmanll2008t iHinshaw et all 2012 ) , requiring 



J?GwV^rad ~ 0-1 on an scales. As I2s(t)//? ra d decays during free turbulent 
evolution, the strongest constraints are obtained on scales of the order of the 
horizon scale at formation, A ~ t* but also there the above constraint just 
requires i?g(t) / 'i2 ra d ~ 0.3. 

Of course for a specified process of magnetic field generation with a fixed 
spectrum, t his can imply much s tronger limits on much larger scales, for an 
overview see ICaprini et a 3 (l2009ah . For example, if the formation is causal and 



evolution proceeds via incompressible MHD turbulence, the magnetic field 
energy on a given scale k at late times is always smaller that its unprocessed 
value which behaves like 

*£^~h?P B (k) = B 2 (k)ock*, (130) 
a log k 

so that on the scale A = 2n/k we obtain the limit 

^^T £(fc/fc * )5 " (fc/10 °^ )5 - (131) 

Here we have assumed an initial correlation scale of 1% of the initial Hubble 
scale. This is probably a reasonable value for first order phase transitions 
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which proceed vial bubble nucleation (see Huber and KonstandiiJ 2008 ). For 



an arbitrary initial correlation scale k* , we have 



k 



5/2 



B(k) ~ B(k*) [ — ) . (132) 



This provides a very strong limit on large scales. We use 

*<fl lW1 i-u(5^) , <0i. (133) 

If this limit is satisfied at the initial correlations scale A* for magnetic fields 
generated at the electroweak phase transition, if* ~ 10 _2 mHz~ (10~ 9 Mpc) _1 , 
on the scale A ~ lOOkpc a magnetic field of the order of 

B(X ~ lOOkpc) ~ 10~ 31 G . (134) 

If we allow for compressible MHD evolution, with a scalar velocity com- 
ponent of similar amplitude as the vorticity, this constraint is significantly 
relaxed. Now the correlation scale increases like A = A»(i/i*) 2 / 5 , and the mag- 
netic power spectrum remains roughly constant at the peak, -Pb(A) = -Pb(A»), 
see Fig. [5J Setting the endpoint of the free turbulent decay at recombination, 
Tfin ~ 0.3eV, we obtain a correlation scale which is about 5 orders of magni- 
tude larger, 

A fin ~ (100GeV/0.3cV) 2/5 10 2 scc ~ lpc. 

The magnetic field power spectrum on this scale has the same amplitude as 
Pb(A*) at the electroweak phase transition. At larger scales it decays like A -2 . 
Inserting this in the condition (|131[) yields 

B{\ ~ lOOkpc) ~ [P S (A)(A/A) 2 A- 3 ] 1/2 = B{K){XJ\f' 2 {X/X r ) . (135) 

Hence the limit is relaxed by about 5 orders of magnitude with respect to the 
one coming from incompressible turbulence. If the magnetic field is maximally 
helical, the limit is relaxed even more. These limits can also been obtained from 
the corresponding endpoints in Fig. [3] by taking into account that on scales 
larger than the correlation scale, A > Xb = A the magnetic field strength of 
causally generated fields decays like (As/A) 5 / 2 . 



5.5 Potential limits from ultra-high-energy cosmic ray observations 

Magnetic fields in IGM can be probed by measuring their effect on trajectories 
of Ultra- High-Energy Cosmic Rays (UHECR), which are charged particles 
with energies Stjhecr ~ 10 20 eV penetrating in the Earth atmosphere. Such 
particles come to the Earth from yet unknown sources which are most probably 
situated outside the Milk Way. This means that UHECR traveling from their 
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source to the Earth cross the IGM. A magnetic field in the IGM deflects 
UHECR trajectories from their straight line by an angle 



9lGMF = * Z [l^eVj [w^g) [soWc) ' (136) 

where Ze is the electric charge of UHECR particle, D is the distance to the 
UHECR source and B± is the strength of IGMF component orthogonal to 
the line of sight. In the above equation the IGMF is assumed to be coherent 
over the length scales larger than the distance D. If the coherence length 
As is much shorter than D, UHECR experiences a sequence of deflections in 
different random directions during the passage through each distance interval 
As, so that the overall deflection angle accumulated over the path D is given 

by 



ZeB ± y/DX B / -guHgciA n „ 7 , 

W " ^uhecr " °' 4 Z {m^eVj (137) 
B ± \ ( D \ 1/2 / \ B 



10~ 10 GJ \50 Mpc/ V 1M P C 



IGMF affect the arrival directions of UHECR particles by displacing them 
from the direction towards the source. This opens a possibility of the mea- 
surement of I GMF with the strength in the rang e B ~ 10~ 9 G using UH ECR 
observations ( Lee et all 1995 ; Lemoine et all 1997 ; Sigl and Lemoine 1998 ) . In- 
deed, if the sources of UHECR would be known, one would be able to measure 
the angular distribution of UHECR events around the source positions. In fact, 
the mere observation of isolated sources of UHECR on the sky would indicate 
that #igmf "C 1, so that IGMF can not be stronger than ~ 10~ 9 G, provided 
that the distance to the UHECR sources is in the D ~ 50 Mpc range. 

Recent attempts to model the deflection of UHECR by the magnetic field 
of the intervening large sca le structu r e, suc h as g alaxy clusters and /or fil- 
aments by two groups (see Sigl et al ( 2004 ) and Dolag et al ( 20051) fled to 
contradictory results, which reflect uncertainties of the structure of magnetic 
fields inside and around clusters and filaments. In addition, if there is a sin- 
gle nearby UHECR source, the host galaxy or galaxy cluster of the source can 
span several degrees on the sky. Significant deflections of UHECR by magnetic 
fields in the host galaxy or galaxy cluster can then produce ext ensions of the 
UHECR emitting region of 1 up tolO degrees (jDolag et alll2009l) . 

There are several obstacles for the measurement of IGMF with this method. 
The first is the uncertainty of the origin of UHECR and of their composition. 
The second is the presence of additional deflection of UHECR trajectories by 
the Galactic magnetic field. Below we summarize the status of these uncer- 
tainties. 
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5.5.1 Sources and the composition of UHECR 



The main problem for the localization and identification of the UHECR sources 
is the low statistics of the signal in the 10 20 eV energy band. Starting from 
the earliest days of UHECR observations and up to the present only about 
~ 10 2 events hav e been collected by different experiments: Pierre Auger Ob- 
servatory (PAO) (Pierre Auger Collaboration et al 2007t Abraham et al 2008: 



The Pierre Auger Collaboration et a ll2010l). HiRes (I Abbasi et aj2008ll2010al 
AGASA (jHavashida et alllf996l) . Telescope Array (TA) (|Abu-Zavvad et all2012 



EL, 

L23Q). 



Abbasi et all ()2010a|) found evidence for a proton-dominated cosmic ray flux 
at the energies above 10 18 eV, based on the analysis of the average depth of the 
maxima of Extensive Air Showers (EAS) of high-energy particles, produced 
by the UHECR penetrating in the Earth atmospher e. If the proton-dom inated 
composition persists until the 10 20 eV energy band (| Abbasi et alll2005l ). a sup- 
pression of the UHECR flux due to the interactions of high-energy protons 
wit h the CMB photons, known as the Greisen- Zatsepin-Kuzmin (GZK) cut- 
off (jGreisenl Il966t IZatsepin and Kuz'mirj 119661 ) is expected be yond this en- 
ergy. This suppression was first observed in the UHECR flux ( Abbasi et all 
2008) in the HiRes da ta and lat er confirmed with better statistics by the PAO 
(jAbraham et alll2008t ) and TA (|Abu-Zavvad et al|[2012bl) experiments. 

An indication for clustering of UHEC R arrival directions on small angu- 
lar scales was first found in AGASA data ( Havashida et all 1996 ). A stronger 
small-scal e anisotropy signal was subsequently found in Peirre Auger observa- 
tory data ( Pierre Auger Collaboration et alll2007 , 2008 ). The arrival directions 
of UHECR events were found to correlate with the sky positions of nearby 
Active Galactic Nuclei (AGN) situated within a distance D < 75 Mpc. The 
correlation was found at angular scale 9 ~ 1° equal to the angular resolu- 
tion of the detector. In the initial analysis ( Pierre Auger Collaboration et all 
20071 ) a significant fraction o f UHECR events was f ound to correlate with the 
AGN positions (see, however iGorbunov et all (|2008l) ). However, the number of 
sources contributing to the signal and, as a consequence, the significance of 
the UHE CR - AGN correlation, has decreased with th e accumulation of event 
statistics ( The Pierre Auger Collaboration et al 2010l) . No significant correla- 
tion with AGN was found in the analysis of the HiRes data with com parable 
event statistics ( High Resolution Flv'S Eve Collaboration et all 2008 ) and in 
the Telescope Array data (|Abu-Zavvad et all20I2a ). 

Detection of small angular scale correlations of UHECR arrival directions 
with particular point sources on the sky would immediately imply that the 
UHECR trajectories are not strongly deflected by the IGMF. Thus, if the 
PAO result on the correlation with AGN holds, it immediately implies that 
the UHECR are protons (Z = 1 in Eqs. 11361 IT371 and that a I GMF can not 
be stronger than - lO" 10 G for X B ~ D (|de Angelis et alll2008l) . 

However, the interpretation that most of the UHECR are protons appears 
to be in confl ict with a recent find ing of heavy composition of the highest 
energy events. Abraham et all ( 2010l) have concluded from the decrease of flue- 
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tuations of the depths of EAS maxima at energies above 10 19 eV that these 
UHECR events are composed of heavier nuclei. A possible heavy composition 
of UHECR flux has to be verified with better statistics and using independent 
analyses fro m the HiRes / TA d etector, which does not confirm this result at 
the moment (jAbbasi et alll2005l ). 

The absence of firm source identifications and the ambiguities in the com- 
position of UHECR at the highest energy end render the use of the UHECR 
techniques for the measurement of the IGMF very uncertain at the moment. 
A dramatic increase of the UHECR event statistics is required to remove the 
above mentioned uncertainties. This requires a dramatic increase (by at least 
an order of magnitude) of the collecting area of the UHECR detector. Existing 
ground-based detectors, like PAO and TA achieve collecting areas in the range 
of (1 — 3) x 10 3 km 2 . A further increase of the area is challenging because this 
would require homogeneous coverage of large surfaces with high-energy par- 
ticle detectors, over the areas in the range comparable to the size of a small 
country like Switzerland. A particular UHECR detection technique, based on 
the measurement of ultraviolet fluorescence emission from the EAS, allows for 
a qualitatively different space-based approach: deployment of a single wide 
ficld-of-vicw downward looking telescope in space. This approach allows to 
monitor a large part of the Earth surface and potentially reach a collect- 
ing area comparable to the Earth surface. A first attempt along this direction 
will be the planned next-gen eration space-based UHECR detector JEM- EUSO 
http://iemeuso.riken. ip (|The JEM-EUSO Collaboration et alll2012l) . It will 
have a collecting area in the range of 10 5 km 2 and an all-sky exposure capabil- 
ity. JEM-EUSO will collect some ~ 10 3 UHECR events in a ~ 5 yr exposure 
on board of the International Space Station. Possible usefulness of UHECR 
observations for constraining IGMF has to be re-assessed when the results of 
the all-sky small- and large-scale anisotropy analysis of JEM-EUSO UHECR 
data will be available. 



5.5.2 UHECR deflections in the Galactic magnetic field 

Apart from uncertainties of the sources and composition of UHECR, another 
major obstacle formcasuring IGMF with UHECR is related to the uncertain- 
ties of the structure of the Galactic magnetic field (see also Section I5T2"]) . The 
magnetic field of the Milky Way is conventionally modeled as the sum of a 
regular and a turbulent component of the field in the disk and in the halo of 
the Galaxy. Therefore, the deflections of UHECR by the Galactic magnetic 
field can be decomposed onto four terms: 

n /^regular . ziturbulont , /jrcgular . ziturbulent finol 

VGal - ^Disk + 6/ Disk + 6/ Halo + 6/ Halo . [L6ti) 

Deflection by the regular and turbulent components of Galactic disk and halo 
can be estimated by substituting the typical disk/halo size at the place of D 
and the typical disk/halo field strength and correlation lengths at the place of 
B and X B in Eqs. (flip]) , (fT37jl . Deflections of UHECR by the regular field in the 
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disk #nS lar w ere studied in many theoretical models starting with a paper by 
IStanevI ( 1997t h Typical values of parameters entering the analog of Eq. (|136|) 
imply -DDisk rcg — 2 kpc (for sources located far from the Galactic Plane) and 
-Boiskrog — 2 /LtG, which give for Eq. (|T36| D 3 g£ lar ~ 4°. Turbulent fields are 
typically assumed to have a coherence scale of ~ 50 pc and B ~ 4 /j,G, which, 
for the same scale -DDisk turn ~ 2 kpc gives # D U isk Ulent — 0.5°. Contributions 
of the Halo fields are less certain, but result in deflections which ar e at least 
of the same order, see the recent discus s ions of all compon ents by Sun et al 
(|2008l ); iHanl (|2009l) ; Ijansson et all (|2009t ) ; iFarrar et all (|2012f ). 

Although the order of magnitude of the deflection angle of UHECR by the 
Galactic magnetic field can be readily estimated, uncertainties in the mea- 
surements of Galactic magnetic field and discr epancies between the existing 
measurements and existing theo retical models ( Sun et al 2008 ; Jansson et all 
20091 : Jansson and Farrar 2012al lb) do not allow to predict the deflection an- 



gle and direction of deflection for particular lines of sight (toward UHECR 
sources). This means that, most probably, the details of the structure of the 
Galactic magnetic field along the line of sight toward UHECR sources will 
have to be deduced from the UHECR data itself, rather than just taken into 
accou nt in the UHECR data analysis ( Sigl and Lemoinelll998t iGiacinti et al 



20f0l ). This, obviously, will introduce large uncertainties into the derivation of 



the properties of IGMF from the UHECR data. 

Deflections of UHECR by the Galactic magnetic field are much stronger 
if UHECR particles are heavy nuclei, rather than protons, as indicated by 
the recent PAO results ( Abraham et al 2010l ). In this case th e UHECR are 
distri buted around their source on angular scales of 10° — 100° ( Giacinti et all 



20111 ). Inhomogeneities of the exposure across the fields-of-view of ground- 



based detectors is a significant obstacle for the recognition of such wide-angle 
patterns of UHECR distribution around their sources. Only all-sky detectors, 
like JEM-EUSO would be potentially able to localize the UHECR sources in 
this case. It is clear that more sophisticated procedure for the identification of 
UHECR sources for heavy UHECR will automatically imply lower sensitivity 
of UHECR for the measurement of IGMF. 



5.6 Constraints from gamma-ray observations 

An alternative way to probe the weakest magnetic fields in the IGM is using 
high-energy (HE, photon energies in the 0.1-100 GeV range) and very-high- 
encrgy (VHE, photon energies in the 0.1-10 TcV range) 7-ray observations. 7- 
rays propagating though the IGM occasionally interact with the abundant low 
energy photons and produce electron-positron pair. Electrons and positrons, 
being charged particles, are affected by the magnetic field. Therefore, sec- 
ondary 7-rays produced by electrons and positrons interacting with CMB 
photons carry information about the properties of IGMF. 
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5.6.1 Absorption of 7 -rays in the IGM 



To produce e + e~ pairs, the energy of the primary 7-ray must be high-enough 
so that the center-of-mass energy in the photon-photon collision exceeds twice 
the electron mass. 7-rays with energy i? 7 propagating through a background 
of soft photon s with energy e can produc e pairs if their energy is higher than 
the threshold ( Gould and Schreder 19661) 



E 1 > 



250 



1 eV 



GcV. 



(139) 



Soft photons with energies in the 0.1 — 10 eV range are abundant in the Uni- 
verse, because they are produced by star formation in galaxies. The homoge- 
neous and isotropic soft photon back ground in this energy range, known as Ex- 
tragalactic Background Light (EBL) ( Madau and Pozzettil2000t Hauser and Dwek 



2001 ; IFranceschini et al 2008; D wek and Krennrichl2012l ). is accumulated through- 
out the history of the Universe, starting from the onset of star formation at 
redshift z ~ 10. The spectral energy density of EBL has a characteristic two- 
bump shape with a near-infrared bump at ~ 1 eV being due to the direct 
starlight emission and the far-infrared bump at the energy ~ 10 -2 eV being 
produced by scattering of starlight with d ust, see Fig. [7] The energy density o f 
the EBL is about p E BL ~ 1CT 2 ' 5 cV/cm 3 (jH.E.S.S. Collaboration et alll2012al) . 
which is a factor of ~ 10 2 lower than the energy density of the CMB. 

The precise value of pebl was somewhat uncertain (by a factor of ~ 2) un- 
til recently, because it is not possible to measure it in the optical and infrared 
bands. The main obstacle for such a measure ment is the much stronger fore- 

et all 
20021 ). A 

proper subtraction of the zodiacal light foreground requires a precise knowl- 
edge of its spectral characteristics and angular distribution in the sky. Un- 
certainties in these parameters introduce an uncertainty in estimates of the 
spectral characteristics of EBL based on the measurements of diffuse emission 
in the optica l and infrared band s. This uncertainty has affected early estimates 
of the EBL ([Hauser et allll998T ). which have turned out to be too high. 

Alternatively, the amount of soft photons in the IGM can be estimated 
from source counts in deep observations by infrared and optical telescopes. 
This approach provides lower bounds on the EBL density (since the contri- 
bution from undetected sources with fluxes below the telescope's sensitivity 
is not taken into account). Along these lines, tight lower bounds on the EBL 
density at different wavelengths wer e derived from the observations with tele 



ground of zodiacal light from the Solar system ( Puget et~al 19961: Hauser 1 
19981 ; IWright and Reesel l2000t iGoriian et a] |2000|; iBernstein et all 12002 



scopes in space and on the ground (lElbaz et all 2002 ; Matsumoto et al 
Fraver et aj 120061: iDole et all 120061 Kashlinskv et all l2012t iBerta et all _ 
Oliver et all 2010l : Keenan et all 201C ) . Dwek and Krennrich (|2012l) have sum 



2005 



2010; 



marized the constraints on the EBL spectrum obtained from the direct source 
count observations. These constraints are shown by the grey-shaded region in 

Fig.m 
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Fig. 7 Measurements, limits and model calculations of the spec tral energy distribution o f 
the EBL. The black solid curve shows the model calculation by [Franceschini et a] (120081 ') ■ 
The lue dotted curve sh ows the calculat ion of lDommguez et all ll201lf) . and the green dashed 
line shows the model b\ lKneiskei J2008t). The grey s haded region is the envelope of direct 
measurements summarized bv lDwek"* and Krennrichl i2012h. The blue shaded region is the 
measurement derived from 7-ray observations bv lH.E.S.S. Collaboration et all [|2012ah . The 
orange curve shows the CMB spectrum. 



Interactions of Very-High-Energy (VHE) 7-rays (photon energies above 
100 GeV) with the EBL photons can be used to constrain or measure the 
EBL density and spectrum. These interactions lead to an energy-dependent 
suppression of the 7-ray flux from extragalactic sources in the VHE band. 
Indeed, the mean free path of 7-rays with energy above the threshold (|139[) is 
given by 



A 77 (.E 7 ) — — 

Jm 2 IE 0-77(^7, e )^EBL(ejae cr 77 n E BL 
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0.8 



E-y 


-l 


Pebl 


-l 


1 TeV 




10 -2.5 eV/cm 3 





Gpc . 



(140) 



Here ^ebl(c) is the number density of EBL photons, tiebl — Pebl/c, cr 77 (e) is 
the pair production cross-section whic h reaches it maxi mum ct 77 ~ 10 -25 cm 2 
at aout 4 times the threshold energy (|Aharomanll2004l ) and we have to insert 
the value of pebl at this maximum, e max ~ 4m 2 /_E 7 . The mean free path of 
7-rays becomes shorter or comparable to typical distances D to extragalactic 
7-ray sources in the TeV band, see Fig. [51 so that the source fluxes are sup- 
pressed by a factor exp(— D/A 77 ). If the intrinsic source spectrum is known, 
comparison of the observed and intrinsic spectrum of the source provides a 
measure of the suppression factor and, via Eq. Q140p , of the density of the 
EBL. 

The main class of extragalactic sources of VHE 7-rays are blazars. These 
are radio-loud A GN (active galact ic nuclei) emitting jets closely aligned long 
the line of sight ( Aharonianl 20041) . Recent rapid development of the ground- 



based VHE 7-ray astronomy has resulted in the dis covery of about 10 2 extra- 



galactic VHE 7-ray sources ( Aharonian et all 2008 ). In a significant fraction 
of these sources, a high-energy suppression of the flux due to the interactions 
with EBL photons is clearly observed. 

In principle, if the quality of these data is sufficiently good, one can derive 
both, the intrinsic shape of the 7-ray spectrum and the energy-dependent 
suppression factor, directly from the 7-ray data. For this one has to use the 
fact that the suppression factor exp(— D/A 77 ) has a characteristic "two bump" 
shape (see Fig. [5]), related to the two components in the EBL spectrum. Fig. 
[5] shows the energ y dependence o f A 77 determined with Eq. p40[) from the 
EBL spectrum by iKneiskd ( 2008 ). which is consistent with all observational 
constraints (see Fig. [7]). The two enhanced absorption energy intervals at E 1 ~ 
1 TeV (absorption on the near-infrared EBL)and E 1 ~ 100 TeV (absorption 
on far-infrared EBL) are clearly seen as deviations of the estimate (|140p from 
a straight line. The shape of the enhancement of absorption around 1 TeV is 
determined by the shape of the EBL spectrum in the near-infrared band. Since 
A"^ is proportional to the logarithm of the suppression factor exp(— D/A 77 ). 
Fig. [8] shows in fact the shape of the EBL suppression factor, which can, in 

principle, be deduced from the 7-ray data. 

This approach has been followed recently using HESS telescope date ( H.E.S.S. Collaboration et all 
The spectrum of the EBL derived in this way from the 7-ray data 



2012al ) 



is shown in Fig. [7] by the light-blue hatched band. A similar analysis using 
data of the Fermi telescope at lower ene rgies has allowed determ ine the level 
of the EBL at higher redshift z ~ 1 ( Ackermann et al 20121 ). Earlier cal- 
culations, based on lower quality data, had to rely on assumptions about a 
range of possible shapes for the intrinsic blazar spectra in the TeV energy 
band. This had then permitted to derive upper limits on the EBL density , 



Funk et all 


1998: 


Aharonian et all 


2006. 
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20081: iMazin and Raud2007tlFinke and Razzaauel2009l:lH.E.S.S. Collaboration et all 
2012bl) . 

The direct measurements of the EBL and the determination using 7-ray 
data are supplemented by models deriving the EBL spectrum and its evolution 
with rcdshift from models of cosmological evolution of opti cal and infrared lu- 
minosities and and of the luminosity function of galaxies (Stecker et all 2006: 
Franceschini et alll2008l : iKneiskeHiooa IPommguez et afcOllUFinke et allhoioi 
Gilmore et all 20121 ) . These models have evolved with time, starting from early 
models which attempted to explain the relatively high level of the EBL de- 
rived from DIRBE data to the currently favored models which predict a EBL 
flux at the level close to that derive d from the source count statistics, see 
Fig. [T] where example model spectra ( Franceschini et al 20081: Kneiskd 120081 
Dommguez et alll201ll ) are shown. 



5.6.2 Secondary j-rays from e + e pairs in the IGM. 



The mean free path of VHE 7-rays Q140p decreases monotonically with the 
energy E 1 . This is because the 7-rays of higher energies can interact with the 
more abundant lower energy photons from the EBL. From Fig. [5] we see that 
already for relatively nearby blazars, such as Mrk 421 (D ~ 150 Mpc), the 
mean free path of 7-rays with energy E 7 > 10 TeV is shorter than the distance 
to the source. 

The absorption of the highest energy 7-rays which due to the pair produc- 
tion on EBL photons injects electron-positron pairs in the IGM. The pairs are 
highly relativistic and are injected along the 7-ray beam from the source, at a 
distance of about A 77 . All the initial source power contained in the absorbed 
VHE 7-rays is transferred to the e + e~ pairs. . This power does not accumulate, 
because high-energy electrons and positrons rapidly loose energy via inverse 
Compton scattering of CMB photons. The dista nce scale on which relativistic 
pairs of energy E e ~ E 1 /2 loose their energy is dBlumenthal and Gould Il970h 



Ac = 



3mi 



0.3 



1 TeV/ 



Mpc, 



(141) 



where pcmb — 0.25 eV/cm 3 is the energy density of the CMB and ctt is the 
Thomson cross-section. This means that all power injected in the e + e~ pairs 
is finally converted into the power of inverse Compton 7-ray emission. The 
mean energy of the inverse Compton photons is 



Eke 



4ecMB-Ee 

3m? 



1 TeV/ 



GcV, 



(142) 



where £cmb — 3Tcmb is the mean energy of CMB photons. Thus, absorption 
of VHE 7-rays in the IGM finally resul ts in the generation of secondary lower 
energy 7-ray emission from the IGM ( Aharonian et all Il994h . This situation 
is illustrated in Fig. where the spectrum of the secondary 7-ray emission 
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Fig. 8 The mean free path of VHE 7-rays as a function of energy. The horizontal green 
line shows the Hubble radius. Dotted lines mark distances to some VHE 7-ray sources. The 
Inclined dashed blue line shows the estimate of Eq. 1 140ft . Labels along the curve mark the 
energies of soft photons with which 7-rays interact. 



generated along t he 7-ray beam from a distant blazar 1ES 0229+200 is shown 
( Vovk et al 2012 ) . One can see that absorption of the primary source flux in 
the TeV band leads to the production of the secondary (cascade) emission 
with spectral energy distribution stretching over several decades in energy, 
from E < 1 GcV to E > 100 GeV. The two panels of the Figure show cal- 
culations adopting different assumptions about the properties of the intrinsic 
7-ray spectrum of the blazar. In both cases the intrinsic spectrum has the form 
of a cut-off power law 



(143) 
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E 7 [eV] 



Fig. 9 The 7-ray spectrum of the blazar 1ES 0229+200, from lVovk et all (120121 s ). Thin solid 
lines show the intrinsic primary source spectrum. Thick solid lines shows the sum of the 
absorbed primary source spectrum and the secondary cascade spectrum after propagation 
of the 7-rays through the IGM from the source to the Earth. Dashed, dotted and dashed- 
dottcd lines show modifications of the propagated spectrum in the presence of IGMF with 
the strength given in the figure legends. 



with slope r and high-energy cut-off E cui . The model shown in the left panel 
assumes a softer slope, a power law with _T = 1.5, while the model in the right 
panel assumes a harder intrinsic spectrum, r = 1.2. 

The secondary 7-ray emission in the range 1-100 GeV is, in principle, de- 
tectable by telescopes sensitive to this energy band, such as Fermi ( Atwood et all 
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Fig. 10 Geometry of secondary 7-ray emissio n from electromagne tic cascade developing 
along the VHE 7-ray beam of the blazer. From iNeronov et all J2010f) . 



2009h and AGILE (ITavani et a which are currently in operation. Sec- 

ondary emission at the energies above 100 GeV can also be detected by the 
groun d-based 7-ray telescopes, HESS, MAGIC and VERITAS (jAharonian et all 
2008h . 

The observed spectrum in the 1-100 GeV band shown by the crosses in 
Fig. [5] contains two contributions: the intrinsic source spectrum in this band 
and the secondary cascade emission from e + e~ pairs generated by interactions 
of TcV photons with the EBL. If only the source spectrum is considered, we 
can in no way to distinguish between the direct and cascade contributions to 
the source flux in this energy band. However, the two contributions can be 
distinguished based on their different imaging and timing properties. 

Indeed, the secondary cascade emission is produced by an extended source 
in the IGM with the linear size given by the mean free path of the VHE 
7-rays (|140p . while the primary emission is produced by a point source at 
the location of the blazar. A telescope with sufficient angular resolution can 
in principle resolve the overall emission in the 1-100 GeV band into a point 
source + extended emission from the IGM. A schematic representation of this 
system is shown in Fig. [TUJ A typical VHE 7-ray source, a blazar, emits most 
of the 7-ray flux into a narrow conical jet with an opening angle (9j c t ~ 0.1°, 
as inferred directly from the observations of blazar jets in the radio band and 
from indir ect observations of relativ istic bulk motions in blazar jets in radio- 
to-7-rays ( Urrv and Padovani 19951 ). Blazar jets are almost aligned with the 
line of sight, so that the direction toward the observer lies within the 7-ray 
emission cone. In genera, however, the jet axis is slightly misaligned by an angle 
&obs < Ojet with the direction toward the observer, as shown in Fig. [TO] VHE 
7-rays interacting in the IGM deposit electron positron pairs in the interior 
of the shaded cone tracing the 7-ray beam in Fig [TO] Secondary cascade 
emission comes from the locations within this cone, rather than directly from 
the primary point source marked by the star in Fig. 1101 

From Fig. [10] it is clear that only electrons and positrons moving in the 
direction of observer would contribute to the secondary 7-rays emission from 
the shaded conical region. A simple geometrical argument tells that if elec- 
trons/positrons would move exactly in the direction of the primary 7-ray, 
no extended emission from the shaded conical region in Fig. [TU] would be 
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observable. However, the pair trajectories are always misaligned with the 
primary 7-ray direction. The minimal possible misalignment angle S is de- 
termined by the kinematics of the pair production collision, S > m e /E y ~ 
5 x 1CT 7 (£ 7 /1 TeV)" 1 . 

Let us now consider am magnetic field which permeates the IGM. In this 
case, electrons and positrons move along curved trajectories whose curvature 
radius Rl is determined by the strength of magnetic field 

E e o / E e \ { B 



R L = — ~10 2 [ — - — Mpc. (144) 

eB \l TeV/ VlO- 17 GJ v v ' 

Deflections of electron/positron trajectories by the magnetic fields lead to an 
additional misalignment of the particle trajectories with the primary 7-ray 
direction. After the propagation over an inverse Compton distance scale Die 
(|14ip . the misalignment angle is 



6 



R L - ' 2 {TTcV) ll(T 17 gJ V0.3 Mpc J ' IC B 



(145) 

For the second line above, we have taken into account many stochastic deflec- 
tions which led the electron or positron perform a random walk. 

A straightforward geometrical calculation based on Fig. IT U] leads to the 



estim ate of the angular size of extended emission from the IGM ( Neronov et all 
201dh 

©ext(-B) = min [S/t, S - 9 ohs ] < 6> oxt , ma x (146) 

where t = D/A 77 is the optical depth for the primary 7-rays. 

The difference in the path between the direct 7-ray signal reaching the ob- 
server from the primary source and the signal from the 7-rays which are con- 
verted to the e + e~ pairs in the IGM and then re-produced as inverse Compton 
emission leads not only to a displacement of the secondary emission from the 
primary source position but also to a time delay of the cascade signal. This 
is illustrated in Fig. [11] in which the results of Monte-Carlo simulations of 
the signal from the casc ade developing in th e IGM with a magnetic field of 



g m th e 

B = 1CP 16 G is shown ( Neronov et all l201Ctl . The characteristic t ime delay 



scales with the off-source angle 6> ox tis given by (INeronov et a 



D ( sin 6* C xt + sin(# bs + ®jet) \ D0 cx t{Q o \> s + (9j C t) 
tdoi(f) ^ — -5 —ET~\ 1 



sin(<9 cxt + 6» obs + ©jet) / 2c 
„| ° V<*t±9-)V«=lU. .147, 



150 Mpc J V 1° / V 1" 

The angular size of the extended emission from the IGM (|146p and, as a con- 
sequence, the time delay of the emission (|147[) depend on the strength and 
correlation length of the IGMF. Thus the IGMF parameters can be deduced 
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Fig. 11 Simulated 1-10 GeV band images of the sky region around a TeV blazar with a jet 
with ©jet = 9 bs = 3° at different times after the instantaneous injection of 1 TeV 7-rays at 
the source. The IGMF parameters are B = 10 -16 G, Xb = 1 Mpc. From left to right: time- 
integrated emission; images after a delay time t da \ < i c i c i < 10 5 yr, 10 5 yr< t<j e i < 10 6 yr, 
10 6 yr < t dol < 3 X 10 6 yr and 3 X 10 6 yr< f. dcl < 10 7 yr after an outburst. From lNeronov et all 
The grid spacing is 2°. The color code indicates the fraction of the total energy 
emitted in this region. 



from imaging (Neronov and Semikozll2007l 2009: Elviv et all 


2009) and /or tim- 


ing (Plaga 19951 iDai et all 20021: Neronov and Semikod 2009 


: Ichiki et al 20081: 


Takahashi et all 


2008) observations of the cascade sienal bv T-rav telescoues 



sensitive in the 1-100 GeV band. 

If the gyro-radius of the high-energy e + e~ pairs is shorter than the inverse 
Compton cooling distance, pair trajectories are randomized by the IGMF be- 
fore the pair energy is lost to the inverse Compton emission. In this case, the 
secondary emission is isotropic and its angular extent is determined only be 
the mean free path of the primary 7-rays. From Fig. [TU]one finds that in this 
case 6>ext = A 77 # b s /-D. ^ n case ^ ext does n °t depend on the IGMF. The 
condition Die < Rl gives the limiting IGMF strength measurable using 7-ray 
techniques 



eDic 



3 x 10" 



1 TeV 



G 



(148) 



Properties of the e xtended 7-ray em i ssion from IGM in the limit B > B max 
were considered by Aharonian et all (1994). The dependence of 6> GX t on A 77 
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leads to a dependence 6> ex t ~ E~ x / 2 of the source size on the energy of the 
cascade photons. This characteristic energy dependence is, in principle, ob- 
servable. Observation of this energy dependence would immediately imply a 
lower bound on the IGMF strength at the level given by Eq. (| 14811 . 

Searches for extended and/or time-delayed emission from the 77 pair cas- 
cade in the IGM were done with data from the Fermi telescope and from 
ground-based Cherenkov telescopes over the recent years. 

Aharonian et all ( 2001 ) searched for extended emission around the blazar 



Mrk 501 i n the energy band above 500 GeV using data from the HEGRA 
telescope. Aleksic et al ( 2010l ) have searched for extended emission around 
bright blazars Mrk 421 and Mrk 501 in the energy band above 300 GeV using 
data from the MAGIC telescope. No signal in excess of the instrument point- 
spread function (PSF) has been detected in both cases. This might impose 
some restrictions on IGMF in the range B ~ 10~ 13 G, provided that the 
cascade signal in the energy band E > 300 — 500 GeV is at a level higher than 
~ 1 — 10% of the prima r y sou rce signal. However, this is not necessary so, 
as shown by iTavlor et~al (|201ll ). Thus the HEGRA and MAGIC data do not 
constrain the IGMF. 

Neronov and Vovkl ( 2010l) ; Tavecchio et all ( 2011 ) have considered data from 



the Fermi telescope and deduced a lower bound on the IGMF at the level of 
B > 10~ 16 G from the non-observation of the cascade signal in the GeV band. 
This limit is independent of the correlation length of the IGMF as long as 
\b 3> -Die- If j to the contrary, As <C -Die, the limit becomes more stringent 

— 1/2 

with the decreasing Xb- The necessary magnetic field increases as X B , be- 
cause of the randomness of deflections of electron/positron trajectories, see 
Eq. (|145|) . This limit is shown by the light-blue hatched region in Fig. [T3] 
Fermi observations are more sensitive to the cascade emission because a num- 
ber of extragalactic VHE 7-ray sources like 1ES 0229+200, 1ES 0347-121 or 
1ES 1101-232 observed with the Fermi telescope have hard intrinsic spectra 
extending into TeV energy band. The intrinsic energy flux of these sources in 
the GeV band is much lower than that in the TeV band. Most of the intrinsic 
power of the sources in the TeV band is converted into cascade emission from 
the IGM, which is released in the GeV band. Thus, the expected energy flux of 
the cascade emission in the GeV band is much higher than the intrinsic flux. 
This makes the cascade emission more easily detectable. 

The limit found by iNeronov and Vovkl (|2010l ) and iTavecchio eTal (|201ll) 
was derived assuming that the cascade signal is suppressed by the large angular 
extension of the cascade source. This is not necessarily so. An alternative way 
to s uppress the cascade e mission is via a suf ficiently large time delay of the sig- 
nal. Dermer et al ( 2011 ); Taylor et all ( 2011 ) used simultaneous observations of 
blazars in the GeV and TeV band by Fermi and ground-based Cherenkov tele- 
scopes to deduce the lower bound B > 10~ 18 — 10~ 17 G, under the assumption 
that the casc ade emission is supp ressed by the long time delay of the signal . 
The limit of iDermer et all (|201ll ) is lower than that of ITavlor et~al (l201lh 
because of the difference in the modeling of the cascade signal. Dermer et all 
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Fig. 12 Constraints on the IGMF from the non-observation of 7-ray cascade emission. 



(|2011l ) used a semi-analytic al model for the pai r production and inverse Comp- 
ton emission spectra, while rTavlor et al ( 2011 ) used Monte-Carlo simulations 
which take into account the detailed di fferential cro s s-sect ions for the two pro- 
cesses. The limit on IGMF derived by Tavlor et all ( 2011 ) is shown in Fig. [T2l 
by the light-blue shaded region. 

An important uncertainty in the lower bounds on IGMF derived from 7-ray 
data stems from the uncertainty of the the EBL measurements. Indeed, the 
overall power of the cascade source is equal to the fraction of initial 7-ray power 
of the primary source, absorbed in the IGM. This fraction is proportional to 
the suppression factor exp(— Z)/A 77 ), which, in turn, depends on the EBL 
density 7Jebl via A 77 , see Eq. (|140[) . The dependence of the power of cascade 
source on tiebl is exponential, so that even a moderate uncertainty of a factor 
2 in the EBL density induces an order-of-magnitude uncertainty in the cascade 
power. This uncertainty affects the prediction for the cascade flux in the GeV 



-18.5 -18.0 -17.5 -17.0 -16.5 -16.0 -15.5 

logio B [G] 



Fig. 13 Dependence of the lower bound for IGMF (at large correlatio n lengths) on the 
assum ed leve l of the EBL den sity, in units of the model predictions by [Franccschin i et"al 
fcOOgT) . From lVovk et all i2012T ). 



band, which is then compared to the data. lVovk et all (120121 ) have investigated 
the influence of this uncertainty on the lower bounds on IGMF and found 
that reduc ing to EBL density to th e leve l of 0.8 of the den s ity as sumed in the 
models of Franceschini et al ( 20081 ) and iDommguez et~al ( 2011 ) reduces the 
lower bound on the IGMF signifi cantly, see Fig. [T3] Similar conclusions have 
been reached by Arlen et"al ( 2012 ) who found that uncertainties of the intrinsic 
source spectra combined with the uncertainties of the EBL measurement might 
ev en wash out the lower bou nd if the EBL is significantly below the value found 

bv lFranceschini et all (|2008l) . 

No te that recent measurement of the EBL spectrum by HESS ( H E.S.S. Collaboration et all 



2012a|) rule s out EBL spectra with n ormalizations lower than that of lFranceschini et all 



( 2008 ) and Dommguez et all ( 2011 ). see Fig. [7] This implies a lower bound on 
IGMF which is somewhat stronger than 10~ 17 G, see Fig. ITUl 

To significantly suppress the cascade signal, the IGMF has to be present 
over a large fraction of t he line of sight toward an extragalactic VHE 7-ray 



source. 



Dolag et~al ( 2011 ) have found that this imposes a bound on the possible 



volume filling factors of the IGMF at the level of > 60%. 

Detection of extended 7-ray emission around extragalactic VHE 7-ray 
sources, with an energy-dependent morphology, would imply a measurement, 
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rather than just a bound on the IGMF. Such a measurement was claimed by 
Ando and Kusenko ( 2010h . based on data from the Fermi telescope. However, 



the excess 7-ray signal co ming; from outside the Fermi PSF in the analysis of 
Ando and Kusenkol (|2010h turned out to be not due to a real sig nal, but due to 
the im perfect modeling of the telescope's PSF, as was shown bv lNeronov et al 
f|201lh . 

An alternative method to detect the cascade emission around blazars is to 
look for the signature of the gamma-ray induced casc ades in the anisotropy 
spect rum of the Extragalactic Gamma- Ray Background ( Venters and Pavlidoul 



2012). 



A promising method for measuring an IGMF with strength close to the 
lower bounds derived from Fermi data is to look for an energy-dependent 
delayed emission following bright flares of TeV blazars. Indeed, from Eq. (|147|) 
one can see that emission from very small off-source angles 6> cx t ~ 1" might be 
detectable in this way. This range of the off-source angles can only be reached 
by imaging observations wit h 7-ray telescop e s whi ch have the PSF smaller 



than ~ 0.1° ~ 3 x 10 arcsec. iTakahashi et all (|2012r> have attempted a search 



for delayed cascade emission in Fermi, following a bright TeV flare of the blazar 
Mrk 501. Unfortunately, the flare flux and the spectrum of the flare in the TeV 
band were not high and hard enough to result in significant cascade emission 
which would be detectable on relatively short time sca les (day-to-mqnth) b y 
Fermi, even if the IGMF would be negligibly small, see iNeronov et a j (|2012h . 
It is clear that next strong TeV flare of a blazar with hard intrinsic spectrum 
extending into multi-TcV band, such as 1ES 0229+200, will provide strong 
constraints on IGMF or will allow us to measure it, if the IGMF level is close 
to the currently existing lower bound. 

However, the known hard spectrum TeV blazars are surprisingly devoid of 
flaring activity. One possible explanation for this might be that the production 
mechanism of VHE 7-rays in these sources is different from that operating in 
flaring sources. In particular, VHE 7-ray emission can be related to the emis- 
sion of UHECR from these sources. A non- variable VHE 7-ray flux can be pro- 
duced in interactions of a beam of pro tons with energies above 10 18 eV during 
their propagation through the IGM (jEssev et a lEssev and Kusenkol 



2012). 



Observations of secondary VHE 7-ray emission from UHECR interactions 
in the IGM is possible only if the UHECR induced cascade develops in an 
anisotropic way along the UHECR beam, and is not isotropized due to the 
deflections of the cascade electrons and positrons by an IGMF. This implies 
that the UHECR cascade scenario for the VHE 7-ray emission wo rks only if 
the IGMF strength is close to the lower bounds (jEssev et alll2011al) . 

The limits on IGMF derived from 7-ray observations stem from the non- 
observation of the cascade inverse Compton emission initiated by the absorp- 
tion of the VHE 7-rays in the IGM. A potential alternative possibility to 
suppress the cascade signal is to dissipate the absorbed 7-ray power not via 
inverse Compton scattering of CMB photons with e + e~ pairs, but through a 
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different channel. Such a possibility was considered by Broderick et all ( 2012 ). 
who considered the possible role of plasma instabilities on the geometry of 
the e + e~ beam in the IGM. In the absence of IGMF, electrons and positrons 
deposited along the primary VHE 7-ray beam, form a collimated beam of 
high-energy particles, with the density 



n-bc 



10" 



FvHE 



D 



10- 12 erg/(cm s)/ V600 Mpcy V 1Te V 



cm 



-3 



(149) 

where we have substituted the flux and the luminosity distance of 1ES 0229+200 
as the reference values for -Fvhe and D. This implies the mean distance be- 
tween the beam particles fbcam ~ "beam — 10 8 cm - 

The collimated electron-positron beam propagates through the ionized 
IGM with free charge density tt-igm = ^bP - 3 x 10~ 7 (f?b/0.04) cm -3 . 
Plasma instabilities due to collective interaction of beam particles can oc- 
cur if ^beam is shorter than the plasma skin depth £ s ki n — 27r/w p ~ 6 x 
10 9 (higm/3 x 10 



7 cm" 3 ) cm and/or the Debye length of the IGM plasma, 



5 c bye 



IGM 



1/2 



— ~10 6 

UJr, 
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-1/2 



?IGM 

10 4 K 
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cm 
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where u> p = Ane^niaM /m P , is the IGM p l asma frequency and Tigm is 
the temperature of the IGM. Broderick et all ( 2012 ) have considered differ- 
ent possible instability modes of the electron-positron beam and found that 
the growth rate of an " oblique" mode of Laingmur waves in the linear regime is 
faster than inverse Compton cooling of electrons and positrons. Thus, plasma 
instabilities developing in the beam can potentially reduce the power of inverse 
Compton emission and instead inject this power into the collective plasma 
motions and the electromagnetic field along the beam path. The c alculation 
of the growth rate of plasma instabilities by iBroderick et al ( 2012 ) was per- 
formed in linear approximation, i.e. not taking into account back reaction of 
the beam perturbations on the growth rate. However, b ack reaction is known to 
be important, especially for r elativistic particle beam s (|Kaplan and Tsvtovich 
1973t ISchlickeiser et"all2012l;lMiniati and Elvivll201^) . In particular, the effect 
of non- linear Landau damping suppresses the growth of plasma inst a bilitie s 
and stabilizes the beam . The calculations by Kaplan and Tsvtovichl (|l973l ); 



Miniati and Elvivl (|2012l ) and ISchlickeiser et all (|2012l ) agree on the fact that 



when nonlinear Landau damping is important, the growth of plasma insta- 
bilities is suppressed. However, they differ in the range of parameters of the 
beam for which nonlinear Landau damping is supposed to be important. Fur- 
ther investigations are required to asses the potential importance of plasma 
instabilities for the electron-positron beam. 

We also notice that development or suppression of plasma instabilities is 
highly sensitive to the angular and energy distribution of thze particles in the 
beam. The presence of cooled (rather than freshly injected) electrons in the 
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beam with energies much below TeV and/or the presence of tiny IGMF with 
the strength much below the limits discussed above, can destroy the narrow 
collimation of the beam and suppress beam instabilities. This suggests that 
plasma effects might potentially be important only at the initial moment of 
the onset of blazar activity, on the time scales about the inverse Compton 
cooling time, t ic ^ -Die - 10 6 (E e /1 TeV) -1 yr. 



5.6.3 Prospects for IGMF measurement with next- generation ^y-ray telescopes. 



If the IGMF strength is close to the lower bound derived from the GeV 7-ray 
data from the Fermi telescope, a measurement of the IGMF can be achieved 
using a combination of the Fermi telescope with one of the ground-based 
Cherenkov telescopes: HESS, MAGIC or VERITAS. Fields with the strength 
~ 1CP 17 G can be measured via the observation of the time delay of the cas- 
cade emission following an exceptionally bright and hard spectrum flares of 
TeV blazars, similar to the brightest flares of Mrk 421 and Mrk 501 in 1997 
and 2000. In the particular case of Mrk 501, the cascade emission in the GeV 
band might be detectable by Fermi only if the 1-10 TeV flux of the source is 
at the level of the h istorically brightest flare observed by HEGRA and CAT 
( Neronov et a 1 2012). Otherwise, the quiescent source flux level of the source 
in the GeV band is high enough to hide the cascade emission. 

The duty cycle of exceptionally bright flares during which the TeV flux of 
a blazar can grow by up to two orders of magnitude is still unknown. After 
an all-sky monitoring of the 7-ray sky by the Fermi telescope of four years, 
no simultaneous GeV- TeV band observations of any of the exceptional flares 
is reported. This suggests that either flares are rare, with typical duty cycles 
much longer than a decade, or that the TeV ban d flares are not nec essarily 
associated to the flaring activity in the GeV band ( Neronov et all 2012 ). If this 
is so, the flares can be observed only by telescopes working in the TeV band. 
However, the narrow fields of views of ground-based Cherenkov telescopes do 
not allow for efficient monitoring of a significant fraction of the sky, similar to 
that done by the Fermi telescope in the GeV band. 

An all-sky monitoring in the TeV band is important for an efficient identi- 
fication of flaring episodes. The possibility of w ide FoV telescope s in the TeV 
band is now demonstra ted by the MILAG RO (lAbdo et alll2007l). Ti bet- AS-7 
(jAmenomori et all l201ol ) and ARGO-A BJ dD'Ettorre Piazzolil 12013^ arrays 
Next generation arrays, such as HAWC ( DeYoung and HAWC CoUaborationl 
(2012), http://hawc.physics.wisc.edu/) and LHAASO (http://english. 
ihep . cas . cn/ ic/ ip/LHAAS0/ ) will have an order-of magnitude better sensi- 
tivity and somewhat lower energy threshold than MILAGRO, Tibet- AS7 and 
ARGO-ABJ. 

A wide field-of-view monitoring of the sky will be also possible with the 
next-ge neration Cherenko v telescopes, like the Cherenkov Telescope Array 
(CTA) (jDubus et all 120121 ) . The monitoring will be possible in the "sky sur- 
vey" mode in which some ^ 20 — 30 sub-arrays of Cherenkov telescopes with 
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individual FoVs of 7° — 10° will be pointing in slightly different directions, 
covering regions of ten(s) of degrees on the sky. 

This will allow an efficient monitoring of blazar activity, rather than just 
source detection. If the Fermi telescope will still be in orbit at the time of 
full operation of HAWC and/or LHAASO, detection and detailed GeV-TeV 
monitoring of exceptional flares will strongly enhance the chance to observe 
the delayed cascade emission from the IGM and thereby to measure an IGMF 
in the range of (1CT 17 - 10" 16 )G. 

The time delay of the cascade emission scales as ( Neronov and Semikod 
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where z is source rcdshift and r is the optical depth for the primary absorbed 
7-rays responsible for the secondary cascade photons at the energy E. 

A "natural" upper limit for the energy of 7-rays produced by blazars is 
the energy at which the VHE 7-rays can not escape from the blazar host 
galaxy because of pair production on the interstellar radiation field (ISRF) 
inside the host galaxy. Taking the nearby radio galaxy M 87 as a proto- 
type host galaxy of TeV blazars, we estimate that photons with energies 
above ~ 30 TeV are efficiently absorbed while propagating through the ISRF 
(jNeronov and Aharonianl 120071 ). This implies that VHE 7-ray blazars should 
have a high-energy cut-off in their intrinsic spectra at ~ 30 TeV energy. Taking 
into account the relation between the primary and cascade photon energies, 



E = 1 



E~. 



30 TeV 



TcV, 



(152) 



we conclude that the delayed cascade emission is observable up to the 0.1-1 TeV 
energy band, by ground-based Cherenkov telescopes. Taking a time scale of 
~ 0.3 — 1 yr as a reasonable span for an observation campaign following an 
exceptionally bright flare, we find that magnetic fields with strength up to 
3 x 10~ 16 G are measurable in this way. Of course, it might be challenging to 
distinguish the delayed cascade emission signal in the 0.1-1 TeV band from the 
intrinsic slowly decreasing flaring emission from the primary source in the same 
energy band. Only observations of multiple flares (from several sources) with 
characteristic energy dependent decay time given by Eq. (|151|) will indicate 
that the delayed emission is due to the cascade, rather than being intrinsic 
to the sources. The range of IGMF measurable via observations of delayed 
emission is shown by the orange shading in Fig. [TU 

If IGMF is significantly stronger than 10 _16 G, the time delay of the cascade 
emission in any energy band is too large to be directly measurable by 7-ray 
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telescopes within a reasonable observation time span. However, in this case the 
angular extent of the cascade emission becomes sufficiently large to be mea- 
surable by the 7-ray telescopes. For a source with a jet a ligned closely with the 
line o f sight, the angular size of the extended source is ( Neronov and Semikozl 
20091) 
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Measurements of moderately extended (degree-scale) cascade emission at the 
highest energies E ~ 1 TeV with CTA and/or HAWC and LHAASO can 
provide a measurement of IGMF, if its strength is below about 10 _13 G (for 
large correlation lengths As). The range of IGMF parameters accessible via 
observations of extended emission around blazars in the 0.1-1 TeV band is 
shown as the lower orange hatched region in Fig. [T?J 

As mentioned above, still stronger magnetic fields are not directly measur- 
able using 7-ray tcchiqucs, because the properties of the extended emission 
from the IGM become independent of the IGMF. However, the presence of 
IGMF stronger than ~ 10~ 12 G in the IGM can still be establi shed from the 7- 
ray data, if the extended emission around blazars is detected ( Aharonian et all 
I994I ; iNeronov and Semikozll2009l) . 



5.7 Constraints from initial seed fields for galactic dynamos 



Cosmological magnetic fields might, in principle, play an important role in 
the cosmic magnetogenesis, because they can be the "seed" fields necessary 



for the action of dynamos in the galaxies ( 


Parkerll955: Ruzmaikin et all|l988 


Kronbers 1994; Beck et all 19961: iKulsrudl 


19991; iGrasso and Rubinstein! 12001 


Widrowll2002l: Brandenburg and Subramanian 2005). Followinethe Qualitative 



arguments used in Section 0] for the description of cosmological evolution of 
magnetic fields, one can estimate the growth time scale of Galactic dynamos 
igai based on the known velocity scale v ~ 10 7 cm/s (v ~ 10 6 cm/s for turbulent 
motions on the scale A < lOOpc in the Galaxy, and v ~ 10 7 cm/s for the large 
scale motions on distances A ~ (1 — 10)kpc). The "eddy processing" time is 
estimated as 

i ga] ~ - ~ 10 8 ( -^—) ( / ; _1 yr . (154) 

This eddy processing time is close to the typical estimates of the growth rate 
(e-folding time ) of the Galactic dynamo in the range of t ~ 10 8 — 10 9 yrs 
(|Kulsrudlll999l) . This means that within some 20-30 e-foldings one would am- 
plify the initially existing field by a factor of exp(to/r) ~ 10 9 — 10 13 . The 
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Fig. 14 The range of IGMF parameters accessible by next-generation 7-ray and UHECR 
telescopes. The orange shaded region shows the range of parameters for which delayed 
emission following bright flares is detectable at the energies below ~ 1 TeV. The lower orange 
hatched region shows the range of parameters for which extended emission is detectable 
below ~ 1 TeV. The upper orange hatched region shows the range of parameters which can 
be constrained by future UHECR telescopes. 



observed amplitude of galactic magnetic fields of about 10 fj,G can, therefore, 
be explained if the pre-existing seed magnetic field has the strength of at least 
(10 -21 — 10 _19 )G, taking into account additional amplification by about 3 
orders of magnitude via compression during the collapse of primordial pertur- 
bations to galaxies, see Eq (|113[) . Also galaxies at high redshift, z ~ 1 — 2 
have similar magnetic fields. Since the age of the Universe at z ~ 2 is more 
than a factor of 2 less than to, the dynamo amplification in these galaxies is 
only about 10-15 e-foldings requiring correspondingly stronger seed fields of 
about 10~ 15 G. The fact that galactic fields at redshifts z = 2 are of the same 
order as this at z — hints to the fact that galactic dynamo amplification is 
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saturated already at z ~ 2 and the amplitude of these fields is not determined 
by the strength of the seed fields. 

If the mechanisms of the action of Galactic dynamos would be well con- 
strained, we could use observations of the structure and strength of magnetic 
fields in different types of galaxies (and galaxy clusters) and at different red- 
shift to gain information about the properties of initial seed fields for the 
dynamo action. However, the efficiency of dynamos in galaxies and galaxy 
clusters are uncertain, so that the existing estimates of the strength and spa- 
tial structure of the seed magnetic fields needed for the Galactic dynamos 
differ in a very broad range. 

The seed fields for the Galactic dynamos might, in fact, originate not from 
the pre-existing primordial fields, but occur locally in the forming galaxies. 
Possible mechanisms for the generation of m agnetic fields during the grav- 
itational collapse ar e the Weibel inst ability ( Schlickeiser and Shukla 2003t 
Medvedev et a 3l2006l) or battery effects (jSubramanian et a lll994t Kulsrud et all 



1997t iGnedin et alll2000h . These mechanisms do not result in significant mag- 



netic fields spread in the IGM. Therefore, in such scenarii, the seeding and 
dynamo amplification of magnetic fields in galaxies is not at all related to the 
properties of IGMF. 

Besides, even if the seed fields for the Galactic dynamos stem form pre- 
existing magnetic fields produced in the Early Universe, the value of the seed 
fields might not be clos e ly rel ated to the primordial field strength. Rvu et al 



(|2f)f)8h : ISchleicher eta] (l2010h argue that the small-scale turbulent dynamo 



operating in the IGM during the collapse of proto-galaxies amplifies initial 
seed magnetic fields up to equipartition with the turbule nt energy on tim e 
scales much shorter than the dynamic al time scale, see also lBeck et a 3 (|l994h : 
Brandenburg and Subramanianl ( 2005 ). This implies that the initial conditions 



for Galactic dynamos are largely independent of the pre-existing seed field 
strength. 

Thus, although cosmological fields produced before the onset of structure 
formation might serve as the initial seed fields for Galactic dynamos, no sen- 
sible constraints on the cosmological fields can be derived from the observed 
properties of Galactic magnetic fields. 



5.8 IGMF from Galactic winds 



Physical processes in the Early Universe are not the only possibility to gen- 
erate magnetic fields in the IGM, in particular in the voids of the large scale 
structure. This means that potential measurement of the IGMF do not auto- 
matically imply the discovery of relic primordial fields. An alternative possi- 
bility is that the magnetic fields can spread through the IGM at late times by 
large-scale outflows from magnetized galaxies. 

The magnetic fie l d might spread into the voids by the outflows from radio 
galax ies (Regs 1987 ; Daly and Loeb 1990t Ensslin et al 1997 ; Kronberg et al 
20011 ). These outflows in the form of jets delivering high-energy particles into 
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Fig. 15 Range of IGMF parameters expected in scenarii where IGMF is spread by super- 
nova driven Galactic winds (light-blue shading and hatching) and/or cosmic rays (orange- 
shading). 



large scale radio lobes are a characteristic feature of radio galaxies (lUrrv and Padovani 



1995). The observed size of the jets reach Mpc scales. iFurlanetto and Loebl 
discuss the possibility that magnetic fields spread by quasar outflows 
and find that by redshift z ~ 3 (i.e. at the maximum of the quasar activity) 
some 5%-20% of the volume of the IGM might be "polluted" by fields from 
quasar outflows. The strength of the field in the polluted regions would be at 
the level of ~ 10% of the thermal energy density of the IGM with temperature 
T ~ 10 4 K, which implies B ~ 10~ 9 G. The correlation length of these fields 
would be of the order of the size of the radio lobes, i.e. in the Mpc range. 
This range of parameters is shown as olive-shaded region in Fig. 1151 In gen- 
eral, the conclusions on the strength of the fields spread by radio galaxies and 
quasars strongly depend on assumptions about the generation of the outflows 
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by quasars: the duration of the period of activity of the quasar central engine, 
its energetics, magnetization etc. 



Alternatively, magnetic f ields might be spread into the IGM by the galactic 
winds driven by supernovae (jKronberg et a Supernova dri ven winds are 

most probably responsib l e for the metal enri c hment of the IGM ( Aguirre et al 



2001 



Theuns et all 12003 iBertone et all 120051) . iBertone et all (I2006D used a set 



of prescriptions for the magnetic field content of supernova driven winds, to 
estimate the volume filling factor of the resul ting IGMF, based o n numerical 
simulations of the wind spreading in the IGM (jBertone et alll2005l) . Their con- 
clusion is that most of the IGM might be polluted with galactic wind fields 
with strengths ranging from 10~ 12 G to 10 _8 G. The correlation lengths of the 
wind-spread fields are expected to be about the size of their footprints in the 
galaxies, i.e. in the A_b ~ 1 — 10 kpc. This estimates strongl y depend not 
only on the "prescriptions" for the magnetization of the winds ([Bertone et al 
20061 ). but also on the "prescriptions" for galactic winds ( Bertone et all 120051 ) 
which are used in numerical simulations. In particular, varying model param- 
eters, opposite results can be obtained: fr om relatively strong IGMF filling 
the voids with high volume filling factor ( Bertone et all [20061 ) down to the 
small fra ction of the volume (only about 10~ 2 ) occupied by the wind blown 
bubbles ( Bertone et a 1 2005). The range of parameters B, Xb of IGMF spread 
by galactic winds is shown in Fig. [15] by a light blu e shaded/hatched re gion. 
Donnert et al ( 20091 ) further developed the model of IBertone et al (120061) and 
performed simulations aimed at testing the hypothesis of the origin of mag- 
netic fields in galaxy clusters from galactic outflows. The results of the simu- 
lations show that the observed /iG scale fields in the cluster cores can be fully 
explained by spreading of galactic winds. 



Mini ati and "Bell (|201ll ) consider the possibility to generate IGMF by cos- 
mic rays escaping from galaxies during the period of rcionization at z ~ 10. 
In absence of strong primordial IGMF, cosmic rays generated by the first su- 
pernovae at the onset of star formation freely stream out of the galaxies and 
produce an current jcr which is compensated by a return current of low- 
energy plasma from the IGM. Spatial variations of Ohmic resistivity of the 
IGM, r\ ~ T IGI ^ caused by the variations of the temperature of the IGM, Tigm , 
provide a source term in the magnetic induction equation, dH/dt ~ (V?y) Ajcr 
wherever the gradient of rj is misaligned with Jcr- This source term might lead 
to magnetic field generation in the IGM which stops as soon as the IGM gets 
heated by the UV radiation produced by the star formation activity. Soon 
after rcionization, the resistivity drops and the source term in the induction 
equation disappears. This mechanism can lead to IGMF with strengths in the 
range ~ 10~ 18 — 10~ 16 G and with comoving correlation lengths in the range 
of typical fluctuations of the IGM temperature, \b ~ 10 — 100 kpc. We show 
this range of IGMF parameters by the orange shaded region in Fig. fTS) 
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6 Constraints on the observationally testable mechanisms of 
cosmological magnetogenesis 

If the magnetic fields present in the IGM are of cosmological origin, under- 
standing of the evolution of the fields from the moment of magnetogenesis 
until the present (see Section H]) can be used to obtain information about the 
physical processes in the Early Universe, based on astronomical observations 
/ limits of the IGMF at z = 0. This provides an attractive possibility of a 
new type of cosmological probe, potentially sensitive to the epochs preceding 
recombination (CMB decoupling) and Big Bang Nucleosynthesis. 

Processes leading to the generation of magnetic fields in the Early Uni- 
verse typically produce magnetic fields on short distance scales (Section [3J . 
Turbulent decay and damping of the fields during subsequent evolution lead 
to the decrease of the field strength and increase of the correlation scale, where 
most of the magnetic field power is located. As a result, relic magnetic field 
which might be present in the IGM today may be relatively weak with short 
correlation scale. 

Observational constraints on the IGMF restrict its strength and correlation 
length (see Section [3]). If the strength of the relic field expected in a cosmo- 
logical magnetogenesis scenario is higher than of the observationally allowed 
IGMF (B, As) plane for a given correlation length, the scenario is not realized 
in the Universe. If it is lower than the lower limit on B at the given correlation 
length, the dominant magnetic field at this scale must come from some other 
mechanism and the model can not be directly tested by the observations. The 
most interesting scenarios from the observational point of view are those which 
predict relic magnetic field parameters consistent with observational bounds 
on the IGMF. 

In this section we review the different models of production of magnetic 
fields in the Early Universe proposed in Section[3]in view of their compatibility 
with observational constraints on IGMF. 

6.1 Inflation 

As it is discussed in Section 13.1.11 the power spectra of the fields generated 
during inflation by coupling the electromagnetic field either to curvature or to 
the inflaton are expected to be blue, with most of the power concentrated on 
short scales. The correlation length of the fields at the end of inflation is of 
the order of the size of the cosmological horizon. This means that the power 
of the fields generated by inflation happening at very high energy scales (in 
principle, up to the Planck scale) is initially concentrated at very short scales. 
This power is subject to fast dissipation, so that the final strength of the relic 
magnetic field at z = is extremely small. 

As discussed in Section [3J the fields generated during inflation have typi- 
cally a blue spectrum peaking at the wavenumbers k* ~ -ffmf where H ln f is the 
expansion rate. The energy density of magnetic field is pb ~ (H- ln {/Mp) p, so 
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Fig. 16 Obscrvationally testable region for inflation generated magnetic fields in scenarios 
with inflaton coupling to the electromagnetic field. The shaded region shows the range of 
magnetic fields and comoving correlation lengths for which are excluded by observational 
bounds on IGMF. Arrows show evolutionary tracks of the field correlation length and field 
strength (see Section |4{. Dashed lines show the power spectrum of the maximally possible 
field generated at inflation occurring at different energy scales. 



that the magnetic field is largely sub-equipartition, unless H m i ~ Mp. The 
generically expected relation between ~ yfpB and A* ~ fc^ 1 is B* ~ A" 1 . 
This range of expected initial values of magnetic field parameters is shown 
by a thick black solid line in Fig. Q1)J The spectrum of the field is generi- 
cally expected to be blue, with the slope n s = 1 in the most optimistic case. 
As an example, we show by the dashed line in Fig. \W\ the initial field spec- 
trum for the case of inflation happening at the Grand Unification energy scale 
E in{ ~ 10 15 eV. 

Evolution of the field strength and correlation length from the inflation till 
the present day is governed by the turbulent decay and proceeds along one 
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of the tracks shown in Fig. [31 The choice between the three possible tracks 
depends on the helicity of the field and on the type of turbulence governing 
the evolution. Evolution of the fields with initial parameters lying on the thick 
solid line in Fig. [1^] and a blue spectrum could not produce a final field con- 
figuration in the observable range (unshaded area in Fig. even in the most 
optimistic case of evolution of maximally helical field via inverse cascade. The 
only hypothetical possibility to have an observable field strength today is to 
generate a field with the scale invariant spectrum n s ~ —3. In this case, the 
possible field parameters B, Xb today satisfy the relation (p?2"f with the field 
strength reaching ~ 10~ 9 G for the field generated at the Grand Unification 
scale, see Fig. [Tl>] However, as it is discussed in Section [3J up to now there 
have been no self-consistent models for the generation of the scale-invariant 
magnetic field at Inflation. 

Another possibility for an observable field generated at inflation is to as- 
sume that the field is significantly amplified e.g. by dynamo action during 
reheating in order to become observable in the present Universe. In this case 
the resulting field energy density could increase up to the equipartition with 
the plasma energy density ps ~ p- If this is the case, the initial field strength 
might be high enough to result in an observable field in the present day Uni- 
verse. For this to happen, the initial field parameters -B*, A* should be in the 
blue-shaded region in Fig. [T|)] The denser-shaded region shows the allowed 
range of parameters of the non-helical fields, while the light-shaded region 
shows the range of possible parameters of the helical fields. 

The lowest possibl e energy scale for infla tion is just above the Electroweak 
energy scale, see e.g. German et all ( 2001 ). The correlation length of mag- 
netic fields which might be generated at this low-energy inflation can reach 
A ~ 0.01 pc (see Section l"3.1.1[) . The maximal possible strength of the initial 
magnetic field, B(k) ~ yj dps /d log k < ^fp 1 from inflation for a typical spec- 
trum, n s = 1 for lowest possible scale E- m i ~ 200GeV is shown in Fig. [16] as 
dashed line. 

The decay of these magnetic fields in the course of cosmological evolution 
can proceed in three possible ways. First, for the case of incompressible turbu- 
lence, no amplification of the field at large scales occurs, the decay just removes 
power from the scales smaller than the scale of the largest processed eddies. 
The evolutionary track of the field then follows the initial power spectrum 
(dashed line in Fig. [T5] up to the point where the time scale of the processing 
of the eddies is comparable to the Hubble time. The locus of the present-day 
strengths and correlation lengths of magnetic fields is along the blue thick line 
in Fig. OH 

Alternatively, the field excites compressible turbulence right after genera- 
tion. In this case, the field strength at large scales is continuously driven into 
equipartition with fluid turbulence. This leads to a mild amplification of the 
field on large scales so that the evolutionary track of the field strength - cor- 
relation length follows the B ~ A - path discussed in SectionHJ This results 
in somewhat stronger relic field and somewhat larger final correlation length. 
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The locus of the possible present day parameters of the relic field is still along 
the inclined thick blue line in Fig. [TBJ but the evolutionary track is along the 

two less steep arrows (not parallel to the black dashed li nes). 

Finally, if the field is helical, as suggested e.g. by Durrer et all ( 2011 ). 
helicity conservation forces it to evolve along the significantly flatter path, 
B ~ A~ 0,5 . Again, evolution ends on the thick blue line fixed by the relation 
Xb/va = t - 

The explanation of these different evolution paths is given in Section 2] 
Tracing back the evolution of _B, Xb from the observationally allowed range 
at z = (shown by the thick blue line in the diagram I16[). we find the range 
initial values -B, Xb for observationally testable inflationary magnetogenesis 
models. This range is shown by the blue hatched region in Fig. [15] Non-helical 
magnetic fields from inflation fall in the observationally allowed range only if 
the energy scale of inflation is below E^f ~ 10 s GcV. This limit is removed 
for helical fields. 



6.2 The electroweak phase transition 



In a similar way one can delimit the range of magnetic field parameters (B, Xb) 
for observationally allowed models of electroweak magnetogenesis. This range 
can be found by back-tracing the evolutionary tracks of magnetic fields in the 
B, Xb plane as it is shown in Fig. 1171 The main difference with the inflationary 
scanarios is that many models discuss the production of helical fields for which 
the field strength and correlation length evolve along the line B ~ A^ 1 (see 
Section[4]) . This results in slower decay and potentially larger correlation length 
at z = 0, provided that the initial correlation length of the fields reaches a scale 
of the order of the cosmological horizon at the elect rowea k phase transition. 
The scenarios considered by Vachaspatl ( 1991 ) and Enavist and Olesenl 



([19931 ) result in fields outside the observationally testable region. Order-of 
magnitude estimates for the field strength which can be produced via such 
process is B ~ m^, with typical correlation length Xb < , where mw is 
the W-boson mass. These initial values of B, Xb are shown in the diagram of 
Fig.[T7J However, taking into account the order-of-magnitude nature of the es- 
timate of the boundaries of the regi on of observationa lly testable models once 
can still suggest that the scenario of lVachaspatil (|199ll ) might result in observ- 
able magnetic fields in the IGM, roughly at the level of the currently existing 
lower bound on the field strength and correlation length. This requires that 
the fields are produced immediately in the maximally helical configuration, 
otherwise, initial faster field decay toward the maximally helical configuration 
reduces the field strength to the range much below the level consistent with 

the observational lower bounds. 

The scenario discussed by Vachaspati ( 1991 ) and Enavist and Olesenl ( 1993 ) 
has a problem in that the initial distance scale on which the magnetic field is 
generated is shorter than the scales of both viscous and Ohmic dissipation in 
the primordial plasma soon after the epoch of Electroweak phase transition. 
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Fig. 17 Obscrvationally testable region for the models of non- helical (green colored region) 
and helical (green hatched region) magnetic field generation during the electrowcak phase 
transition. Arrows show the evolution paths of non-helical (compressible turbulence) and 
helical magnetic fields. Points with labels show various model predictions. For comparison, 
the blue thick line shows the range of relic magnetic field strengths and correlation length 
from testable inflationary magnetogenesis scenarios. 



Indeed, the mean free path of particles at the moment of the Electroweak 
phase transition is oA m f p ~ g _4 T~ 1 where g ~ 0.1 — 0.3 is the coupling con- 
stant. The effective viscosity in the gas of relativistic i s related to the m ean 



free path of the least coupled particle is v = A m f p /5, see IWeinbera (|197lf ) and 
Appendix El so that the viscous dissipation scale Xd — vjv is close to the par- 
ticle mean free path for the characteristic plasma/fluid velocities v. Assuming 
that the characteristic velocity scale is v ~ va, one finds that in the case of 
magnetic field with strength B ~ m^, v ~ va ~ 1- Furthermore, at these 
energies the magnetic diffusivity is <t _1 ~ g _2 T _1 , is of the same order of 
magnitude as the viscosity and the field is also damped by magnetic diffusion, 
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see Section |U However, on these small scales the MHD equations might be 
not applicable right after magnetic field generation. Instead, the evolution of 
the coupled magnetic field-plasma system on the distance scales comparable 
to particle mean free path has to be modeled using the Boltzmann equation 
for the particle distributions. 

The scenario of Diaz- Gil et al (|2008allbh considers topologically non-trivial 
configurations of the inflaton and Higgs fields at the end of the post-inflation 
preheating as a source of magnetic fields. This operates in the temperature 
range of the electroweak phase transition. The estim ates for the cha r acteris tic 
correlation scale and t he magneti c field strength in iDiaz-Gil et a 3 (l2008allbT ) 
are similar to those of IVachaspati (1991), with a correction for the fact that 
more elaborated numerical calculations (rather than qualitative arguments) 
result in magnetic field energy density which is two orders below equipartition 
with the energy density of the Universe at the temperature T* ~ 100 GeV. 
The expected initial parameters for the magnetic field in this scenario are also 
shown in Fig.[T71 With its much lower initial field strength, the fields generated 
in this model are clearly outside the observationally testable range. 

A larger initial correl a tion l engt h of magnetic field is expecte d in the mod- 
els of Grasso and Riottol (1998) and Ahonen and Enqvist (1998) in which the 
electroweak phase transition is assumed to be first order and the generation of 
magnetic field is related to the collision of bubbles of new phase which leads 
to macroscopic classical vortex-like configurations of th e gaug e fields forming 
at the Electroweak phase transition. iGrasso and Riottol (119981 ) found that the 
initial correlation length of magnetic field might be larger than A ~ T~*, while 
the field strength can still be somewhat below the equipartition. This larger 
correlation length (a numerical estimate A ~ lO 4 ! 1 ^ 1 in the case of a first 
order phase transition and A ~ 10T" 1 for the second order phase transition is 
given in the paper) relaxes the problem of immediate dissipation of the mag- 
netic field energy via plasma viscosity, which operates at a smaller scale. From 
Fig. 1 171 we se e that the strength and ini tial correlation length of magnetic field 
estimated by Grasso and Riotto (jl998l) are also outside the region of directly 
testable models. 



Vachaspatil ([20011 ) has considered the generati on of helical magnet i c field s 



during the electroweak phase transition, see also ISemikoz et a 3 (|2009L 120121) 
A generation of non-zero helicity is related to the changes of Chern-Simons 
number for the SU(2)xU(l) gauge field which accompanies the generation 
of non-zero baryon number in electroweak baryogenesis. In this setting the 
helicity of the magnetic field generated simultaneously with baryons is directly 
related to the baryon number, h ~ 10 2 ns. Taking into account that the field 
is generated with a characteristic scale L ~ (e 2 T*) _1 , one finds that the field 
strength is roughly in equipartition with the energy density of baryons, which 
is 77b/9 ra d, where 776 ~ 10 -8 ' 5 is the baryon number of the Universe. This means 
that the magnetic field strength at generation is by a factor of ^/t% ~ 10~ 4 
below equipartition at the moment of generation. The initial values of the 
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magne tic field strength and correlation length in the scenario of Vachaspatil 
(|200lh . shown in Fig. [T7]lie also outside the testable range. 



Sim ilarly to the scenario of Vachaspatil (1991) an d Enqvist and Olesen 
(|l993l ). fields generated in the scenario of IVachaspatil (| 2 lh are initially on 



the distance scales close to the viscous damping scale. A potentially signifi- 
cant difference between the two scenarios might be in the fact that the field 
strength in the model of I Vachaspatil (]200ll ) is strongly sub-equipartition. The 
Alfven velocity is, therefore strongly non-relativistic va ~ \/B ~ 10~ 4 . Com- 
paring the viscous damping scale A^ ~ vv ~ uA m f p ~ 10 VT~ X with the field 
correlation length Xb ~ lO 2 ! 1 " 1 one can find that if the velocity scale of the 
dissipating modes is v ~ va, the field energy is not directly damped by the 
viscosity. However, it can be damped my magnetic diffusion. 

A much larger initial c orrelation length of the magnet ic field A ~ 10 7 /T* 
is possible in the model of I Joyce and Shaposhnikovi (|1997| ). who consider the 
production of helical magnetic field due to an imbalance of right- and left- 
handed electrons and positrons. This imbalance leads to a non-zero chemical 
potential /i for right-handed positrons. The chemical potential results in a 
dynamo-like term in the induction equation for the magnetic field. This term 
is responsible for the field amplification at the scale A ~ /.t -1 ^> T^ 1 . The 
dynamics of the field generation and evolution in this model in the temperature 
range T < T* and at higher temper atures T < 80 TeV has recent l y bee n 
considered by Bovarskv et all ( 2012bl lal) and Dvornikov and Semiko j ( 2013 ). 
The presence of relatively strong hypermagnctic fields at the moment of the 
electroweak phase transitio n can modify the d ynamics of the phase transition 
itself making it first order ( Elmfors et all 1998 ) and can further modifying the 
magnetic induction equa tion describing the evolution of the magnetic field 
(|Semikoz and Vallell2008h . 



Bavm et al consider much larger bubbles of new phase in a first- 

order electroweak phase transition, with sizes reaching ~ (10~ 3 — 10~ 2 )£h. 
Developed MHD turbulence with largest eddies of the size of these bubbles 
leads to equipartition between the energy density of magnetic fields and the 
energy density of fluid motions. Taking into account that the bubble walls 
can reach mildly relativistic speeds, the energy density of fluid motions can, 
in fact, be comparable to the energy density of the Universe. The magnetic 
fields are, therefore, also amplified to an energy density comparable to the 
energy den sity of the Univers e. The initial magnetic field parameters in the 
scenario of Bavm et~al (1996), shown in Fig. [T7] are clearly well inside the 
observationally testable region. In fact, if a measurement of IGMF today with 
parameters (B, Xb) which, after back tracing to the Early Universe will imply 
initial magnetic field strength close to equipartition with the rest of matter in 
the Universe at scales close to the Electroweak horizon scale would provide a 
strong argument in favor of electroweak magnetogenesis and in favor of a first 
order electroweak phase transition. 
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Fig. 18 The observationally testable region for models of non-helical (orange colored re- 
gion) and helical (orange hatched region) magnetic fields generated during the QCD phase 
transition. Arrows show the evolution paths of non-helical and helical magnetic fields. Points 
with labels show various model predictions. 



6.3 QCD phase transition 

For the QCD phase transition, the range of possible initial correlation lengths 
extends to larger values up to the Hubble radius at the temperature T* ~ 
100 MeV. This somewhat broadens the range of possible final field strengths 
and correlation lengths. However, similar to the case of the electroweak phase 
transition, most of the models consider field generation at very short scales, 
which lead to initial field parameters outside the observationally testable range, 

see Fig. HU 

In the model of iForbes and Zhitnitskv (|2000h the field is generated at the 



domain walls of the bubbles of a first-order QCD phase transition. The ini- 
tial correlation length of the field is aA„ ~ T^ 1 and its strength is close to 
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cquipartition, B ~ T% . The domain walls subsequently coalesce increasing the 
correlation length and decreasing the field strength as B ~ A" 1 . Since the 
exact value of the final correlation length is not known, the estimates of the 
field strength and correlation length in this model can only be represented in 
the form of a line in the diagram of Fig. [T5] This line lies outside the region in 
which the magnetic field might leave a detectable imprint in the IGM today. 



Kisslinger ( 2003 ) has proposed that, contrary to the assumption of lForbes and Zhitnitskv 



1)20001 ) . the magnetic field strength for the large domain walls might not be 
suppressed by a factor A _1 T*. He assumes that a field strength B ~ right 
at the scale of cosmological horizon at the moment of QCD phase transition 
(no special motivation for this assumption is provided though, although am- 
plification of the field up to equ ipartition by turbulence, similar to the scenario 
discussed bv lBavm et a 1 (1996) in the context of Electroweak phase transition 
is, in principle, possible). This is clearly the most optimistic scenario which is 
logically possible. In this setting, the field can produce measurable effects on 
the polarization of the CMB and also induce a gravitational wave background. 
The model predictions f or the field strength and correlation length for the fields 
obtained by iKisslingerl ([20031 ) fall within the range of observationally testable 
models, see Fig. [TS] The effects of helicity on the evol ution of magnetic field 
in this scenario was discussed by Tevzadze et al ( 20121) . 



In the scenario of lBoeckel and Schaffner-Bielichl (|2012|) the magnetic field is 
generated together with the baryon asymmetry of the Universe in result of the 
period of a "little inflation" accompanying a first order QCD phase transition. 
The energy density of the magnetic field generated in this way can reach the 
equipartition with the baryon energy density, so that B can be estimated 
as B ~ \/2rjb(had- The characteristic correlation length is determined by the 
diffusion of baryons across the walls of bubbles of new phase, occurring at a 
first order QCD phase transition. It is estimated to be in the range of ~ 10 cm 
physical length which is about 10~ 5 £h,qcd at the moment of the QCD phase 
transition. These reference values of the field strength and correlation length 
are shown in Fig. 1181 The model prediction lies within the observationally 
testable region, if the field is helical. 



Cheng and Olintol ( 1994 ) also consider the generation of magnetic fields by 



the propagating domain walls separating the two phases. The magnetic field is 
produced by an electrical current arising as a result of uncompensated charges 
on different sides of the wall. The field strength produced via this mechanism 
is significantly below equipartition with either baryon or total energy density 
and falls below the observationally testable range, see Fig. [T5] Subsequent 
a mplification of th e fields by MHD instabilities and turbulence was considered 
by Sigl et all ( 1997 ) who found that the fields can, in principle be increased up 
to a strength approaching equipartiton with the thermal energy of the plasma. 
This moves the initial field parameters in the range of observationally testable 
models, see Fig. [T5] 
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7 Conclusions 



In this work we have discussed the generation and evolution of cosmological 
magnetic fields (IGFM) and we have discussed their possible observation. 

Inflation only generates magnetic fields if conformal invariance is broken. 
But also then it must either happen at relatively low energy, E- m t ~ 10 5 GeV or 
generate helical fields in order to yield relic fields which are observable today. 
Otherwise the correlation length is too short and the fields are damped by 
viscosity during the subsequent MHD evolution below the amplitude required 
by observational lower bounds. If the magnetic field generated during inflation 
is not amplified subsequently to near equipartition, we expect a very low initial 
magnetic field amplitude determined by ps/p — (Hi n f/Mp) 2 which is well 
below the observable range for all inflation scales. 

We have also considered magnetic field generation by cosmological phase 
transitions, in particular we have discussed the electroweak and the QCD tran- 
sition. There we have seen that if the transition is of first order and proceeds 
via bubble nucleation, the correlation length can be sufficiently large and the 
amplitude can be sufficiently high for the fields to be observable today. 

This is very interesting as standard model physics predicts a cross-over for 
both, the QCD and the electroweak transition. Therefore, relic magnetic fields 
from the electroweak phase transition would be a signal of physics beyond 
the standard model. A first order QCD phase tr ansition is possible only if 
the chemical potential of the neutrinos is large, see lSchwarz and Stukd(l2009T l. 
In this case, we expect strong hclicity in the leptonic sector which modifies 
the evolution of magnetic fields via a significant contribution from the elec- 
tr oweak anomaly. This is expected to induce magnetic field helicity as outlined 
bv lBovarskv et all (|2012ah . 

Observationally accessible magnetic fields could be generated in the stan- 
dard model electroweak phase transition if the correlation len gth of the field is 

much larger than the inverse temperature scale, as proposed bv lJovce and Shaooshnikov 
(|l997h . 

A summary of the existing observational constraints and different model 
predictions for the IGMF is shown in Fig. Q1JJ This figure allows us to as- 
sess possible future measurements of the IGMF. Unfortunately cosmologi- 
cally produced fields and field ejected by the galaxies span the same range 
of field strengths, from the existing lower bounds from 7-ray observations 
(~ 10~ 17 G at large distance scales) up to the upper bounds from radio obser- 
vations (~ I0~ 9 G at large distance scales). Measurement of the field strength 
alone cannot provide a clue on the origin of the IGMF. In order to distinguish 
between early (cosmological) and late (galaxy formation) origin of IGMF, a 
combined measurement of the field strength and the correlation length is nec- 
essary. The detection of fields with correlation length shorter than ~ kpc favors 
a cosmological origin of the IGMF. If the IGMF correlation length is larger 
than ~ kpc, measurement of its strength significantly below ~ I0 -9 G seems 
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Fig. 19 A summary of observational constraints and model predictions for IGMF. Grey 
shaded region shows the range of parameters excluded by observations and thcroctical ar- 
guments presented in Section \5\ Solid lines mark the locus of possible present day strength 
and correlation length of relic magnetic fields produced at the phase transition in the Early 
Universe, as discussed in Seciton[4] Green and blue hatched regions show possible ranges of 
IGMF spread by the outflows from galaxies fScction l5.8ll . 



to imply that the field is produced by galactic outflows, rather than by the 
processes in the Early Universe. 

It is, therefore, crucially important to develop observational methods which 
allow a measurement of not only field strength (or of a particular combination 
of field strength and correlation length) , but also separate measurement of the 
field strength and the correlation length. 

Another possibility to distinguish cosmological IGMF from the fields spread 
by galactic winds is to search for turbulence in the voids of LSS. Indeed, 
cosmologically produced IGMF decays by transferring its power to turbulent 
motions of the plasma in the Universe. The latest episode of the magnetic 
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field driven turbulence might have happened recently, after reionization of the 
Universe by the star formation activity. The strength and correlation length 
of the relic cosmological magnetic fields today is such that the time scale of 
turnover of the eddies of the size comparable to the correlation length is just 
about the Hubble time, \b/va = to (this is the relation which defines the 
thick colored line in Fig. IT9|) . Therefore, the IGM in the voids of the LSS 
might be turbulent today if cosmological magnetic fields are present. A search 
of turbulent plasma motions in voids may therefore provide an alternative way 
to the discovery of primordial magnetic fields. 
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Appendix 



A Viscosities and Reynolds numbers 

In this work the Reynolds number of a given scale which is inversely propositional to the 
viscosity plays an important role since it defines the time when the velocity field and with it 
th e magnetic field on t he given scale is damped into heat. We closely follow the treatment 
of lCaprini eTal <2009d ). 



A.l Kinematic viscosity 

The kinematic viscosity is given by 

v = -^-, (155) 
p + p 

where rj is the shear viscosity. The kinematic viscosity characterizes the diffusion of trans- 
verse momentum due to collisions, and is given roughly by the mean free path £ m f p of the 
particles. In this Appendix ^ m f p is the physical mean free path, while A m f p appearing in the 
main text is the comoving mean free path. The relation is s imply l m f D = oA m f p . 
A more precise expression for the shear viscosity is, sec IWeinbereNl97lft 



V= TE ^ 9 * T * Vfp so that v = ^? ' (156) 

The largest viscosity comes from the weakest interactions. However, non-interacting particles 
do not contribute to the viscosity. For this reason simple analytical approximations to the 
viscosity have unphysical jumps whenever a species decouples from the plasma. 

Estimates from kinetic theory show that the shear viscosity of highly relativistic parti- 
cles, T»m, behaves as (to leading-log accuracy): 

V = C 4] T3 _ 1 , (157) 
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where g is the appropriate coupling constant (depending on the temperature and the length 
scale at which one wants to compute the Reynolds number) and C is a numerical coefficient 
that can only be obtained from a detailed analysis. 

At temperatures larger than the electroweak phase transition, neutrino interactions are 
not supp ressed. The shear v iscosity is dominated by right handed lepton transport and is 
given by l lArnold et alll2000t) 

/5\ 3 /12\ 5 3/2 T 3 

r '"UJ (i5)2 {~) 9^ + 224(5 + 1/2) g '4l ogg ,-i ( 158 > 
where g' is the hypercharge coupling. This leads to 

v(T > 100 GcV) — . (159) 

After electroweak symmetry breaking, neutrino interactions arc suppressed by a factor 
(T/Mw) 4 . In th is regime, neutrinos have the longest mean free path and dominate the 
viscosity. We use lHeckler and Hoganl <1993ft 

tmtp ~ (3GIT 6 )- 1 , (160) 

leading to 

u{T < 100 GcV) « 4.9 x 10 s . (161) 

At temperatures smaller than 100 MeV, after the QCD phase transition, the remaining 
relativistic particles in the cosmic plasma are electron/positrons, neutrinos and photons 
and the neutrino mean free path increases to 

^nfp-ytGlT 5 )- 1 (162) 

such that 

w(T < 100 MeV) rj 1.6 x 10 9 ■ (163) 

The neutrino mean free path determines the viscosity until neutrinos decouple at T ~ 1.4 
MeV, after which photons take over. Below lMeV, when electrons and positrons annihilate 
and the remaining electrons become non-relativistic, the viscosity can be approximated by 

u(T < 0.5MeV) w ((ttii,)- 1 ~ 0.5 x 10" 22 GeV" 1 f"^"") ■ ( 164 ) 

After neutrino decoupling, the viscosity drops by about 30 orders of magnitude and 
the Reynolds number increases correspondingly. Therefore, all scales on which turbulence is 
maintained until T ~ lMeV will the remain turbulent until decoupling, T ~ leV. 

But even if turbulence is lost before neutrino decoupling, as long as the magnetic field 
survives, it will become turbulent again after T ~ lMeV and we expect equipartition between 
the magnetic field and the velocity field to be re-established. 

Summing up all the results we find 

21.6GCV- 1 (-%^) ifT>100GcV 



"CO 



4.9 x 10 8 GeV~ 1 (-%^)° if lOOGeV > T > 100 MeV 
1.6 x 10 9 GcV" 1 (-%^)° if lOOMeV > T > lMeV 
0.5 x 10- 22 GeV _1 (^f^) ^ if 0.5 MeV > T > leV . 



(165) 



The unphysical jumps come from regions where our approximations are invalid. Nevertheless, 
at neutrino decoupling viscosity is significantly reduced and turbulence resumes on the 
relevant scales. This is different after matter and radiation equality since then the Alfven 
speed decays and the coupling of the magnetic field to the velocity field soon becomes 
negligible. 

The evolution of v with temperature is plotted in Fig. 1201 
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Fig. 20 Evolution of the kinematic viscosity v(T) as a function of temperature for T > leV. 
The unphysical discontinuities and kinks come from our crude approximation. 



A. 2 Magnetic diffusivity 

Here we derive expressions for the magnetic diffusivity also called resistivity for rclativistic 
electrons in the cosmic plasma with te mperatures 1 MeV < T < 100 GeV. Again, we follow 
the treatment of lCaprini et a 3 d2009cT) . 

To determine the magnetic diffusivity, we derive an expression for the conductivity c(T), 
which is the inverse of the diffusivity. The Lorentz force acting on an electron is 

m e = eF^Uv . 

dr 

If we average this equation over a fluid element containing many electrons, the magnetic field 
term is sub-dominant. Even though the electrons are highly relativistic, the average fluid 
velocity is small. Furthermore 7 = 1/ Wl—v^ ~ T/m e is nearly constant and we may neglect 
the contribution d-y/dr from du l /dr = d( / yv l )/dT above. With dr = , y~ 1 dt = (m e /T)di, 
this yields the following equation for the mean velocity of the electron fluid: 

dv e _ 
— = - E. 
dt T 

If we denote the collision time for the electrons by t c , they can acquire velocities of the order 
v ~ |;Et c between successive collisions. Hence the current is 

e 2 n e 

J ~ en e v ~ t c — E = <tE 



so that the conductivity becomes 



t c : 

T 



We now derive an estimate for t c from Coulomb interactions. For a strong collision between 
the electron and another charged particle we need an impact parameter b such that e 2 /b > 
E e ~ T. Hence the cross section becomes at ~ 7rfe 2 ~ ire 4 /T 2 (this simple argumentation 
neglects the Coulomb logarithms which enhance the cross section by ln(l/a m i n ) where a m i n 
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is the minimal deflection angle (see [Landau and Lifschitz 199(J)). With v e = 1 the time 
between collisions is therefore t c = l/(atn e ) ~ T 2 /(ire 4 n e ) and 

v s . (166) 

Note that this result is independent of the electron density. This is physically sensible as n e 
enhances the current on the one hand but it reduces in the same way the collision time. 
With (1 1 66 II wc obtain for the magnetic diffusivity 

1 e 2 (T) 1CT 1 1CT 2 

- ~ — ~ . 167) 

a AT T T K ' 

The first value applies close to T ~ lOOGeV, where a = e 2 /4n ~ 0.1, while the second value 
corresponds to low energies, T ~ lMeV. For non-relativistic electrons we obtain the standard 
result for the conductivity by simply rep lacing T by the electron mass and multiplication 



)r tnc conductivity by simpl y rep lg 
by u 3 ~ (m e /T) 3 / 2 so that iSpitzerlll97Sh 

.2 m V 2 

(168) 



1 2 1/2 

~ T 3/2 



(J 



A. 3 Reynolds numbers and Prandl number 

The kinematic Reynolds number is given by 

H*(T) = ^ , (169) 

where v is the kinematic viscosity, tlx = ^/ k^P v (kx) / (27r 2 ) is the mean velocity which is 
roughly the velocity at the integral scale \k = 2n/ki(, and v = u/a, see Section 2| 
Correspondingly, the magnetic reynolds number is defined by 

R m (T) = ^L, (170) 

Inserting the resistivity from Eq. (11671) and the kinematic viscosity from Eqs. 1161> or 
11631 . assuming cquipartition so that va = an( i As = Xk i wc obtain for the Prandl 
number 

This number is larger than 1 for all temperatures 1 MeV< T <g 100 GeV where the deriva- 
tion applies. 

The non-linearities in the Euler and induction equation are stronger than the damping 
term whenever the Reynolds numbers are larger than unity. In this regime MHD turbulence 
develops. 



A. 4 The Prandl number at very high energy 

Finally let us consider the situation at very high temperature assuming that all particle 
interactions are given by the same coupling strength g 2 and all particles are relativistic and 
in thermal equilibrium. This approximation is roughly valid above the electroweak scale. 
(We neglect strong interactions in this picture.) The cross section then is of the order of 
cr c ~ g 4 T~ 2 and 

tc = Amfp = (ven)- 1 ~ -j- ~ v , T > lOOGeV . (172) 
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This qualitatively reproduces Eq. (11 59 II - For the conductivity we have with the same ap- 
proximations 

<7 = t c — T>100GeV. (173) 
T g l 

In this case the Prandl number becomes 



10, 



T > lOOGeV . 



(174) 



It will be important for our discussions that at the elcctroweak phase transition T ~ 
lOOGeV, both, the kinetic viscosity and the magnetic diffusivity are actually of the same 
order. At significantly lower temperatures, the magnetic diffusivity is always much smaller 
than the kinetic viscosity. This is due to the fact that the kinetic viscosity is governed 
by the most weakly interacting particles, the neutrinos while the conductivity is of course 
determined by the stronger electromagnetic interactions. 



B Maxwell's equation in curved space times 

We consider the 4- velocity « M with u^u^ = — 1 in an arbitrary curved spacetime and define 
the electric and magnetic fields as in Section [2] eqs.[3] 

£ M = F^u" , B M = e M!/7 F"T/2 , (175) 

such that 

F^v = u^Ey - UvE^ + e M „ 7 .B 7 . (176) 

We define the expansion rate Q, the shear cr^y, the vorticity ujy^ and the acceleration a M of 
the 4-velocity by 

6 = , (Tftv = i (u^-v + Uv-fx - ^Oh^S) , (177) 

w M „ = - {Ufi-u — , a M = u"u M; „ . (178) 

Here h M „ = g M „ + u^Uy is the projector to the tangent space normal to u. 

In terms of th ese qu antities the homogeneous Maxwell equations, Fr^y^ = become 
(see iBarrow et al l l2007ft ) 

„y + + ~6h,tv) B v - e^a" E a (179) 

B" ;v = epvcU^El*. (180) 

The 3-component e-tensor is given in terms of the totally antisymmetric tensor Tja^va by 
= uPr/p^vc see Section [2] 
Introducing the 4-current j M , the charge density p e = —u^j^ and the 3-current J M = 
ju — peU/_i we obtain for the inhomogeneous Maxwell equations, F^ v ' ;u = j^, 

E v . v = p a - e^ a w va B>* , (181) 

-h^u a E v . a + e^B^ = - (a^ + + ^Oh^ E v - e^ va a v B a + J" . (182) 

The energy momentum tensor of the Maxwell field in terms of E B and u is 

J-IIV — -fjiaJ' V — -r a /3^ 9fiV 

= ^(E 2 + B 2 ) Uf ,Uy + ^{E 2 + B 2 )h^ - E^Ey - B^By + P^Uy + u^Py ,(183) 
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where P M = ^^ a pE a B@ is the Poynting vector, i.e. the energy flux seen by an observer with 

4-vclocity u. The energy density is p< cm ) = T^u^u" = (E 2 + B 2 )/2 and T^ cm) M = as 
we expect it from the electromagnetic energy momentum tensor. 

In a Fricdmann universe, for a co- moving observer with u = a~ 1 dt, these equations 
simplify considerably. Since uj m „ = cr M „ = a M = and 6 = 3a/a 2 = 3H, we obtain with 
= o(0, B) and = o(0, E). 

d t (a 2 B) + a 2 V AE = (184) 

V-B = 0, (185) 

V ■ E = ap e , (186) 

-d t (a 2 E) + a 2 V A B = a 3 J . (187) 

Note that B and E scale like 1/a 2 . In term if the re-scaled quantities a 2 B, a 2 E and a 3 J, 
the Maxwell equations assume the same form as in Minkowski space, the expansion factor 
can be 'scaled out'. 
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